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The Raman spectrum of benzene has been photographed using a fine slit and a high dis- 
persion glass Littrow spectrograph. The faint companion at 984 cm™ is found to possess an 
intensity ratio of 6: 100 when compared with the principal Raman line at 992 cm™. The 
view that this satellite arises from a benzene molecule in which one of the carbon atoms is 
replaced by the heavier isotope of atomic weight 13 is supported by this result. The distribu- 
tion of intensity in the 3048 and 3062 Raman lines has been determined and the contours of 
the lines are drawn separately. The peaks of these two lines are found to be in the ratio of 
1: 1.35. The relatively sharp nature of 3062 and the diffuseness of 3048 are correlated with 
the facts that the former is a well-polarized Raman line whereas the latter is a fully depolarized 


line. 





INTRODUCTION 


HE author is engaged in a systematic study 
of the fine structure exhibited by the vari- 
ous Raman lines under high dispersion and some 
of the results obtained with benzene are of special 
interest as they confirm the recent findings of 
Cheng, Hsueh and Ta- You Wu.! The results given 
in this paper relate to the 992 and 3060 Raman 
lines of benzene and are obtained independently 
by the author during the course of a more exten- 
sive investigation. They are now reported here in 
view of the publication referred to above. 


EXPERIMENTAL 


A quartz mercury arc is kept close on the top 
of a Raman tube cooled by circulating water. 
Polished aluminum reflectors are used to increase 
the illumination. A Littrow glass spectrograph 
with a dispersion of about 10A/mm in the \4358 
region has been used to photograph the scattered 
spectrum. By using a fine slit (0.03 mm) the 


(1938) Hsueh and Ta-You Wu, J. Chem. Phys. 6, 8 
90). 


Raman spectrum of benzene has been photo- 
graphed in which the component 984 is clearly re- 
vealed. In the same picture 3062 and its com- 
panion line 3048 are also recorded very strongly. 
A set of intensity marks are then put on the same 
plate by the method of varying slit widths. The 
lines at 984 and 992 and the lines at 3048 and 
3062 along with the intensity marks in the corre- 
sponding regions are run through a Moll micro- 
photometer and the relative intensities are de- 
duced in the usual manner with the help of 
density —log intensity curves. 


RESULTS AND DISCUSSION 


Figures relating to 992 and its component are 
not reproduced as they have already been given 
by a number of earlier investigators. The ratio of 
intensities of 984 and 992 lines, obtained as the 
average of three measurements comes out in this 
investigation as 6:100. This is in excellent 
agreement with the recently reported value of 
Cheng, Hsueh and Ta-You Wu.! No other com- 
ponents are recorded alongside this line. 
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Various investigators? have reported results 
on a study of the components of the 992 cm™ line 
in benzene. Of these authors Gerlach was the 
first to suggest that the faint companion at 984 
may be due to the principal oscillation of the 
benzene molecule in which one of the carbon 
atoms is replaced by its heavier isotope of atomic 
weight:13. Grassmann and Weiler, Bhagavantam, 
Specchia and Scandurra and Pattabhiramayya 
are among those who held that this view was 
correct. Subsequently Ananthakrishnan  ex- 
pressed the view that the origin of this companion 

2 J. Weiler, Zeits. f. Physik 69, 586 (1931). L. F. Howlett, 
Can. J. Research 5, 572 (1931). P. L. Mesnage, J. de phys. 
2, 403 (1931). W. Gerlach, Ber. d. Bayr. Akad. d. Wiss. 1, 
39 (1932). L. and E. Bloch, Comptes rendus 196, 1787 
(1933). P. Grassmann and J. Weiler, Zeits. f. Physik 86, 
321 (1933). O. Specchia and G. Scandurra, Nuovo Cimento 
12, 129 (1935). S. Bhagavantam, Proc. Ind. Acad. Sci. 2A, 
86 (1935). R. Ananthakrishnan, ibid. 3A, 52 (1936). 


P. Pattabhiramayya, Thesis submitted to the Andhra 
University (1937). 
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line may not be isotopic. The fact that the recent 
quantitative study of Cheng, Hsueh and Ta-You 
Wu already referred to and the author’s work 
reported here both lead to an intensity ratio of 
6 : 100 for the companion line in relation to the 
main line lends strong support to the view that 
the origin of 984 is isotopic. 

The continuous curve in Fig. 1 is a graph show- 
ing the intensity distribution in the 3048-3062 
lines. The broken lines represent the separate in- 
tensities of 3048 and 3062 which together by par- 
tial overlapping go to make up the observed in- 
tensity curve. These are drawn on the assumption 
that the lines are symmetrically broadened and 
it is easily seen that the line at 3062 is very sharp 
while that at 3048 is broad and diffuse extending 
on either side to about 25 cm~. These intensity 
characters are in perfect agreement with their 
observed depolarization factors. Veerabhadra 
Rao’ and Simons‘ have respectively given 0.88 
and 0.89 for the depolarization factor of 3048 and 
the line should therefore be expected to exhibit 
a strong wing on either side. On the other hand 
they have given 0.30 and 0.34, respectively, for 
the 3062 line and it is accordingly comparatively 
much sharper. It may easily be seen from the 
graph that the ratio of the peak intensities of 
these lines is 1 : 1.35 which agrees with the 
figures given by Veerabhadra Rao, namely 
1 : 1.28. The figures given in this paper are likely 
to be slightly in error on account of the fact that 
illumination from the top is adopted and the in- 
tensities of two lines whose states of polarization 
are different are compared.® Recent work of 
Angus and others ® has shown that both these 
lines are distinct and they represent two separate 
fundamental frequencies of the benzene molecule. 

The author is highly thankful to Mr. S. Bhaga- 
vantam for his kind interest in the work. He is 
also thankful to the Syndicate of the Andhra 
University for granting him a research fellowship. 


3 Veerabhadra Rao, Zeits. f. Physik 97, 154 (1935). 
4Simons, Comm. Soc. Scien. Fenn. 6, 13 (1932). 

5 A. Veerabhadra Rao, Proc. Ind. Acad. Sci. 7, (1938). 
6 R, Angus and others, J. Chem. Soc. 912-987 (1936). 
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Exchange Reactions with Deuterium 
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The exchange reaction between hydrogen chloride and deuterium, initiated photochemically, 
is found to be a chain reaction with quantum yields ranging from 2 to 19 under the experimental 
conditions employed. The rate determining step is the reaction Cl+D.—DCI+D. The pro- 
portions at the photostationary state are shifted from those at thermal equilibrium in a di- 
rection which indicates a somewhat stronger absorption by HCl than by DCI in the region 


2000—2150A. 





HEN mixtures of hydrogen chloride and 

deuterium, or of deuterium chloride and 
hydrogen, are exposed to light absorbed by the 
hydrogen chloride or deuterium chloride, an ex- 
change reaction is observed which proceeds 
rapidly at first, then approaches a photosta- 
tionary state (Fig. 1) at which the equilibrium 
proportions are near to, yet definitely different 
from, the proportions corresponding to thermal 
equilibrium at the temperature concerned. Using 
the apparatus previously described,' we have 
studied this exchange with respect to both the 
initial rate and the photostationary state. 


EXPERIMENTAL 


The hydrogen chloride-deuterium or deuterium 
chloride-hydrogen mixtures were exposed, in the 
quartz reaction vessel R (Fig. 1 of reference 1), 
to the light of a zinc spark; the effective radiation 
consisted of the four strong lines in the region 
2000-2140A. Quantum yields were obtained by 
actinometric comparison with the rate of de- 
composition of ammonia. Measurements of the 
absorption of hydrogen chloride and ammonia? 
show that for the system used (average light path 
~25 mm) absorption of the zinc lines at 2000- 
2140A may be taken as complete (>95 percent) 
at ammonia pressures of ca. 50 mm or over, while 
HCl absorbs even more strongly. 

The rate of ammonia decomposition was ob- 
served through the pressure change, with the as- 


1936)" Cross and P. A. Leighton, J. Chem. Phys. 4, 28 

*S. W. Leifson, Astrophys. J. 63, 73 (1926). H. E. Bacon 
and A. B. F. Duncan, j. Am. Chem. Soc. 56, 336 (1934). 
H. S. Taylor and J. C. Jungers, J. Chem. Phys. 2, 375 
(1934). L. Farkas and P. Harteck, Zeits. f. physik. Chemie 
B25, 257 (1934). 


sumption that each molecule decomposing pro- 
duces an increase of one molecule in the gas. The 
rate of the HCI—Dz» exchange was determined 
through the analysis, using the gas density 
balance, of the hydrogen-deuterium mixtures ob- 
tained. The quantum yield of the exchange 
reaction (number of H atoms exchanged per 
quantum absorbed) over any exposure interval is 
then given by 


2(F2—F;)Pp,+np+n,Pnu, 
(t2—t) (Inci/Inu,)(dPwu,/dt) 


where F, F: are the fractions of hydrogen in the 
hydrogen-deuterium mixture at the beginning 
and end of the exposure interval, ¢; to te (min- 
utes); Pp,+nHp+n, is the total pressure of the 
D:+HD+He: mixture; @yy, is the quantum 
yield of the ammonia decomposition at the pres- 
sure used; dPyu,/dt is the observed rate of pres- 
sure increase, in mm/min., when ammonia is 





(1) 


obs 

















30 
EXPOSURE (MIN.) 


Fic. 1. Hydrogen-deuterium analyses as a function of 
exposure time. Curve A, for a starting mixture of H.+DCl, 
curve B, D.+HCIl. 
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TABLE I. 








A. Initial rates in HCI+Dz2 mixtures 























IT “ 
Pp, +HD+H:| Pucr | HC |a-n 
a +H: mm | Jsta |min.| Fi | Fe | Foo | Pops} # 
31 79 | 1 3 | 0.153] 0.270}0.72 | 1.4] 1.9 
50 81 | 1 3 | 0.153] 0.268] 0.61 | 2.2] 3.1 
53 126 | 1 3 | 0.095/ 0.230| 0.68 | 2.8] 3.5 
54 330 | 1 5 | 0.095] 0.295] 0.875| 2.5] 3.2 
56 49 | 1 4 | 0.095] 0.232|0.42 | 2.2} 3.2 
62 234 | 1 5 | 0.095| 0.286| 0.775| 2.8| 3.6 
63 176 | 1 10 | 0.111] 0.432! 0.725) 2.3] 3.7 
75 85 | 0.33 5 |0.095| 0.198] 0.48 | 5.4] 7.0 
76 84 | 0.063 | 10 | 0.095] 0.200] 0.475] 14.7 | 19.0 
78 137. | 0.118 | 14 | 0.095] 0.275|0.60 | 9.9 | 13.5 
80 53 | 1 5 | 0.095] 0.234| 0.36 | 2.6] 4.1 
113 203 | 1 5 | 0.210] 0.344| 0.64 | 3.5] 5.4 
161 82 | 1 3 | 0.125] 0.203] 0.335} 4.9 | 7.0 
B. Initial rates in DCI+H2 mixtures 
Pp,4HD+H:| Ppci | (PC! |a=n 
mm mm Istq | min.| Fi F2 Fo | Pops} Po 
90 175 | 1 3 | 1.000] 0.853) 0.62 | 5.1] 6.5 
91 90 | 1 3 | 1.000] 0.907/ 0.76 | 3.3] 4.2 
175 74 11 3 | 1.000] 0.92110.85 | 5.3] 7.6 
































placed in the reaction cell; and Juci, Jnu, are the 
incident intensities of 2000-2140A when HCl 
and NH; were in the cell, with complete absorp- 
tion assumed in both cases. 

Since the quantum yield of ammonia decompo- 
sition varies with pressure,’ a standard ammonia 
pressure of 600 mm was adopted, at which @yu, 
was taken as 0.18. Runs with ammonia were made 
at a standard light intensity at which dPyu,/dt 
=0.31+0.01 mm/min. When using HCl—Dz 
mixtures, intensities other than this were ob- 
tained by varying the distance of the source and 
assuming the inverse square law. Introducing 
these modifications, (1) becomes 


1.16( F2— F;) Pp,+np+Hn, 
Pops = ’ (2) 
(t2—t) (Inci/Tsta) 


where 1.16=(2 X0.18)/0.31, Ista is the intensity 


3 E. O. Wiig, J. Am. Chem. Soc. 57, 1559 (1935). 
4 With a cell volume of 150 cc this corresponds to the ab- 
sorption of 1.410!” quanta/sec. 
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for which dPyu3/dt=0.31 mm/min. and Jc) is 
the intensity used in producing the exchange. 

As the exposure time is lengthened, or as the 
photostationary state is approached, 15 will 
decrease, and to obtain the true initial quantum 
yield it is necessary to extrapolate to zero expo- 
sure time and zero exchange. The most satis- 
factory way of doing this was found to be that 
based upon the empirical relationship 


},«(F,—F,), 
which leads to the equation 
1.16Pp,4np+n,/ sta( Po — Fo) - F.— sal 
Tuci(1 — Fo) (t2—t) Fo— Fe 





Po= (3) 


©,, By are the quantum yields at the times ¢=/ 
and t=0. F., and Fy are the fractions of hydrogen 
in the hydrogen-deuterium mixture at the photo- 
stationary state and at the time ¢=0, respec- 
tively. The factor (1— Fo) corrects for the hydro- 
gen content of the deuterium at =0. This equa- 
tion gives satisfactory constants over different 
time intervals of a given run, and the ®po values 
agree approximately with those obtained by 
cumbersome semi-theoretical methods. The extra- 
polations are from concentration ranges in which 
the reaction De+HCI—HD+DCI predominates 
over the second-step reaction, HD+HCI—H: 
+DCl, and the precision of the results does not 
warrant any attempt to apply a correction for 
the effect of the latter. Values of ®o are given in 
Table I. Since 49 was estimated from the first 
time interval of each series, F, (of Table I) = F). 

Experimental determinations of the amount of 
exchange at the photostationary state, F.., are 
compared, in Table II, with the value of F for 
the calculated thermal equilibrium at the tem- 
perature of the experiments, 25°C. The photo- 
chemical values were obtained by exposing the 
starting mixtures shown until theanalyses became 


TABLE II. Comparison of photochemical with thermal equilibrium. 














FRACTION OF H In D2+HD-+H2 at 
INITIAL ANALYSES EQUILIBRIUM (25°C) Ratio 
F . 
Pp.4+HD+H: | PHCI+DCI eee. 
MM MM D:+HD +He HC1+DCl PHOTOCHEMICAL THERMAL F thermal 
53 126 90.5% D 100% HCl 0.68 0.636 1.07 
113 203 79.0% D 100% HCl 0.64 0.610 1.05 
90 175 100% H 93% DCI 0.62 0.587 1.06 
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constant; values of the equilibrium constants 
used in estimating F for the thermal equilibrium 
were 
LHD ?/CH2 J[D2]=3.27, 
CHD JCHC1)/CH2 J[DCI] = 1.27, 
[HD J[DC1)/[D2 J[HC1]=2.58. 


DISCUSSION 


The initial quantum yields of greater than 
unity, together with the effect of light intensity 
(Fig. 2) indicate without question that the ex- 
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Fic. 2. Quantum yield as a function of light intensity. 
Solid line is for inverse proportionality between ®o and the 
square root of incident light intensity. Points are 
experimental. 


change is propagated by atomic chains which are 
stopped by recombinations. Using 102.7 and 
101.7 kcal. per mole as the dissociation energies 
of H, and HCI, and taking account of the zero 
point energies (HCI=4.2, DCI=3.0, H2=6.2, 
HD=5.4, and D.=4.4 kcal. per mole) one ob- 
tains the following energetically favorable reac- 
tions for starting mixtures of HC]+Dz: 


HCl+hyH +Cl, (1) 
H+HCI—-H.2+Cl+1.0 keal., (2) 
Cl+D.—-DCIl+ D — 1.6 kcal., (3) 
D+HCI—-HD+CI+1.8 kcal., (4) 
H+D.—-HD+D—1.0 kcal., (5) 
D+HCI—DCI+H+1.2 kcal., (6) 
H+H+M-—H:2+M +102.7 kcal., (7) 
H+D+M—HD+M+103.5 kcal., (8) 
D+D+M—D.+M-+104.5 kcal., (9) 


H+Cl+M—HCI+M +101.7 kcal., 
D+Cl+M-—DCI+M-+102.9 keal., 


(10) 
(11) 
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C1+Cl+M-—-Cl.+M+56.9 kcal., (12) 
H+Cl.—HCI+ Cl+ 44.8 kcal., (13) 
D+Cl.—DCI+Cl1+46.0 kcal. (14) 


For long chain lengths, these reactions give rise 
to an expression of the following form for the 
quantum yield: 


B[D:2][HCI o[ HCl 
ent [D2 ][HC!]((D2]+e[H ) (4) 
PM YL BmLD2}"[HCl}")! 


m+n=4 





where A varies between 0 and 2, depending upon 
the relative velocities of the chain carrying reac- 
tions above. B, a, and the 8,,,’s are also functions 
of the rate constants of the above reactions. The 
complicated dependence of #) upon the concen- 
trations of deuterium and hydrogen chloride as 
given by Eq. (4) may be approximated by simpler 
forms in case certain discrepancies in order of 
magnitude of the chain carrying rate constants 
prevail, and provided one may assume that all the 
recombination rate constants are of the same 
order of magnitude. For example, if (ke, ks, ke) 
> (ks, ks), ko>>(Rs, Rs, Re)>>ks, (Rs, ke) >> (Rs, Rs) 
>ke or ks> (Re, Ra, ke) >ks, Py a [De], whereas if 
the preceding orders of magnitude are reversed, 
Py a [HCI]. 


From the observation that the quantum yield 
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Fic. 3. Quantum yield as a function of the concentration 
of D2 and HCl. Solid curve for #*=0.073[D2]; O, experi- 
mental values of the quantum yield, corrected to the 
hypothetical conditions of 100 mm total pressure and 100 
mm HCI pressure, using ®o « [HCI ]*4([D2 ]+[HC1])~°**. 
Broken curve for #* = 1.15[HCI ]°*; ©, experimental values 
corrected to 100 mm total pressure and 100 mm D, pres- 
sure using )«[D»]({[D2]+[HCI1])~**. Corrections are 
based upon the empirical formula @9* =0.115[D»2 ][THCI1]°4, 
where #,* is the quantum yield o, reduced to the basis of 
total pressure = 100 mm with the assumption that M=[D, ] 
+[HC1], i.e., that deuterium and hydrogen chloride are 
equally effective third bodies for the recombination reac- 
tions, and that recombinations at the walls are relatively 
unimportant. 



























Fic. 4. General scheme for the reaction sequence, show- 
ing the various possible chain mechanisms for the exchange. 
Figures correspond to numbers of the reaction equations 
(see text). 


is pronouncedly more dependent upon the deu- 
terium concentration than upon the hydrogen 
chloride concentration, as shown in Fig. 3, it is 
concluded that either (3) is slow compared to 
(4), or (5) is slow compared to (6), or both; de- 
pending upon which of the numerous chains is 
contributing most to the reaction under the pre- 
vailing experimental conditions. 

Several investigators have found values near 
6.0 kcal. for the activation energy of the reaction 
Cl+H:2—HCI+H.$ Since this reaction is 1 kcal. 
endothermic and is the reverse of (2), the activa- 
tion energy of (2) may be taken as very close to 
5.0 kcal. The activation energy of (5) has been 
reported as 6.5 kcal.,* so it is at once indicated 
that most of the H atoms react by (2), and that 
chains involving (5) may, with reservations, be 
excluded from further consideration. The activa- 
tion energies of the Cl+DH and Cl+Dz ex- 
changes have been estimated to be 0.6 and 1.2-1.4 
kcal. higher, respectively, than the activation 
energy of Cl+H2.°:7 These estimates, combined 
with the zero point energies, give 7.2 kcal. for 
the activation energy of (3) and 4.6 kcal. for (4).° 
Having as yet no estimate of the activation 
energy of reaction (6), which competes with re- 
action (4) for the deuterium atoms, one must ad- 
mit the possibility of the two chain processes, 
(2)L(3)(4) ]+(3)(4) and (2)(3)(6)+(3)(6)(2), 
which remain after the elimination of reaction 
(5) from the general scheme of reaction sequence 
shown in Fig. 4. Experimentally, reaction (3) is 

5 J. C. Potts and G. K. Rollefson, J. Am. Chem. Soc. 57, 
1027 (1935) give 5.8 kcal., and calculate 5.88 kcal. from the 
data of Hertel, Zeits. f. physik. Chemie B15, 325 (1931). 
W. H. Rodebush and W. C. Klingelhoefer, J. Am. Chem. 


Soc. 55, 130 (1933), give 6.1 kcal. 
( : 35) Farkas and L. Farkas, Proc. Roy. Soc. A152, 124 
1 ‘i 
7H. W. Melville, Science Progress, 123, 499 (1937). 
§ Arbitrarily assuming that the DHCI complex has 0.2 
kcal. lower zero point energy than the HDCI complex. 
(Reference 6.) 
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the rate determining step, in keeping with the 
known activation energies for the first chain 
process. If (6) were slower than (3), the second 
chain process would be of negligible importance 
in the exchange, and the dependence upon [D, | 
would still maintain. 

For the exchange between DCI and Hg the 
chain carrying reactions are 


DCl+hv—-D+Cl, (ir) 
D+DCl—-D.+Cl+1.6 kcal., (2r) 
Cl+H.—-HCI+H —1.0 kcal., (3r) 
H+DCI—-HD+C1+0.6 kcal., (4r) 
D+H.—-HD+H-+0.8 kcal., (5r) 
H+DCI—-HCI+D —1.2 kcal. (6r) 


Taking the activation energy of Cl+Dz as 7.2 
kcal.,° one obtains an activation energy of 5.6 
kcal. for (2r). The observed activation energy of 
(3r) is 6.0 kcal., while the same value is found for 
(4r) from Cl+DH=6.6 kcal.6 The activation 
energy of (5r) has been estimated as 4.8 kcal.° In 
this case the deuterium atoms should react chiefly 
by (5r) rather than by (2r). The chain carrying 
steps obtained by eliminating (2r) from the gen- 
eral scheme of reaction sequence analogous to 
that shown in Fig. 4, and assuming (5r) rela- 
tively faster than (3r) become chiefly (5r)(6r), or 
(5r)[(4r)(3r) ]+(3r)(4r), or both, depending 
upon the relative speeds of (4r) and (6r). 
Although our results on initial mixtures of 
DCI+He are unfortunately meager, they defi- 
nitely indicate a dependence upon the DCI con- 
centration nearly as great as that upon the H: 
concentration, and a quantum yield of the same 
order of magnitude as for mixtures of HCI and D>. 
This can be interpreted as favoring the second 
chain process given above. For the first chain 
process to give nearly equal dependence upon 
[D2] and [HCI] would require (6r) to have an 
activation energy near to that of (5r), namely, 
4.8 kcal. This would indicate that the exchange 
between DCI and He should be much faster than 
that between HCl and Ds, for which the rate 
determining step has an activation energy of 7.2 
kcal. It appears more plausible that the second 
process, carried by chains of (3r)(4r), both having 
activation energies of 6.0 kcal., accounted for the 
observed exchange. Hence (6r) must have an 
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activation energy >6 kcal. and (6) >4.8 kcal., 
indicating that reaction between H (or D) and 
DCI (or HCl) tends to produce HD and Cl rather 
than HCI (or DCl) and D (or H). 

Of the large number of possible mechanisms 
for the exchange reaction between HCI and Dz, 
all of which are of nearly the same order of im- 
portance, the results reported here appear to 
favor chains of reactions (3) and (4) as having 
the greatest influence on the kinetics of the photo- 
chemical exchange under the experimental con- 
ditions employed. Similarly, for DCl+Hz2 mix- 
tures, chains of (3r) and (4r) are indicated. Before 
these conclusions can be taken as final there are 
factors which may require further study. These 
include the effect of wall reactions, the effect of 
variations in the rates of the chain terminating 
reactions, the possible difference in the efficiency 
of HCl and Dy as third bodies in the atomic re- 
combinations, and the possibility that the activa- 
tion energies may be in error by an amount suffi- 
cient to alter the conclusions drawn from their 
relative magnitudes. Also, the entire treatment 
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has assumed sufficient diffusion to insure that the 
local and over-all rates are the same. 

At the photostationary state the forward and 
reverse reactions given above, together with the 
reactions of the H, D, and Cl atoms with HD 
molecules, reach a state of balance in which the 
differences in rates should result in the attain- 
ment of thermal equilibrium proportions. Actu- 
ally the photostationary states are displaced in 
the direction of higher hydrogen content in the 
D.+HD-+He mixture than that corresponding to 
thermal equilibrium (Table II). Although it is 
uncertain just what significance should be at- 
tached to this difference, it probably arises from a 
stronger absorption by HCI than by DCI in the 
region 2000—2150A, which is to be expected from 
the lower zero point energy of DCI. A difference 
of about 10 percent in the absorption of the zinc 
lines by HCI and DCI should be sufficient to pro- 
duce the observed shift of the equilibrium. 

The authors are indebted to the Carnegie Insti- 
tution of Washington for a grant in support of 
this research. 
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It is pointed out that chemiluminescent reactions are 
necessarily nonadiabatic, since emission of radiation can 
only occur by transition from an upper to a lower potential 
energy surface. The restriction imposed by the second 
law of thermodynamics on the photon yield is, 


Z 
eK AF/E, 

where ¢ is the ratio of quanta evolved to molecules reacted, 

AF is the free energy decrease in the reaction, and E is 


INTRODUCTION 


HE purpose of the present paper is to treat 

chemiluminescence as a particular example 

of the general theory of chemical reactions. 

A general survey of the field is given along with 

. Sufficient examples to illustrate the point of view 
and the method of treatment. 


the energy of the photon emitted. Upper potential energy 
surfaces are reached either by activation by radiation, 
which may be that of a black body at the temperature 
of the reaction, or by the ‘‘crossing”’ of potential energy 
surfaces. From Audubert’s data on the radiation produced 
by the thermal decomposition of sodium azide approxi- 
mate surfaces which account for the observations are 
constructed. 


Chemiluminescent reactions are noteworthy 
because of the wide variation of the ratio of 
the number of quanta evolved to the number of 
molecules reacted. Thus reactions are found 
experimentally with photon yields ranging from 
about 0.3 to 10-'®. The upper limit of 0.3 for the 
photon yield is reported for the reaction of sodium 
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vapor with chlorine, a type of reaction studied 
by Polanyi and his co-workers.'! An intermediate 
value for the yield has been observed by Harris 
and Parker’ for the oxidation of 3-aminophthal- 
hydrazide, while reactions with very low photon 
yields have been studied by Audubert and his 
associates.? Those with the very high photon 
yields are plainly visible to the unaided eye, 
while the rather elaborate and sensitive photon 
counters employed by Audubert and his associ- 
ates must be employed to detect the radiation 
emitted by those with very low yields. There is 
no apparent reason why reactions with photon 
yields anywhere between zero and one should not 
be found, subject to the thermodynamic re- 
striction discussed below. 


EFFICIENCY OF RADIATIONLESS TRANSITIONS AND 
THE SECOND LAW OF THERMODYNAMICS 


In this section are discussed certain thermo- 
dynamic restrictions applying to the efficiency 
with which normal or excited atoms or molecules 
can be further excited by ordinary thermal 
collision processes. Consider a mole of mercury 
atoms in the metastable 6*P» state dispersed in a 
foreign gas. Then the energy, E,; per mole, 
corresponding to the transition from the excited 
state, 6*P», to the normal state, 6'So, is almost 
exactly equal to the free energy difference 
between the two states. This is true because if 
the energy of excitation were emitted, the work 
which it could be made to do would simply be 
equal to the energy difference between the two 
levels. Now suppose that a number, Nj, of the 
excited mercury atoms are raised by thermal 
collision from the state, 6*Po, to the still higher 
state, 6°P:, where their energy referred to the 
ground state is E2. In such a process an amount 
of energy, Ni(Z2—£)/N, is converted from the 
thermal energy of the molecules of the inert gas 
to the free energy of excitation of the mercury. 
Clearly this increase in free energy for part of 


1 Beutler and Polanyi, Zeits. f. Physik 47, 379 (1928); 
Zeits. f. physik. Chemie, Bl, 3 (1928). Bogdandy and 
Polanyi, Zeits. f. physik. Chemie B1, 3 (1928). Polanyi and 
Schay, ibid. B1, 30 (1928). 

2 Louis Harris and Almon S. Parker, J. Am. Chem. Soc. 
57, 1939 (1935). 

(ssey’ for example, R. Audubert, J. Chim. Phys. 33, 507 






the atoms must be balanced by a decrease, at 
least equal in amount, caused by radiationless 
transitions from the first excited to the ground 
state of the mercury atoms. The decrease in free 
energy, caused by Np» such transitions, is given 
by, N2(Ei/N), and we have the inequality, 


Ni(E.—E1) | NoEs 
< , 
N N 





(1) 


Since N,; plus Ne equals N, this leads to the 
relations, 


N, ky No (E2— EF) 
—<€— and —>———. 
E2 N Es, 


Thus we see that the second law of thermo- 
dynamics places an upper limit on the efficiency 
of excitation by radiationless transitions, and it 
follows that the number of quanta of energy, 
E:2/N, that can be emitted is restricted in the 
same way. In reactions of this kind which do 
not involve the formation of a new kind of 
molecule no concentration terms should be in- 
cluded in the free energy expression since normal 
and excited atoms or molecules cannot be re- 
garded as separate components. If they could be, 
then energy rich molecules could conceivably be 
separated from those of low energy by a mem- 
brane permeable to only one of the species. 
Such a membrane, however, would be an exam- 
ple of Maxwell’s demon. 

The above example is a special case of the 
more general law, 


(Ni/N) < (AF/E), (3) 


where AF is the decrease in free energy in any 
chemiluminescent reaction, E is the energy perf 
Einstein of the quanta which are emitted, and 
N,/N is the photon yield. It may be seen from 
(3) that in those cases where AF>E these 
thermodynamic considerations place no restric- 
tion on the quantum yield. 

The AF of Eq. (3) is the decrease in free 
energy that actually takes place in the reaction 
vessel ; i e., it includes a term for the free energy 
of mixing. Thus, even if the equilibrium constant 
is unity, and consequently the free energy 
difference of the pure products and reactants is 
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zero, there is still a free energy decrease in the 
reaction due to the mixing of reactants and 
products. As an example let us consider the 
problem of how much radiation can be evolved 
if a mole of A reacts to form 3 mole of its optical 
isomer, B, in ideal solution with A. Here the 
equilibrium constant is certainly unity and the 
free energy change is due entirely to the entropy 
of mixing. This is equal to R In 2, which corre- 
sponds to a free energy decrease of 0.415 kcal. at 
300°K. The maximum photon yield from such 
a reaction is found by application of Eq. (3) to 
be about 1 percent for a wave-length of 6900A 
and correspondingly different for other wave- 
lengths. Since actual reactions in general have 
much larger free energy decreases, as well as 
photon yields much below this value, it is 
obvious that the restriction imposed by Eq. (3) 
will rarely, if ever, be the factor limiting the 
yield. When pure A first forms B the free energy 
decrease per molecule is extremely large so that 
initially there is no thermodynamic restriction on 
the quantum yield. As the concentration of B 
approaches equality with that of A the thermo- 
dynamically allowed gy drops to zero. Our 
definition has been couched in terms of the over- 
all AF and ¢. 

The question has sometimes been asked 
whether or not it is possible for an endothermal 
reaction to be chemiluminescent. The answer 
given here is that if the reaction is spontaneous, 
ie., if it involves a free energy decrease, then it 
can be chemiluminescent with a maximum 
photon yield given by Eq. (3). How nearly this 
maximum efficiency is approached is purely a 
question of mechanism and can only be solved 
by reaction rate theory. 


ACTIVATION BY BLACK-Bopy RADIATION 


Since chemiluminescent reactions must neces- 
sarily involve the transition of molecules from 
upper to lower potential energy surfaces, the 
theory of these reactions must provide some 
mechanism for reaching an upper level. One such 
possible mechanism is activation by means of 
absorption of black-body radiation originally 
Proposed by Perrin and W. C. McC. Lewis as 
the mechanism for the activation of molecules 


CHEMILUMINESCENCE 














Reaction Coordinate —» 


Fic. 1. Schematic energy diagram on which a chemilumi- 
nescent reaction could take place by means of activation by 
black-body radiation. 


undergoing ordinary chemical reactions. There 
is abundant evidence showing that the reaction 
rates via this mechanism are ordinarily slower by 
a factor of about 10° than the rates resulting 
from activation by thermal collision when the 
activation energies are the same. However, 
activation by black-body radiation may be of 
importance in chemiluminescent reactions. 

Figure 1 illustrates a potential energy surface 
on which it would be possible for a chemi- 
luminescent reaction to take place by means of 
the black-body mechanism of activation. 

Here A; (lower) and A, (upper) refer to the 
species that are reacting in the normal and 
excited states respectively; B; and B, have the 
same significance for the products of reaction. 
A,* and A,* indicate the activated molecules for 
the lower and upper surfaces respectively. When 
enclosed in parentheses, these symbols refer to 
concentrations. The k’s indicate the specific 
reaction rate constant for a process proceeding in 
the direction indicated by the arrow. When the 
reaction rate constants refer to directions oppo- 
site to that indicated by the arrow, they will be 
primed. 

At equilibrium between A, and A), the rates 
in the two directions will be the same and a 
discussion applying to k1’ will be applicable to k:. 
An expression for k;’ can be written in terms of 
the Einstein coefficient for spontaneous transi- 
tion, A,:, from the upper state corresponding to 






















































A,, to the lower state, A;. We have* 





hy’ Ez _ 
_= Wy 
7 a7 1—exp (—hv.1/kT) 


Xexp (—«/T)} {>ow, exp (—e./kT)}—. (4) 


Here w,, is the usual statistical weight, v,: is the 
frequency of the emitted radiation, and e, is 
the energy in the excited state referred for con- 
venience to the lowest state of the normal 
molecule, A;. If the radiation density is equal 
to that of a black body at the temperature of 
the reaction, the states u and / will be in equi- 
librium and we have the relation, 


Equilibrium constant 
Vw: exp (—e€2/ kT) 
kj # 


=——=2 mone, 8) 
ky Yow, exp (—€./RkT) 





Equations (4) and (5) may be combined to 
give an expression for 3. 


Aui 
ky _ ; i Wu 
ul 1—exp (Avyi/RT) 


X {Lwrexp (—e/kT)}. (6) 
l 


exp (—e,/RkT) 





The expressions for k,’ and k, are identical 
with those for k;’ and k; except that in this case 
the summations are made over the energy levels 
of B, and B,, respectively. 

Since expressions (4) and (6) represent the 
total transition probability between states we ex- 
pect that a satisfactory approximate expression 
can be secured by taking a transition probability 
averaged over all states, and an averaged energy 
of activation secured in the same way. This 
gives for the rate expression, 


ki=(1/r) exp (—E/kT). (7) 


Here 1/7 and E are interpreted as being the 
reciprocal mean lifetime and energy difference 


4See for example, Condon and Shortley, The Theory of 
Atomic Spectra, (Cambridge Press, 1935), Chap. IV, or 
Tolman, Statistical Mechanics with Applications to Physics 
= — (Chemical Catalog Co., 1927), Chaps. 16 
an , 
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between the two states, respectively, and the 
difference in entropy between the upper and 
lower states has been neglected. When sufficient 
information is available the exact expression (6) 
will be used. For many atomic transitions, 1/7 
is of the order of 10° and is probably not greatly 
different for transitions in molecules. For for- 
bidden transitions it may drop to values as small 
as 10-?.5 

We now go on to develop a general theory of 
chemiluminescent reactions involving as one step 
the transition from a lower to an upper level by 
means of black-body radiation. This can be done 
conveniently by referring to the schematic 
energy diagram in Fig. 1. When there is not too 
much difference between the free energy of the 
activated complex, A,* on the upper surface 
and A,* on the lower surface, an appreciable 
fraction of the reaction will proceed by way of 
the upper surface and will therefore emit radia- 
tion in the transition, B,—B,;. At equilibrium 
this will just be balanced by the absorption 
B.-B,,, but in a system containing more than 
the equilibrium amount of reactant, A, more 
energy of this wave-length will in general be 
emitted than absorbed. 

As an example we set up the rate and photon 
yield expressions for the particular case where all 
the reactions are first order with respect to the 
substances indicated in Fig. 1. Any dependence 
on other substances is implicitly included in 
the k’s. We can write in the usual way. 


d(A u)/dt _ k(A 2) +h, (B.) — (A u) (Ry +s) = 0, (8) 
d(B,)/dt=ks(Au)—(Bu)(Ri +kst+ks)=0. (9) 


Here the processes with specific rate constants, 
k;’ and k;’ have been neglected. Combining 
Eqs. (8) and (9) so as to eliminate (A ,,) we have, 


ksk1(A 1) 


(B.) = (10) 
(ki +h,) (Ra +hs+hs) — aka’ 


Then if dN/dt is the rate of emission of light we 
have, 





dN k3k4k1(A1) 
pease (11) 


dt (k1' +h,) (Ra! +hs-+hs) — Rake’ 


5 Condon and Shortley, reference 4, p. 282. 
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Reaction Coordinate —» 

Fic. 2. Schematic energy diagram illustrating the 
“crossing’’ mechanism of activation. The separation of the 
two surfaces at the point corresponding to A,* is exag- 
gerated. 
and for the photon yield, ¢, we obtain, 

(dN/dt) 
—d(A,)/dt 
ks 
Rol (Ri’ +s) (Ra +Rs+ks) — Rak’ ] 
kik, 








if we note that d(A,)/dt is given by 
—d(A;)/dt=dN/dt+k;(Bu)+k2(A1). (13) 


We will consider three limiting cases for the 
rate dN /di, in each case neglecting the quenching 
process involving ;. 

Case I.—When the rate determining step is the 
transition A;—A, then k,’<k, and 


g=ki/(kitke). (14) 


Case II.—When the rate determining step is 
the transformation A,—B,, then ki<&ky', ka’ Kk; 
and, 


ea(ks/ky’) 
PT Ralka/ka!) + Re 


Case IJI.—Finally, when the rate determining 
step is the transition B,—B, then kk,’ and 


ks(Rika/Ri'ky’) 
les . 
k3(Riks/ki'k,’) +ke 





(15) 


(16) 





The quantities in parentheses in Eqs. (15) and 


(16) are equilibrium constants between upper 
and lower levels. 

From Fig. 1 it may be seen that a photon 
yield approaching the thermodynamically allow- 
able maximum given by Eq. (3) is not to be 
expected unless the upper surface at the point 
marked A,* approaches close to the lower sur- 
face. For higher temperatures where the rate 
corresponding to k; is the slow step, a high yield 
will be obtained if the energy difference between 
A, and A; is very small. 


ACTIVATION WITHOUT BLACK-Bopy RADIATION 


In this case chemiluminescence will arise be- 
cause of the nonadiabatic ‘‘crossing’’ of energy 
levels. This is illustrated in the schematic 
diagram (Fig. 2). 

The essential difference between Figs. 1 and 2 
is that in the latter case the two potential sur- 
faces come very close together and to a first 
approximation they may be said to cross making 
possible a nonadiabatic jump of the system from 
the lower to the upper surface without absorbing 
radiation. The approximate expression for the 
probability of this process happening has been 
given by Zener® and by Landau’ and the ex- 
pression has been applied in a treatment of the 
nonadiabatic decomposition of N.O.* This ex- 
pression is 


P=exp | (17) 


—4r’ e210 


hV|si—Se| 


Here |si—s2| is the absolute difference in slope 
between the two dotted lines in Fig. 2 while 2¢). 
is the difference in energy between the upper and 
lower surfaces at the activated point correspond- 
ing to the species, A,*. V is the velocity along 
the reaction coordinate at this point. 

The transmission coefficient,? x, may be re- 
lated to P in the following way. If, of the number 
of molecules that reach the activated point 
corresponding to A,*, a fraction P jump to the 
upper surface, then 1—P will be the fraction 
that remains on the lower surface. Those on the 

6 C, Zener, Proc. Roy. Soc. A137, 696, (1932); A140, 660 
a Physik. Zeits. Sowjetunion 2, 46 (1932). 

8A. E. Stearn and H. Eyring, J. Chem. Phys. 3, 778 


(1935). 
®H. Eyring, J. Chem. Phys. 3, 107 (1935). 





upper surface may return to the lower surface in 
two ways, one corresponding to reaction and 
the other corresponding to a return to the initial 
configuration, A). x, the probability that the 
system will react after having once reached the 
activated configuration, is given by, 


=1—P+P{(1—P)+(1—P)! 
+(1—P)5++++]=2(1—P)/2—P, (18) 


which reduces to 2—2P when P is almost equal 
to unity. Further, when this is the case, the ex- 
ponential in Eq. (17) may be expanded and the 
higher powers of (47°e12")/(hV|si1—s2|) dropped 
so that we have as an approximate expression 
for xk, 


k= (87e12")/(AV|s1—S2|), (19) 


which is valid for values of P that are not too 
far from unity. 

One other possible mechanism for reaching the 
upper surface is that corresponding to the transi- 
tion with the specific rate constant, k;. This will 
in general, however, not be of importance for the 
case of activation by black-body radiation unless 
the upper surface comes very close to the lower 
one at a point directly above the one correspond- 
ing to A; on the lower surface. However, if 
activation is brought about by illuminating the 
system with radiation of the proper wave-length 
this mechanism for reaching the upper surface 
may become the important one. This would be 
a case of photoreaction which might be accom- 
panied by fluorescence. The general considera- 
tions discussed above also apply in this case. 

Considering now the rate determining step for 
the luminescent reaction we see that it can be 
either the process with specific rate constant, 
ky, or the one with specific rate constant, k;. 
When the mean life of B, is of the order of 10-8 
sec., then the rate corresponding to k, will be 
the slow step if the energy of activation for k,’ is 
less than about 237 cal. or 6900 cal. at room 
temperature. For energies of activation greater 
than 23T cal., ks will become the slow process. 
These approximate computations neglect any 
difference in entropy between B, and A,*. 

There is, however, an important group of 
reactions which involve the possibility of emitting 
radiation, but where the mean life of the excited 
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molecules may be of the order of minutes. This is 
the case for those instances where the ‘“‘crossing”’ 
of potential energy surfaces is due to a change in 
multiplicity. The general expression for the 
Einstein transition probability A(a, b) between 
level a and b can be written as the sum of three 
terms, i.e., one for the electric dipole, one for the 
magnetic dipole and one for the quadrupole 
radiation, respectively. The transition proba- 
bility can be written in the form :'° 


A(a, b) = (641403/3h)[ | (a|P! 0d) |? 
+ | (a@|M| 6) |?-+(3/10)x?o?| (a | N| 5) |?]. 


Here o is the wave number of the emitted 
radiation; 


=z edn; = (e/2uci) 2 (Lit28;) ’ 
= —eLriti; 


where e is the charge on the ith electron whose 
vector distance is r;; wis the mass of the electron; 
L; the orbital angular momentum, and §; the 
spin angular momentum. For forbidden lines 
|(a|P|b)| equals zero due to a change in 
multiplicity and A(a, b) instead of being about 
10° is much smaller. For example the interpre- 
tation of the nebular lines of OJJI as given by 
Condon" leads to an Einstein transition proba- 
bility of A=0.018 sec.—' for the D2—3Pz2 transi- 
tion and an A= 0.006 for the transition D.—3P,. 
The transition probabilities for quadrupole radia- 
tion are found to be about 1000 times smaller. 
For transitions involving a change in multiplicity 
the magnetic dipole transition probabilities are 
known to be proportional to (e12/(£,—£:))’ 
where ei is the spin orbital interaction between 
the upper and lower state and E,—£z2 is the 
energy emitted in the transition. Now this e12 is 
formally the same as the e12 occurring as a factor 
in the transmission coefficient for ke discussed 
above. The difference in the two ¢12’s is that they 
are the spin orbital interactions calculated for 
different configurations but an examination of 
multiplet widths (which enables one to calculate 
£12 for atoms)" indicates that the value for ex 

10 Condon and Shortley, reference 4, p. 90 etc. 

1 Condon, Astrophys. J. 79, 217 (1934). Condon and 


Shortley, reference 4, p. 282. 
12 See Condon and Shortley, reference 4, p. 268. 
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often changes only slightly when an atom 
becomes a part of a diatomic hydride molecule. 
There are notable exceptions to this rule among 
nonhydride molecules.'* Thus for forbidden mag- 
netic dipole transitions of B,—B; we expect 
(xo(kT/h)) often to be proportional to 1/73, the 
former being about 10!° times the latter. Thus if 
ks is to be of the same order as ko we must then 
have (Ei—E2)/RT ~ 10"; or E;—E2~46T cal., 
i.e., if the reactions are to proceed at about the 
same rate at room temperature then the energy 
of B, should be about 14 kcal. lower than A,*. 
From this one sees that there are types of 
“crossing’’ potential surfaces which can lead to 
high photon yields even when the emission 
process involves a change of multiplicity and is 
consequently forbidden in the electric dipole 
approximation. For. ‘‘crossing’’ potential surfaces 
which, however, do not involve forbidden transi- 
tions for B,—B,(xe(kT/h)) will probably be more 
nearly 10° than 10'° times 1/73. This case has 
already been discussed. 

The quenching of radiation can in general be 
expected to proceed by reaction over the potential 
hills marked k; in Figs. 1 and 2. If this is the case 
the amount of quenching for a given activation 
energy can be expected to be proportional to the 
mean life of B,. Quenching can thus be expected 
to be a much more important factor for cases 
involving forbidden transitions than for other 
instances where the mean lifetimes are very 
short. 

In Fig. 2 a particular manner of “‘crossing’’ of 
the two surfaces is shown. Clearly the photon 
yield for radiation will be higher the higher the 
activation energy hill for the nonluminescent 
reaction is in comparison to the luminescent one. 
If, instead of ‘‘crossing’’ before the system 
reaches the activation hill, A,*, the ‘‘crossing”’ 
occurs afterwards, the chemiluminescent effi- 
ciency should approach the thermodynamic 
maximum. In this case there would be no 
difference in activation energy for the lumi- 
hescent and nonluminescent reactions, and the 
entire difference in rate would arise from the 
difference in the transmission coefficients for the 
two alternative processes. 

It is of some interest to note that the thermo- 
dynamic considerations discussed in Section II 


* See R. S. Mulliken, Rev. Mod. Phys. 4, 37 (1932). 
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place a lower limit on the transmission coefficient 
corresponding to the reaction rate constant, ke 
of Fig. 2. Thus, x cannot become so small that 
the condition imposed by Eq. (3) will be violated. 


CHEMILUMINESCENT REACTIONS OF METALLIC 
AZIDES 


As an example of the interpretation of data 
that may be obtained, we consider a series of 
chemiluminescent reactions studied by Audu- 
bert,'* the decomposition of metallic azides. It 
appears likely that the potential energy surface 
for these reactions is similar to that given in 
Fig. 1, although a surface of the type illustrated 
by Fig. 2 is not entirely precluded. However, 
evidence will be presented tending to show that 
if the luminescent reaction proceeds by way of 
“crossing” of potential energy surfaces, the 
crossing cannot be due to a change of multiplicity. 

Audubert finds that the azides of silver, lead, 
sodium and potassium all emit radiation lying in 
the neighborhood of 2000A, about 140 kcal. per 
mole. This we interpret as corresponding to the 
transition, B, to B, of Figs. 1 or 2, and the fact 
that all the different azides emit light of about the 
same frequency seems to indicate that the same 
molecule is radiating in all cases. 

The decomposition of NaN; furnishes an 
example of what might be expected to be a fairly 
common phenomenon, a chemiluminescent reac- 
tion in which the activation energy for photon 
emission show two quite different values. The 
rate determining step at low temperatures is 
due to the reaction, A,—B,, proceeding with an 
activation energy of about 48 kcal. The activa- 
tion energy of 23 kcal. given in column four of 
Table I may be either that for the reaction of B,, 
going to B; shown in Figs. 1 or 2 or that for the 
transition, A; to A,, of Fig. 1. However, the 
following consideration will throw some light 
on the question. If the mechanism for the 
reaction is that illustrated by Fig. 2, then the 
“‘crossing’’ of the two surfaces cannot be due to 
change of multiplicity. For transitions involving 
change of multiplicity the value of 7 in Eq. (7) is 
greater by a factor of about 10! than for other 
transitions which are not forbidden. But the 


temperature at which the two rates corresponding 


14 Audubert, J. Chem. Physique 34, 29 (1937). 
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to the two activation energies become equal is the 
point at which the two straight lines obtained by 
plotting In Rk against 1/T intersect. By referring 
to Fig. 3 of Audubert’s paper this temperature is 
found to be 640°K for the case of NaN;. We can 
then set up the equation, 


kT —48+2 kcal. 1 —23+2 kcal. 
— exp ( ) =— exp ( ) 
RT RT 


e 








and solve for 7. For this case we get values for 7 of 
approximately 3X 10-!° and 10-* sec. for upper 
and lower limits, respectively. These values are 
even smaller than the value of about 10-® sec. 
that would be expected for most allowed transi- 
tions, and are very much smaller than would be 
expected for transitions involving change of 
multiplicity. We accordingly draw the conclusion 
that if the potential energy surface for the 
decomposition of the azides is similar to that 
illustrated in Fig. 2, the “‘crossing’’ does not 
involve a change in multiplicity. 

It is not possible to decide from the evidence at 
hand whether the process with the activation 
energy of 23 kcal. corresponds to the reaction of 
B,, going to B; of Figs. 1 and 2 or to the reaction 
of A; going to A, of Fig. 1. This point could, 
however, be decided experimentally by increasing 
the density of radiation at the wave-length 
corresponding to 23 kcal., about 12,000A. This 
would speed up the transition A; to Az, but 
would alter the velocity of B, going to B; by only 
a negligible amount. 

In confirmation of our suggestion that the 














TABLE I. 
ENERGIES OF ACTIVATION FOR THE 
Reactions, (See Fig. 1 or 2) 
ki 
A l Azle 
or 
Sus- HEAT OF ko ks ka 
STANCE REACTION!5¢ | A,—+B)1% By--Bi Al B,'e 
AgNs |— 67.3 kcal. 47+1 kcal. 
Pb(Ns3)2|— 103.0 38 kceal.| 25+1 kcal. 
NaN; |—_ 5.0(?) 34.4 2342 48+2 
KN; |—. 0.5(?) | 36.1 48+2 




















1a Data quoted by Audubert. For AgN3; and Pb(N3)2 the reference to 


the original work is, Wéhler and Martin, Zeits. ges. Schiess Spreng- 
stoffen, 12, 1 (1918). For the alkali azides Audubert gives no reference 
and we have been unable to find one. 

18 Garner and Mark, J. Chem. Soc. 657 (1936). 

use Audubert apparently omitted the factor of 2.3 in calculating his 
activation energies from the slope of the curves in which logi k is 
plotted against 1/T. 
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decomposition reactions of the azides all proceed 
by the same mechanism, it will be noted that the 
activation energies for the reaction A; going to 
B,, are the same for the three cases recorded in 
Table I. Likewise the activation energies for the 
other step in the luminescent reaction are 
essentially the same for the two cases recorded. It 
appears probable that if the decomposition 
reaction of Pb(N3)2 could be followed to a low 
enough temperature an activation energy for the 
reaction of A; going to B, would be found. 
Conversely, if the decomposition reactions of 
AgN; and KN; could be studied at high enough 
temperatures another activation energy should 
be found with the value of about 25 kcal. 

We now attempt to construct as complete a 
potential energy surface for the decomposition of 
NaN; as is possible by using the data at hand. A 
possible surface that accounts for most of the 
facts is shown in Fig. 3. 

The amount of free energy, 26+2 kcal./mole, 
by which the maximum of the upper curve lies 
above the maximum point of the lower one is 
obtained from Audubert’s value of the photon 
yield. He states that the order of magnitude of 
this quantity is 10-® to 10—'° for the azides which 
leads to the equation: 


exp — (AF.t—AF2t)/RT=10~ to 10-", 
AF w=AF.t+26+2 kcal./mole. 


To the approximation to which we can neglect 
T(ASit—ASet) this is also the difference in 
energy of activation for the two reactions. 

The difference in energy between the A,* and 
A,* can be obtained by a second method. 
This energy difference is found to be (48+2) —34 
= 14+2 kcal./mole from the temperature coeffi- 
cient of the photon yield. The disagreement 
indicates either that the value of about 10~° for 
the photon yield is too low, or that it is not a 
good approximation to neglect the difference in 
entropy between the activated complexes on the 
lower and upper surfaces. 

The emission of radiation of a frequency corre- 
sponding to an energy of 140 kcal. per Einstein 
presents an interesting problem. The energy 
difference between A; and B, according to the 
data cannot be greater than about 23 kcal., and 
may be less than this depending on whether the 
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transition A; to A, or B, to B; is the slow process. 
Audubert quotes the figure 5 kcal./mole for the 
heat emitted during the decomposition of NaN; 
into solid sodium and gaseous nitrogen. If this 
figure is correct and if there is no energy differ- 
ence between the stable final product and the 
product obtained immediately after emission of 
the photon, then 23 molecules of NaN; would 
each have to contribute their energy of reaction 
in order to build up an excited molecule capable 
of emitting a photon of the required energy. 
On the other hand, this difficulty is not en- 
countered in the case of AgN; and Pb(N3)o. In 
these cases it is only necessary to assume that 
two azide groups are involved in the steps leading 
to the emission of each photon in order to have 
sufficient energy available. 

Since all the reactions proceed with about the 
same quantum yield it does not appear to be 
probable that the luminescent reaction for the 
azides of silver and lead involve only two azide 
radicals per quantum emitted while the same 
reaction for the azides of the alkali metals in- 
volves twenty molecules or more. What appears 
to us to be more probable is that the quoted 
values for the heats of decomposition of the alkali 
azides are incorrect. 

The work of Sklar'® on the theory of color of 
organic compounds makes it possible to offer a 
tentative suggestion concerning a possible struc- 
ture of the excited molecule that emits the 
quanta of wave-length 2000A. Sklar expresses 
the energy levels corresponding to the absorption 
of light of different frequencies in terms of the 
resonance energy of the benzene molecule. He 
finds the first strong absorption band to be at a 
wave-length of about 2000A, just the same 
wave-length as the light that is emitted during 
the decomposition of the azides. This fact, along 
with the observation that in the case of the 
azides of lead and silver the energy of two azide 
radicals is necessary to furnish the energy 
emitted in one photon, leads us to suggest that a 
cyclic nitrogen compound with six nitrogen atoms 
arranged in a ring may be the molecule emitting 
the radiation. Although the total energy differ- 
ence between cyclic Ng and 3N2 must be very 
much less than that between CsH¢ and 3C:Hz, 
it is not unreasonable to assume that the reso- 


* Alfred L. Sklar, J. Chem. Phys. 5, 669 (1937). 
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Fic. 3. Schematic energy diagram for the chemilumi- 
nescent decomposition of sodium azide. The energy units 
are kcal./mole. A surface similar to Fig. 2 is not completely 
excluded. 


nance energy of the cyclic N¢ would be about the 
same as that of benzene. If this were the case 
the absorption bands according to Sklar’s theory 
would be expected to be about the same. 

It should be pointed out that the energy of 
the reacting system immediately after emitting a 
photon is not necessarily greater than the energy 
of the final products. This, of course, is not true 
of the free energy which must be lower for the 
stable products than for any intermediates. In 
this connection it may be pointed out that if 
the value of 5 kcal. for the heat of decomposition 
of NaN; is accepted and if the cyclic Ng molecule 
is the immediate product after the radiation is 
emitted, then the Ng molecule would be expected 
to be thermodynamically stable up to about 
1300°K. This estimate is based on the assump- 
tion that the difference between the entropy of 
N, and 3N¢z is about the same as the difference 
in entropy between CsHg and C2He. Even if 
the reaction mechanism is quite different from 
that proposed here, an energy difference of the 
order of magnitude of that shown by the dotted 
curve of Fig. 3 seems to be entirely unreasonable. 


CONCLUSION 


Although quenching processes play an im- 
portant role in most reactions involving radia- 
tion, they have been only briefly mentioned 
above. There are many data available on the 
quenching of excited atoms in the gas phase by 
collision with foreign atoms or molecules. These 
















































data have generally been interpreted in terms of 
an effective cross-sectional area, but a more 
rational explanation is to be found in the 
“crossing” of potential energy surfaces. For 
many simple cases the surfaces could be con- 
structed for direct comparison with experiment, 
and work along these lines is planned. Aside from 
these cases, there are quantitative data on 
quenching for only a few cases. Anderson!® has 
obtained data on the effect of dissolved salts on 
the chemiluminescence of cypridine luciferin, and 
has found increases in the photon yield in some 
cases. This result is understandable either as a 
change in the activity coefficients of the re- 
actants or activated complexes for some of the 
reactions or as a change in the important 
mechanisms. One possible example of the latter 
type would be the removal of a quenching agent 
originally present. 

The phenomenon of mitogenetic ray emission 
has been very extensively studied with rather 
contradictory results, and a review has been 
written by Hollaender and Claus." Although we 
do not attempt to express an opinion regarding 
the physiological effects of such radiation, the 
occurrence of many weak chemiluminescent re- 
actions in living organisms is to be expected in 
the light of the mechanism for their occurrence 
discussed above. 

Audubert and Victorin'® have investigated the 
ultraviolet luminescence occurring at aluminium 
anodes showing an electric valve effect. While it 
is difficult to assign the precise mechanism from 
the available experimental material it apparently 
is possible to understand these results in terms 
of ‘‘crossing”’ potential surfaces. ; 

The following general considerations regarding 
chemiluminescence in flames may be mentioned. 
Figs. 1 and 2 show only a single upper potential 
energy surface. For complicated molecules there 
will be many such surfaces and if the activated 
complex on any such surface is denoted by 
A,*(n), the concentration of A,*(m) is deter- 

16 R, S. Anderson, J. Am. Chem. Soc. 59, 2115 (1937). 
17 A, Hollaender and W. D. Claus, J. Opt. Soc. Am. 25, 


270 (1935). 
18 Audubert and Victorin, J. Chim. Phys. 34, 18 (1937). 
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mined by the temperature and by the free energy 
difference between A; and A,*(m) and for 
sufficiently high temperatures many of these 
upper levels will contribute chemiluminescent 
reactions. On the other hand the frequency of 
the emitted light is determined by the energy 
difference between B,(m) and B,;. Here A,*(n) 
and B,(n) are denoted by simply A,* and B,, on 
the lowest excited surface represented in Figs. 1 
and 2. Since reactions proceed with a free energy 
decrease, the free energy of B; will be lower than 
that of A;. We thus expect that in general the 
energy of the light emitted will be greater than 
the free energy of activation of A, going to 
A,,*(n), but that the amount of radiation of the 
corresponding frequency will, in the absence of 
quenching processes, be governed by the energy 
of activation. Quenching processes may, how- 
ever, destroy any correspondence of this kind, 
especially if the excited molecule, B,(m), has a 
long mean life. The black-body spectrum corre- 
sponding to the temperature of the flame will, 
of course, form a background on which is super- 
imposed the radiation from chemical reactions. 

In a recent paper dealing with the interpreta- 
tion of alkali halide-thallium phosphors, Dr. 
Frederick Seitz'® has concluded that it is possible 
for a particular phosphor containing an excited 
electron to pass less readily over a potential 
barrier by an adiabatic process than it does over 
a higher barrier which can only be reached by 
absorption of radiation. This case, if correctly 
interpreted, is of great interest in providing the 
only known type of reaction in which the rate 
by the radiation mechanism is actually faster 
than the rate resulting from a collision mecha- 
nism involving a lower energy of activation. 

We want to thank Doctor Rubert S. Ander- 
son, Professor M. Polanyi and Professor Hugh 
S. Taylor for helpful discussions while this 
work was in progress. One of us (H.E.) is glad 
to have this opportunity of expressing his 
appreciation to the Rockefeller Foundation for 


the opportunity of working in England. 


19 F, Seitz, J. Chem. Phys. 6, 150 (1938). 
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The thermal and photochemical decompositions of mercury dimethyl have been studied in 
presence and absence of hydrogen. The thermal process occurs measurably above 290° in 
absence of hydrogen and is complex, yielding ethane, methane and carbonaceous products. 
In presence of hydrogen there is no carbonaceous deposit and methane is predominantly 
formed with consumption of hydrogen, which accelerates the decomposition. The influence 
of surface effects has been noted. The photo-decomposition has been studied from 50-300°C. 
In absence of hydrogen, ethane is practically exclusively formed up to 200°. In presence of 
hydrogen, methane is formed at all temperatures, in increasing amounts with increase of 
temperature. A mechanism based on photo-decomposition to free methyl radicals with asso- 
ciation to yield ethane and reaction with hydrogen to yield methane has been quantitatively 
determined. For the activation energy CH;+H2=CH,+H a value of 9+2 kcal. has been 


deduced. 





HE reaction CH3+H2—-CH,+H has re- 

cently been shown to be of importance as 
an intermediate process in exchange reactions 
involving the simpler saturated hydrocarbons," 2 
and any independently derived information re- 
garding its energetics would be of value in 
explaining the results obtained. The present 
investigation was undertaken in an attempt to 
provide such information, the metal alkyl, mer- 
cury dimethyl, being used as a source of free 
methyl radicals. It was first necessary to establish 
the conditions under which the decomposition of 
mercury dimethyl proceeded at a measurable 
rate: then, from the effect of added hydrogen 
upon the kinetics of the decomposition, some 
conclusions were drawn regarding the process 
under consideration. 

Direct evidence of the momentary existence of 
free organic radicals in the decomposition of 
metal alkyls was first obtained by Paneth and 
Hofeditz,? who passed lead tetramethyl vapor 
through a quartz tube heated at one point to 
about 700°, and detected free methyl! radicals by 
causing them to combine with a number of 
metals deposited as mirrors in the cool parts of 
the tube. It was found‘ that free ethyl could be 
prepared from lead tetraethyl and identified in a 
similar manner, but the higher analogues were 


on Morikawa and Taylor, J. Chem. Phys. 5, 203 
/). 


—_— Benedict and Taylor, J. Chem. Phys. 5, 212 
7). 


(1 
* Paneth and Hofeditz, Ber. 62, 1335 (1929). 
‘Paneth and Lautsch, Ber. 64, 2702 (1931). 


apparently unstable at the temperatures neces- 
sary for the decomposition of the corresponding 
metal alkyls. Indirect evidence of the thermal 
liberation of free radicals from mercury dimethyl 
has been furnished by Taylor and Jones® who 
showed that such postulated radicals induced the 
polymerization of ethylene at reaction tempera- 
tures of 200-300°. Frey® also observed that the 
addition of 1 percent of the compound to butane 
brought about the decomposition of twenty 
equivalents at 525°, at which temperature the 
hydrocarbon itself would not have been affected 
in the time of the experiment ; this result can only 
be explained as due to chain initiation by free 
methyl radicals. Rice, Rodowskas and Lewis’ 
obtained similar results with acetone containing 1 
percent of mercury dimethyl in the temperature 
range 350-400°: the relatively short chains 
observed were attributed to more effective chain- 
breaking processes than in the case of butane. 
The photochemical decomposition of mercury 
dimethyl has been somewhat more extensively 
investigated, notably by Terenin and co- 
workers. The absorption spectrum shows a 
continuum beginning at 2800A and extending 
further into the ultraviolet, upon which is 
superimposed a group of diffuse bands in the 


5 Taylor and Jones, J. Am. Chem. Soc. 52, 1111 (1930). 

6 Frey, Ind. Eng. Chem. 26, 200 (1934). 

7 Rice, Rodowskas and Lewis, J. Am. Chem. Soc. 56, 
2497 (1934). 

8 Terenin, J. Chem. Phys. 2, 441 (1934); Acta Physico- 
chim. U.S. S. R. 1, 762 (1935); Terenin and Prileshajeva, 
Trans. Faraday Soc. 31, 1483 (1935). 
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region 2100—-1970A. Terenin postulated a mecha- 
nism of decomposition involving free methyl 
radicals, and reported their identification by the 
metal mirror technique, the maximum concen- 
tration of radicals being observed when light in 
the 2100A region was used. Linnett and Thomp- 
son,® using radiation of about 2537A, reported 
rapid conversion of mercury dimethyl, largely to 
mercury and ethane, at room temperature, the 
quantum yield of the process being very nearly 
unity ; these investigators first concluded that at 
least 90 percent of the metal alkyl decomposed 
directly to molecular ethane in the primary act, 
only a small fraction decomposing by way of free 
radicals. In their second paper they reversed 
their position and developed a mechanism of 
decomposition via free methyl radicals which 
leads to quantum yields higher than unity at 
elevated temperatures. In studying the analogous 
decomposition of lead tetramethyl, Leighton and 
Mortensen!® held the radical view; in this case, 
free methyls were detected by the use of lead 
mirrors, the method being rendered much more 
sensitive than heretofore by the introduction of a 
radioactive indicator. 

In several instances, the gaseous products 
resulting from the decomposition of metal alkyls 
have been analyzed, and the intermediate forma- 
tion of free radicals has been suggested by the 
results obtained. Simons, McNamee and Hurd" 
and Paneth, Hofeditz and Wunsch,” postulate 
an initial dissociation of lead tetramethyl into 
free methyls in order to account for such 
products as methane, iso-butylene, and a solid 
carbonaceous deposit. Meinert!* was led to a 
similar conclusion on examining the products 
from the thermal decomposition of lead tetraethyl. 

The effect of added hydrogen in the decompo- 
sition of metal alkyls has received only brief 
consideration. It was disregarded by Paneth and 
Hofeditz in their original investigation,’ in which 
hydrogen Was used as a carrier gas, but in their 
later publication” interaction between free 


®Linnett and Thompson, Trans. Faraday Soc. 33, 501, 
874 (1937). 

10 Leighton and Mortensen, J. Am. Chem. Soc. 58, 448 
(1936). 

ul _— McNamee and Hurd, J. Phys. Chem. 36, 939 
(1932). 

% Paneth, Hofeditz and Wunsch, J. Chem. Soc. 372 
(1935). 

13 Meinert, J. Am. Chem. Soc. 55, 979 (1933). 
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methyls and molecular hydrogen was postulated 
to account for the removal as methane of a large 
fraction of the radicals produced at 900°, and 
subsequently passed through a tube at 20°. It is 
significant that the relative extent of methane 
production was greatly increased by raising the 
temperature of this tube to 350°. Simons, 
McNamee and Hurd" mention that the addition 
of hydrogen to lead tetramethyl in a single 
experiment at 700° increased the proportion of 
methane produced at the expense of the ethane, 
and suggest the same explanation. 

Hartel and Polanyi" have studied the reaction 
CH;+H.—CH.+H by the method of “highly 
dilute flames,”’ free methyl radicals being liber- 
ated in the presence of hydrogen by the action of 
sodium atoms on methyl halides. From observa- 
tions of the rate of reaction at two temperatures, 
they deduced an activation energy of approxi- 
mately 8 kcal./mole for the reaction. Trenner, 
Morikawa and Taylor,' from their results on the 
reaction between ethane and atomic hydrogen 
produced by Wood’s method, have suggested 
11 kcal./mole as an upper limit. 


EXPERIMENTAL DETAILS 


Apparatus 


Preliminary experiments showed that the de- 
composition of mercury dimethy], alone or in the 
presence of hydrogen, did not involve any change 
in the total pressure, so the course of the reactions 
could only be followed by analysis of the products. 
For reasons outlined below, it was necessary to 
operate with small samples of gas or vapor, so a 
special apparatus was constructed. The essential 
parts are represented diagrammatically in Fig. 1. 

Pyrex glass was used throughout, capillary 
tubing being inserted wherever possible to limit 
the dead space. V; and V2 were all-metal valves 
with copper sylphon bellows, their use being 
necessitated by the readiness with which mercury 
dimethyl dissolves in grease of any kind, and the 
experimental procedure was so arranged that 
these were the only stopcocks with which the 
vapor came into more than fugitive contact. The 
mercury dimethyl was stored in A, and because a 
small amount of Picein was used in the ground 


44 Hartel and Polyani, Zeits. f. physik. Chemie B11, 97 
(1930). 
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Fic. 1. Apparatus for study of decomposition of mercury 
dimethyl. 


glass joints between the valves and the tubing, it 
was kept frozen out by means of a carbon 
dioxide-toluene mixture except for short intervals 
of time when the samples were being taken. The 
reaction vessel B was of Pyrex for the thermal 
experiments, and of quartz with quartz-to-Pyrex 
graded seals for the photochemical ; a number of 
independently calibrated vessels were used, those 
of Pyrex ranging from 220 to 250 cc, and those of 
quartz from 115 to 150 cc in volume. 

Attempts were first made to heat the sample in 
vessel B in situ, but this was soon found to be 
impracticable due to the readiness with which the 
mercury dimethyl distilled into the cooler parts of 
the connecting tubing and manometer, even with 
a small temperature gradient. In all the later 
runs, consequently, the reaction vessel was sealed 
off after introduction of the sample. After 
heating or irradiating for a definite length of 
time, the vessel could then be sealed in place by 
means of the 10 mm tube shown to the right of 
the vessel in the diagram. After evacuation, the 
small glass tip could be broken at will electro- 
magnetically with an iron plunger, and the 
reaction products thus readmitted to the 
apparatus. 

In the thermal runs, vapor baths containing 
the following substances were used: aniline, 
naphthalene, diphenyl, alpha-naphthol, dipheny]- 
amine, and anthracene. For the photochemical 
experiments, a furnace was constructed out of a 
cylindrical block of aluminum about twice the 
length of the reaction vessel, and of such diameter 
that a hole drilled to fit the reaction vessel left 
walls of 1” thickness. An aperture the length of 


the quartz vessel and 3” in width was cut in one 
side, covered by a quartz window, and the 
furnace was wound in such a fashion that the 
heating element did not cross this window. [he 
furnace was then well covered with asbestos 
paste, and it was found that, at the highest 
temperatures used, the temperature did not vary 
by more than 5° throughout. 

The arc used was of the mercury-in-quartz 
type, operated ‘‘hot”’ at 50 volts and 4 amp. d.c., 
under which conditions no resonance radiation 
could be detected; it was operated for fifteen 
minutes preceding every run in order to attain 
reproducible intensity, and then placed in a 
vertical position against the quartz window in the 
furnace for the duration of the experiment. 

Hand Tépler pumps C and E of approximately 
400 and 250 cc capacity, respectively, were used 
to transfer gas samples from one portion of the 
apparatus to another. Into a similar vessel D 
were sealed at different levels three glass pointers ; 
by raising the confining liquid to any of these 
tips, and at the same time setting the level in the 
connecting arm of manometer M, at a standard 
reading, it was possible to confine a gas sample in 
a definite volume; these volumes were calibrated 
with mercury and found to be 33.6, 104.6, and 
204.6 cc. Mercury was also used as the confining 
and manometric fluid throughout. 

The bulb G of about 50 cc capacity, contained 
between 1.5 and 2 g of an active copper hydro- 
genation catalyst, and the U-tube at F contained 
about 5 g of partially-reduced Kahlbaum copper 
oxide. 

The open arms of the manometers were actu- 
ally connected to the vacuum line (not shown), 
thereby eliminating the effect of barometric 
variations; each manometer and hand Topler 
pump was provided with a leveling bulb con- 
nected through an air-trap as shown for C, so 
that all gas pressures could be read at constant 
volume. 


Materials 


The mercury dimethyl used was prepared by 
the method of Marvel and Gould,'® and further 
purified by repeated fractionation at low pres- 
sures and in the presence of freshly fused calcium 


18 Marvel and Gould, J. Am. Chem. Soc. 44, 153 (1922). 
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chloride. The hydrogen was obtained by the 
electrolysis at 2-3 amp. of 10 percent potassium 
hydroxide solution, and passed in succession 
through calcium chloride, activated charcoal, 
platinized asbestos at about 250°, and a trap 
chilled in liquid air. It was not considered neces- 
sary to purify the tank helium used to increase 
the total pressure in a few experiments. The 
active copper was prepared from Kahlbaum 
copper oxide by reduction and oxidation at 
successively lower temperatures, the last reduc- 
tion being made below 150°; tests on tank 
ethylene showed rapid hydrogenation below 100°. 


EXPERIMENTAL PROCEDURE 


Before each experiment, the volume of the 
particular reaction vessel used, together with the 
connecting tubing to a set level of the mercury in 
manometer M, was determined by measuring the 
pressure of a sample of air and then transferring 
it, by means of C, to the burette D, where its 
pressure was again measured in a previously 
calibrated volume. This was usually repeated at 
the close of the experiment, to determine the 
increase in volume due to the introduction of the 
10 mm tubing containing the “‘inner seal.” 

After the sample of air had been pumped out 
of B, a sample of mercury dimethyl was admitted 
from A. Due to the readiness with which the 
material distills and condenses, it was considered 
advisable to set as an upper limit of mercury 
dimethyl pressure its vapor pressure at room 
temperature, which is in the vicinity of 50 mm; 
in very few cases did the pressure of the metal 
alkyl significantly exceed 40 mm. With valve V2 
closed and V; opened very slightly, the tempera- 
ture of the vessel A was carefully raised to within 
a few degrees of that of the room, and the vapor 
thus allowed to diffuse slowly into B. After 
standing for some time in a water bath at room 
temperature, the pressure and temperature of the 
sample in B were recorded. A sample of hydrogen 
could similarly be admitted through V2, the 
pressure of the hydrogen being maintained high 
enough to prevent back diffusion of the mercury 
dimethyl. 

To avoid loss or decomposition of the mercury 
dimethyl while sealing the reaction vessel, it was 
first immersed in liquid air ; this expedient did not 
prevent a small loss of hydrogen, but this could 
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be estimated by comparing the total pressures 
before and after the run, as the reaction proceeds 
without change in the total pressure: moreover, 
as the hydrogen was always in large excess, the 
importance of such an error would be minor. The 
disconnected vessel was then transferred to the 
vapor bath or furnace, depending on the experi- 
ment. In the former, the time of heating was 
measured from the introduction of the vessel, the 
time required for it to reach the bath temperature 
being considered to be offset by that necessary 
for cooling after its withdrawal; in the photo- 
chemical runs, however, the vessel was placed in 
the furnace for fifteen minutes prior to illumi- 
nation ; it will be seen that the thermal reaction 
during that period of time at the temperatures 
used is negligible. 

After heating or illumination, the reaction 
vessel was again connected to the apparatus and 
the product gases readmitted. 


Analysis 


The vessel B was immersed in liquid air, and 
the fraction gaseous at that temperature pumped 
over, by repeated use of the hand Tépler pump C, 
into D where its pressure and temperature at 
known volume were recorded. This gas fraction 
could consist only of the so-called ‘‘permanent 
gases” methane and hydrogen. The possibility of 
imperfect fractionation due to solution of these 
gases in other products in B was investigated by 
closing V2, melting and refreezing these products, 
and again pumping any residual gas over into D: 
in no case, however, was there observed to be any 
gas evolved. The hydrogen was removed from the 
mixture by passing the gases repeatedly over the 
copper oxide in F at 270—300°, with the mercury 
leveling bulbs attached to D and E;; the residual 
gas after this operation was measured as methane. 

A second gas fraction was transferred from B 
to D in the same manner, the vessel B being 
immersed in a carbon dioxide-toluene mixture. 
This could consist of ethane, propane, and 
possibly butane, together with the lower members 
of the unsaturated series. Numerous attempts 
were made to detect unsaturated hydrocarbons 
by passing this second gas fraction, together with 
added hydrogen, over the active copper in vessel 
G maintained at 100°, but in no case was there 
any evidence of the pressure decrease which 
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would attend hydrogenation. Moreover, the ob- 
served stoichiometry of the decomposition in the 
majority of the experiments indicated that 
propane and butane, if present at all, appeared in 
negligible amounts, so that the second gas 
fraction might well be designated as ethane. 

The residue in B, which probably consisted 
entirely of undecomposed mercury dimethyl in 
the majority of the runs, was not transferred 
from the reaction vessel, but was estimated by its 
pressure measured on manometer M, while the 
vessel was immersed in a water bath at room 
temperature. By using this value in conjunction 
with the directly determined quantities of gase- 
ous products, a satisfactory carbon and hydrogen 
balance could normally be obtained. 


EXPERIMENTAL RESULTS 


Thermal decomposition of mercury dimethyl 


No measurable decomposition of the metal 
alkyl alone could be detected after long heating 
at 184, 215 and 255°. At 290°, gaseous products 
corresponding to about 2 percent decomposition 
were observed after five hours heating; at 302°, 
the rate had increased only very slightly, but 
when the temperature was raised to 348°, it was 
relatively rapid. 

The course of the decomposition at 348° was 
followed by heating separate samples for definite 
periods of time, and the results have been 
tabulated in Table 1. All the pressures given have 
been calculated to 0°, only a short extrapolation 
being involved, as all the actual measurements 


TaBLe I. Thermal decomposition of Hg(CHs)2 in absence 





























of 2. 
INITIAL 

Temp. TIME Hg(CHs)2 CH, C2He 
s RuN HR. MM MM MM 
302 25 3 21.3 -— 1.0 
8 21.6 1.0 1.4 

213 Za2 2.4 3.2 

30 22.8 5.1 2.9 

413 23.3 7.7 2.9 

| 

348 106 2 37.0 14.2 4.8 
; 107 4 38.6 22.5 4.0 

| 105 6 43.3 29.0 | 5.1 

| 102 11 47.7 39.1 8.3 

| 103 21 47.1 39.8 7.3 

| 108 4 14.8 4.5 2.9 

| 1091 4 46.7 44.3 8.4 








' Surface increased 8-fold by packing with Pyrex tubing. 





were made at room temperature. The pressures 
for the products listed in the last two columns 
have been recalculated to the volume of the 
original reaction vessel in each case, for com- 
parison with the values in the fourth column. 
Appended to the table are two additional runs at 
348°, which serve to throw some light on the 
nature of the decomposition: run 108 shows the 
effect of decreasing the initial mercury dimethy] 
pressure, run 109 that of increasing the total 
surface. 

A single run at 302° is also reported for 
comparison, although the accuracy of the data is 
relatively low. This experiment was performed 
without detaching the reaction vessel from the 
apparatus; rough analyses were performed by 
interrupting the heating, and measuring the 
total pressure when the vessel was immersed in a 
water bath, a CO, ether mixture, and liquid air. 
The pressures given, therefore, are obtained by 
fairly long temperature extrapolations, and serve 
only as a rough indication of the effect of 
temperature. It is definitely apparent that the 
process producing methane shows much greater 
dependence upon the temperature than does that 
forming ethane. 

In every case in which measurable decompo- 
sition took place, it was observed that a dark 
brown deposit formed on the walls of the reaction 
vessel. This could be removed to a greater or less 
extent by washing out with cold water, on the 
surface of which it formed an oily scum. On 
gentle heating in the presence of air, the deposit 
darkened in color, and became definitely carbo- 
naceous. Reference to runs 102 and 103 in Table 1 
shows that although decomposition is un- 
doubtedly complete, less than 60 percent of the 
carbon has been found in the product gases; it 
would appear that for every carbon atom ap- 
pearing as methane there is one lost in the solid 
deposit on the walls. Added evidence for the 
belief that the deposit might well be a CH» 
polymer is supplied by the experimental fact that 
no free hydrogen was detected in any run. 


Thermal decomposition of mercury dimethyl in 
presence of hydrogen 


In the runs listed in Table II, the initial 
mercury dimethyl pressures were 40.0+0.2 mm, 
and those of hydrogen (corrected for loss in 
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sealing the reaction vessel) 200+4 mm. As in 
Table I, all pressures are given for 0° and the 
volume of the reaction vessel. The values in the 
fourth and fifth columns were obtained by a 
method of difference, and so are less accurate 
than those for methane and ethane, which have 
been measured directly and in smaller volumes. 
However, the stoichiometry of the reaction pro- 
ducing methane is apparent, one mole of mercury 
dimethyl and one of hydrogen disappearing in 
the production of two moles of methane. This 
indicates that, with the introduction of hydro- 
gen, a new methane-producing process occurs. 
The absence of any dark deposit at the con- 
clusion of an experiment when hydrogen is 
present is corroborative evidence. 

The presence of hydrogen markedly accelerates 
the decomposition, which is complete in one hour 
at 348° instead of eleven, as in the absence of 
hydrogen. That this effect is not due merely to 
the increase in the total pressure is shown by the 
runs using helium. In addition, the methane- 
producing process in the presence of hydrogen is 
much less strongly temperature-dependent than 
that in the case of mercury dimethyl alone, 
indicating a much lower energy of activation. 


Effect of initial pressure 


The data of Table III show the influence of 
reactant pressures at two temperatures. They 
yield the rate expressions: 
at 302°: d(CH,)/dt=k(Hg(CHs)2)**(He)*’, 
at 348°: d(CH,)/dt=k(Hg(CHs)2)"7(H2)”’. 


TABLE II. Thermal decomposition of Hg(CHs)2 in presence 











0 2. 
He 

TEMP. TIME Hg(CHs)2 USED CH, CoHe 
~~ RuN HR. USED MM MM MM MM 
302 125 4 9.3 8.8 | 13.2 — 
126 i 8.1 6.4 | 12.0 | — 

124 1 13.8 11.2 | 22.0 — 

122 2 23.2 23.6 | 42.8 —_ 

123 4 40.1 44.8 | 80.2 0.7 

130! 43 40.0 —- 2.4 — 

348 115 3 35.3 31.2 | 65.0 1.2 
114 1 39.8 39.8 | 79.3 2.4 

119 1 39.6 38.4 | 74.9 Be 

112 2 37.8 37.2 | 72.0 1.6 

110 4 40.0 40.0 | 77.6 2.4 

118! 3 40.0 — 18.3 1.2 


























1200 mm helium substituted for hydrogen. 


P. CUNNINGHAM AND H. S. 





TAYLOR 


TABLE III. Variation of initial pressures. 











INITIAL INITIAL 

TEMP. TIME Hg(CHs)2 He CH, 

i RUN HR. MM MM MM 
302 122 2 40.1 199.5 42.8 
127 40.4 100.6 26.9 
129 18.7 185.5 36.6 

128 19.6 98.4 22.2 
348 115 3 39.9 198.2 65.0 
120 40.0 100.4 40.1 
121 19.7 195.5 38.6 
117 19.7 98.4 24.0 























There is a strong temperature dependence with 
respect to mercury dimethyl but none with 
respect to hydrogen. 

Experiments at lower temperatures showed 
erratic behavior indicating surface effects, meth- 
ane formation increasing with increase of mercury 
dimethyl, the effect of hydrogen remaining 
unchanged. 


Photochemical decomposition of mercury di- 
methyl 


The data of Table IV record the results of 
photodecomposition by the unfiltered radiation 
of a hot mercury arc from 50 to 300°, with an 
initial pressure of 40 mm of mercury dimethyl. In 
no photochemical experiment, with or without 
hydrogen, was any deposit observed to form on 
the walls of the reaction vessel, with the possible 
exception of a thin gray film of mercury. No 
methane was observed in the decomposition of 
mercury dimethyl below 160°, and it seemed 
probable that the reaction was free of the 
undesirable heterogeneous characteristics of the 
thermal decomposition. 


TABLE IV. Photochemical decomposition of Hg(CHs)2 im 
absence of He. 











Hg(CHs)2 
TEMP. TIME DECOMP. CHa CoHe CoH¢/MIN. 

<— RuN MIN. MM MM MM MM 
50 | 142 30 11.2 — 10.4 0.35 
144 60 20.4 — 19.7 0.33 

146 90 28.7 — 28.0 0.31 

90 | 152 60 23.2 — 24.6 0.41 
160 | 154 30 20.1 0.1 | 20.8 0.69 
200 | 168 10 9.3 0.3 7.9 0.79 
250 | 173 6 6.4 0.6 6.2 1.03 
275 182 5 5.7 _ 4.1 0.82 
300 | 184 3 5.1 1.2 3.4 1.13 
185 3 — — 2.8 0.93 
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TaBLE V. Photochemical decomposition of Hg(CHs)2 in 
presence of Ho. 











Hg(CHs)2 3 CHa 

TEMP. TIME | DECOMP. | CHa/CeHe| +CoHe | CH4/min.|CoHe/MIN. 
°C | RUN] MIN. MM MM | MM MM MM MM 
50 | 143] 30 13.8 4.0 |10.4| 12.4 0.13 0.35 
145 60 23.2 6.8 |19.5 22.9 0.11 0.33 
147} 90 (28.7) {11.1]29.3 | 34.9 0.12 0.33 
90 | 153] 60 33.4 {10.8 /30.2| 35.6 0.18 0.50 
160 | 155} 30 21.3 8.9 116.1] 20.6 0.30 0.54 
159 28.1 13.6 |23.8| 30.6 0.45 0.79 
160 28.6 |16.2|20.9| 29.0 0.54 0.70 
200 | 169} 10 14.2 10.9 | 7.8 13.3 1.09 0.78 
250 | 174 6 11.1 12.1} 4.1 10.3 2.02 0.70 
275 | 183 5 11.8 11.8] 5.0 10.9 2.36 1.00 
300 | 180 3 10.1 12.4} 3.2 9.4 4.13 1.07 
181 3 9.9 12.2] 2.8 8.9 4.07 0.93 
































Photochemical decomposition of mercury di- 
methyl in presence of hydrogen 


Addition of hydrogen to the mercury dimethy] 
was found to result in the production of measur- 
able amounts of methane even at the lowest 
temperature employed. It was observed, more- 
over, that the process producing methane showed 
a much greater dependence upon temperature 
than that forming ethane, indicating a process of 
higher activation energy. The rate of ethane 
production was not affected by the presence of 
the hydrogen, indicating that the processes 
producing methane and ethane proceed quite 
independently of each other. This observation 
must be given due importance in postulating a 
mechanism for the decomposition. 

In all runs reported in Table V very nearly 
200 mm of hydrogen has been added to the 
mercury dimethyl. The values in the seventh 
column have been given for comparison with 
those of the fourth column, to indicate how 
satisfactory is the carbon balance throughout. It 
will be observed also from Table IV that the 
mercury dimethyl decomposed in each run is 
accounted for within the experimental error. 

Special tests, using quartz vessels containing 
mercury and hydrogen as filters, showed that the 
effects produced were not due to mercury 
resonance radiation. It was also shown that a 
mixture of He+Dz did not react to form 2HD 
under the experimental conditions. When a 
“cold” arc was used, the character of the decom- 
position changed, because of the resonance 
radiation, methane predominating even at room 
temperature, because of interaction of the mer- 
cury dimethyl and atomic hydrogen produced by 
photosensitisation with excited mercury. 
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DISCUSSION OF RESULTS 


The primary photo-process, in agreement with 
Terenin® and the later conclusions of Linnett and 
Thompson,® must occur in either of the two 
following ways 

HgCH;+CHs; 


Hg(CHs3)2 


\ 
Hg+2CHs(h). 


The preponderance of methane in presence of 
hydrogen at the highest temperatures is definite 
evidence that dissociation to yield mercury and 
molecular ethane is a minor or negligible result of 
the absorption act. A choice between the two 
radical mechanisms cannot be made from the 
kinetics of the over-all reaction. The complex 
HgCH; will, in any case, be of relatively short 
life. 

The ethane formation which preponderates in 
the low temperature range may be ascribed to the 
recombination process 


CH;+CH3=CeH¢(k,), 


the rate of production of ethane in absence of 
hydrogen increases by a factor of 3 in the range 
50 to 300° (Table IV). There is no evidence that 
this can be ascribed to temperature dependence 
of light absorption, since this is very strong, even 
at room temperature. The temperature coefficient 
indicates, therefore, an increase of quantum yield 
and therefore a chain mechanism, for example, 


CH;+Hg(CHs)2=CoHe+Hg+CuHs, 


the chain-breaking mechanism being the recom- 
bination of two methyl groups. This latter 
reaction must have a small activation energy 
since it occurs less rapidly than the recombination 
of CH;+I in the photolysis of methyl iodide as 
indicated by the low quantum efficiency of the 
iodide decomposition and by the thermal stability © 
of alkyl iodides.'® The small temperature coeff- 
cient of ethane formation must be composite of 
these two processes. 
The data of Tables IV and V show that the 
ethane producing reactions are not measurably 
affected by the presence of a large excess of 


16 Ogg, J. Am. Chem. Soc. 56, 533 (1934). 
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hydrogen. However, the over-all rate of de- 
composition of mercury dimethy] is progressively 
increased, by a 5-fold addition of hydrogen, and 
approximately in the ratio 1 to 8 between 50 and 
300°. An increasingly large proportion of methane 
appears in the products with increase of tempera- 
ture. We ascribe this to the following new 
secondary steps 


CH3+He=CH4+H (ke) 
and 
H+Hg(CHs)2=CHi+Hg+CHs(ks), 


this chain being terminated by the recombina- 
tion of two methyl radicals (R,). 
For a photochemical stationary state 


d[CHs ]/dt=0=kyJavs— Rol CHs |[ He } 
+k3[H J[ Hg(CHs)2]—Rs[CHs ? 
and 


d(H ]/dt=0=k.[CHs Jl He ]—AsLH JLHg(CHs) |. 
By subtraction, 
Rilavs=kaLCHs ? 
or [CHs3 ]= (Rilabs/Rs)?. 
Substituting these values in the rate expression, 


d[CH,]/dt=ko[CHs ][H2]+s(H ][THg(CHs)2] 
=2k-[CH; ][H2] 
= 2ko( Ril avs/ks)*L He |. 


If the light absorption is constant, the rate of 
methane production is independent of mercury 
dimethyl concentration. The data at 50° in 
Table V are in agreement with this result. In this 
formulation of reaction mechanism it is seen that 
ethane formation is also independent of the 
hydrogen concentration as observed. 

The production of methane in the absence of 
hydrogen must be ascribed to processes of the 


types 
CH;+Hg(CHs)2 
CH.+HgCH; 


=CH,+CH:2.Hg.CHs 


\ 
CH.+Hg+CH; 


and CH;+CH;=CH,+CHsz. 





P. CUNNINGHAM 


AND H. S. TAYLOR 





The large temperature coefficient of thermal 
reaction between 302 and 348° cannot be 
attributed to the reaction CH;+Hpe since this 
has, as will be shown, a much smaller coefficient. 
The dark deposit on the reaction walls is evidence 
of the production of polymethylenic polymeriza- 
tion products. Evans and Lee!’ on studies of the 
electrolysis of methyl mercuric bromide and 
Shukla and Walker’ in electrolysis of fatty acids 
and their salts postulate methane production by 
the second of the above processes. 


Activation energies of the reactions 


If the logarithm of the rates of methane 
production from Table V be plotted against 1/T 
in the range 160-300° a straight line is obtained 
whose slope is 8.1 kcal. per mole. As is evident 
from the expression 


d[ CH, ]/dt=2ko[fCHs |[He | 


this value is a summation of the effect of temper- 
ature on the stationary state concentration of 
methyl radicals and the activation energy Ey, of 
the process CH;+H2=CH.+H. 

A plot of ethane production, in presence or 
absence of hydrogen, against 1/T gives a single 
straight line over the entire temperature range 
the slope of which indicates an activation energy 
of 1.5 kcal. per mole. As suggested previously, 
this energy requirement may be attributable to 
the reaction between free methyl and mercury 
dimethyl. In such case the effect of temperature 
on the recombination of methyl radicals whether 
bimolecularly, or at the wall, must be very small 
indeed and Ey would be 8.1 kcal. per mole. This 
latter value would be reduced, if the value of Ez 
for ethane formation were to be attributed to an 
activation energy of the recombination process, 
by the whole amount, 1.5 kcal., in case of a wall 
recombination, or by 0.75 kcal. for a bimolecular 
association in the gas phase. These alternatives 
would set the value for Ey between 6.6 and 8.1 
kcal. with the greater weight attaching to the 
upper figure. 

The experimental figures of Hartel and 
Polanyi led them to a value Ey=8 kcal., 
although a series at 10 mm hydrogen pressure 


17 Evans and Lee, J. Am. Chem. Soc. 56, 654 (1934). 
18 Shukla and Walker, Trans. Faraday Soc. 27, 35 (1931). 
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gave a value of ~6 kcal. In a parallel investi- 
gation by Taylor and Rosenblum'® on _ the 
photolysis of acetone in presence of hydrogen the 
value for Ex, deduced on the basis of bimolecular 


19H. S. Taylor and C. Rosenblum, J. Chem, Phys. 6, 119 
(1938). 
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recombination with negligible activation energy 
is 8.7 kcal., rising to 11 kcal. if the recombination 
process were unimolecular at the wall. Upon these 
bases therefore we must conclude that the most 
trustworthy value now available for the activa- 
tion energy, Em, is 9+2 kcal. 
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An Improved Calculation of the Energies of Metallic Li and Na 


J. BARDEEN* 
Harvard University, Cambridge, Massachusetts 


(Received April 13, 1938) 


The method of calculation used in the present work is 
a modification of that used by Wigner and Seitz in their 
original calculations of the energies of metallic Li and Na. 
The main difference lies in an improved method for the 
calculation of the Fermi energy. In each polyhedral cell, 
the wave function of an electron is taken to be of the form 
ve=([uo(r) +7(k-x)(v1(r) —uo(r)) ] exp [ik-x], where 1o(r) 
and v;(r) are radial s and p functions, respectively, which 
depend only on the distance from the ion in the center of 


HE most accurate method for the calculation 

of metallic binding which has yet been 
suggested is that of Wigner and Seitz.! The 
method is especially suitable for the monovalent 
metals, and was first applied to the calculation of 
the binding energy and lattice constant of 
sodium. A rough value was also given for the 
compressibility. Lithium was considered in a 
later paper by Seitz.? His calculation of the 
energy did not extend to values of the lattice 
constant smaller than the observed, so that it 
was not possible to obtain the compressibility. 
The main object of the present work is to 
determine the energies of Li and Na and their 
dependence on volume with sufficient accuracy to 
obtain the compressibilities and their variation 
with pressure. A comparison of the results with 
recent experimental values of Bridgman, which 


* Society of Fellows. 

‘E. Wigner and F. Seitz, 43, 804 (1933); 46, 509 (1934); 
E. Wigner, Phys. Rev. 46, 1002 (1934). A survey of this 
work, together with references to applications of the method 
to other problems may be found in Mott and Jones, Theory 
of the Properties of Metals and Alloys (Oxford, 1936), 
Chapters II and IV. 

*F. Seitz, Phys. Rev. 47, 400 (1935). 


the cell. Both functions are determined explicitly from 
the differential equations, and the energy of the electron 
is expressed in terms of the boundary values of the wave 
function. Values of the Fermi energy, and of the total 
energy of each metal, are tabulated as a function of the 
lattice spacing. Calculated values of the lattice constants, 
heats of sublimation, and compressibilities are in fair 
agreement with experiment. 


extend to pressures of the order of 45,000 kg/cm?, 
will be given in the following paper.* 

The two main factors which determine the 
binding energy of a metal are the energy of an 
electron in its lowest state and the Fermi energy. 
Only two electrons can occupy the lowest state ; 
the Fermi energy represents the average addi- 
tional energy of electrons in higher states of the 
filled band. We will first give the general theory 
of the method used for the calculation of the 
electronic wave functions and energies, and then 
give the detailed numerical calculations of the 
binding energies. 


ELECTRONIC WAVE FUNCTIONS AND ENERGIES 


The wave function of an electron in its lowest 
state is determined as follows. By passing planes 
perpendicular to and bisecting the lines joining 
near neighbors, the whole space of the crystal is 
broken up into a set of regular polyhedra, one 
surrounding each ion. It is a good approximation 


3 A discussion of an approximate method for the calcu- 
lation of the energies of the alkali metals, due to Frohlich, 
is also given in this paper. 
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to take for the potential in each polyhedron that 
of the ion in the center of the polyhedron.‘ Since 
the wave function of an electron in its lowest 
state has the symmetry of the lattice, its normal 
derivative must vanish on the surface of the 
polyhedron. Wigner and Seitz replace the 
polyhedron by a sphere of equal volume (s 
sphere), and satisfy the boundary condition on 
the surface of this sphere. The wave function is 
then the radial solution of the Schrédinger 
equation : 


— (h?/2p) Apot V(r) Yo= Evo, (1) 
subject to the boundary condition 
(do/dr) rer, =0, (2) 


where 7, is the radius of the s sphere. In Eq. (2), 
V(r) is the potential of a free (positive) ion, and 
E¢ is the energy of the electron in its lowest state. 

The wave function of an electron in a higher 
state is of the form 


Ve=Ux(x) exp [ik-x], (3) 


where u;(x) has the translational period of the 


lattice. The equation for u,, obtained by substi- 
tuting the above in the Schrédinger equation, is: 


— (h?/2u) (Au, +2tk- grad ux) 
+ Vu,= (E; —hPk?/2y) ux, (4) 


which differs from the equation for uo (the wave 
function of an electron in its lowest state) in the 
presence of the cross term, —i(f?/y)k-grad u,. 
This term will be small if 1 is fairly flat through- 
out most of the volume (free electron case). In 
any case, it will be small for small k. We may 
treat this term as a perturbation, and obtain the 
energy, E;, to the second order in k, by means of 
the ordinary Schrédinger perturbation theory.® 


4It might be better to take instead the potential 
V(r) + (3e?/2r,) — (e*r?/27,3). The last two terms represent 
the potential due to the space charge of the valence elec- 
trons in the polyhedron. In order to get the proper value 
for the interaction energy of the electrons, it would then 
be necessary to subtract approximately 1.2¢?/r, from the 
expression for W, given in Eq. (29). The final value for the 
total energy would probably be little changed. 

5 An alternative method for the calculation of the energy 
of electrons in excited states has been given by J. C. Slater, 
Phys. Rev. 45, 794 (1934). Slater applied the method to 
Na; it has been applied to Li by J. Millman, Phys. Rev. 47, 
286 (1935). The entire wave function, U;(x) exp [ck-x] is 
expanded in a series of spherical harmonics in each poly- 
hedron, and the boundary conditions are fitted at a number 
of points on the surface of the polyhedron. The disadvan- 








BARDEEN 


To this order, the energy may be written: 
E,=Eotah?k?/2y (5) 


since the linear term in k is zero. For perfectly 
free electrons, a=1. In general, a may be either 
greater or smaller than unity. 

It is our purpose to determine the value of a. 
We attempt to solve (4) by a series starting 
with uo: 


Uz, =Ugtuituet::s, (6) 
where — (h?/2p)Auo+ Vuo= Eons, (7) 
— (h?/2p)Auy+ Vuy= Eou,+i(h?/u)k- grad uo. (8) 


We consider only the first-order correction to the 
wave function, and thus determine the energy to 
the second order in k. The function 1; is odd, and, 
since it is periodic, it will change sign as we pass 
from one polyhedron into the next. The wave 
function will therefore be very small on the 
surface of the polyhedron. It is a good approxi- 
mation to again replace the polyhedron by a 
sphere, and to require that uw; vanish on the 
surface of the sphere. 

Wigner and Seitz find 2; by the usual pertur- 
bation theory, i.e., they expand 1, in terms of the 
characteristic functions of the unperturbed 
equation (7).° The solution may be obtained in a 
much simpler way, if we note that a particular 
solution of (8) is: 


uy= —i(k-x)uo(r), (9) 


as may be verified by direct substitution. To the 
particular solution, we must add a solution, 
i(k-x)v;(r), of the homogeneous equation in 
order to satisfy the boundary condition. Thus we 
take 

u;=1(k-x)v(r) =i(k-x)[01(r7) —uo(r) ] (10) 


and require that 
vi(r)=uo(r) at r=fz. (11) 


There is one solution of the homogeneous equa- 
tion which is finite at the origin; Eq. (11) 
determines the scale factor of this solution. Let 


tage of this method lies in the fact that it is necessary to 
expand the factor exp [sk-x]. This may introduce an error 
even though the potential is flat (cf. W. Shockley, Phys. 
Rev. 52, 866 (1937)). The method we have used is, of 
course, most satisfactory when the wave functions are 
nearly the same as those of free electrons, which is the 
usual case for the alkali metals. 
6 See especially F. Seitz, reference 2. 
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us set 
u(r) = R(r)/r, 
v(r) = P(r)/??. (12) 


The equations for R(r) and P(r) are then: 
@?R/dr?+(2u/h?)(Eo—V)R=0, (13) 
dP /dr?—(2/r?)P+(2u/h*?)(Eo—V)P=0. (14) 


The first is the radial equation for the s function, 
the second is the radial equation for the p 
function. These equations must be integrated (in 
general numerically) to find R and P. 

We now want to find the mean value of the 
energy using the wave function: 


ve =[uo(r) +1(k-x)o(r) ] exp [i(k-x)]. (15) 


We have 


E.= f vet —(H2/2u)A+VWidr. (16) 


eell 


The integration is over one polyhedral cell. The 
wave function is normalized so that 


Vi*Pidr=1. (17) 
cell 


Substituting (15) into (16), and making use of the 
Schrédinger equation and the normalization 
condition, we find that we may express E,, to the 
order k, in the form (5), with 


a=1+ f woL0(r) + (e-x)*(br)-*(do/dr) Wr. (18) 


After integration over the angular coordinates, 
we obtain: 


a=1t4e { uo(r)[o(%)+ (7/3) do/dr) War. (19) 


Integrating by parts, we may express a in the 
form: 


a=1 — (82/3) { v4r) (duo/dr)r-dr. (20) 


The integration may be carried out as follows.’ 
Let us put (cf. Eqs. (10) and (12)) 


v(r)=P/r?—R/r (21) 





‘ The possibility of carrying out the integration in Eq. 
(20) in terms of the boundary values of the wave function 
was pointed out to the author by C. W. Ufford. 
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AND Na 


and u(r) =R/r. 


With these substitutions: 


a=1—(8x/3) f "(P—»R)(dR/dr—R/*dr. (22) 
Now " 


[ rRGR/dr)ar= (rR?) =r, 
0 
— [(R—rRUAR/dr yar, 
0 
so that, to the approximation we are considering, 
(8x/3) f "pR(dR/de)dr = (4n/3)r,3Cuolr,) P—}. 
0 


Thus 
a=(4r/3)r,3[uo(r.) 2 


—(8r/3) {- P(AR/dr—R/rydr. (23) 
0 


This result may be expressed in the following 
form, after a further integration by parts: 


a= (47/3) | r3[ uo(rs) }?>—P(rs)Uo(1s) 


+ f “TR(dP/dr) —P(dR/dr)+2PR/r]dr}. (24) 
| 


From the differential equations which are 
satisfied by P and R, we have: 


R(@P/dr*) —P(@R/dr?) =2PR/r’, 
which may be transformed to 


d dP | dP dR 2PR 





—{ rR——rP—} = R—— P—+ 


. (25) 
dr dr dr dr dr r? 


The integration in Eq. (23) can therefore be 
carried out in terms of the boundary values of the 
wave functions. Making use of the boundary 
condition (11), we have, finally, 


a=y[(r/P)(dP/dr)—1Jrry (26) 


where y= (47/3)r,*[uo(r,) ? is the ratio of the 
square of the wave function at the boundary of 
the s sphere to its mean value. 
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Equation (26) provides a very simple method 
for the calculation of the Fermi energy. But a 
single integration of the differential equation 
(14) is required for each value of the lattice 
constant. 


NUMERICAL CALCULATION OF THE ENERGIES OF 
Li AND Na 


According to Wigner and Seitz,! the energy of 
a monovalent metal may be expressed in the 
following form: 


W=E.+ F+W.+W1. (27) 


This is the negative of the energy per atom 
required to break up the metal into positive ions 
and free electrons, and thus is the sum of the 
atomic heat of sublimation and the ionization 
potential. The first two terms represent the value 
of E; averaged over all the occupied states; Ep is 
the energy of an electron in its lowest state, and 
F is the additional energy required by the Pauli 
principle. The next term, W,, is a small correction 
to the interaction energy of the electrons (not 
taken adequately into account in the first two 
terms) which will be discussed at more length 
later on. The last term, W;, is the very small 
nonelectrostatic interaction energy of the ions 
(the overlap and van der Waals energies) which 
we will neglect in the present work. 

In order to obtain Eo, it is necessary to find the 
radial solution of Eq. (1) (or, equivalently, of 
Eq. (13)), subject to the boundary condition (2). 
The solutions, Yo(r, Eo) which are finite at the 
origin, depend on two variables, the radial 
distance r and the energy parameter, Ey. The 
actual value of Eo for a given value of r, is that 
for which (dyo/dr),-.,=0. It was desired to find 
E, for a fairly wide range of values of r, near the 
energy minimum. In this region, Eo does not vary 
rapidly, so that it was possible to obtain accurate 
values of ~o(r, Eo) covering the desired range by 
integrating Eq. (1) numerically, using but three 
different values of Ey for each metal. The values 
chosen for Li were —0.63, —0.65, and —0.67 Ry, 
and for Na —0.60, —0.62, and —0.64 Ry. The 
derivative, dy»/dr, for each of these solutions, 
was obtained by numerical methods. The energy, 
E,(r.), for a given value of r,, was then obtained 
by finding, by quadratic interpolation between 
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the three different energies, the value of Eo for 
which (dy/dr),-.,=0.8 The same ion-core fields 
used by Wigner and Seitz were also used in the 
present calculation.° The results are given in 
Tables I and II for Li and Na, respectively. 

Neglecting terms in k‘ and higher order, the 
Fermi energy is: 


F=3ahh?/10u=2.210/r2Ry, (28) 


where kp is the magnitude of & for the highest 
occupied state of the Fermi distribution. We have 
determined a from Eq. (26) after finding P(r, Ey) 
by numerical integration of Eq. (14). The ener- 
gies chosen were again —0.63, —0.65, and 


TABLES I AND II. 7, ts the radius of a sphere of atomic 
volume; Eo is the energy of an electron in its lowest state; y is 
the ratio of the square of the wave function at the boundary of 
the s sphere to its mean value; a is the effective number of free 
electrons per atom; F is the Fermi energy =2.21la/r.2; and W 
is the total energy per atom. Atomic units (reference 12) have 
been used. 








Energy of lithium 








's —Eo ¥ a F -W 
2.36 | 0.5727 1.179 0.2034 0.0807 0.4501 
2.52 .6219 1.165 .3282 .1142 4717 
2.68 .6509 1.151 4322 .1330 .4869 


2.84 .6646 1.134 5171 1417 4964 
3.00 .6682 1.114 5843 1435 5022 
3.16 6653 1.094 .6391 .1414 5048 
3.32 6582 1.072 .6839 1371 5049 
3.48 .6484 1.050 7211 1317 5031 














Energy of sodium 








's —Eo ¥ a F -W 
3.16 0.6619 1.1294 1.1392 0.2521 0.3907 
3.32 -6552 1.1033 1.1165 .2238 4152 
3.48 .6458 1.0779 1.1016 .2010 4312 
3.64 .6346 1.0531 1.089 .1816 4417 
3.80 .6225 1.0287 1.079 1651 4482 
3.96 .6099 1.0047 1.069 .1507 4520 
4.12 .5972 .9793 1.059 .1379 4537 
4.28 .5846 .9549 1.050 .1266 4540 
4.44 .5722 .9305 1.041 .1166 4530 























8 An alternative procedure would be to find the values of 
r, for which (dj/dr)r-r,=0 for each of the three energy 
values. This would yield six values of Eo(r,), there being 
two values of r, for each energy. One could then inter- 
polate between these values to obtain Ep for intermediate 
values of r,. This procedure is not as satisfactory as the one 
adopted, because, in the region of interest, interpolation 
between the different energies for a fixed value of 7; 'S 
more accurate than interpolation between different values 
of r, for fixed Eo. ne 

° The Prokofjew field was used for Na; the semi-empirical 
field of Seitz was used for Li. 
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—0.67 Ry for Liand —0.60, —0.62,and —0.64Ry 
for Na. From these solutions, values of rdP/Pdr 
were obtained by numerical methods. Interpo- 
lation between these gave the corresponding 
values for the actual energy, Eo, going with a 
particular value of 7,. To obtain a, values of 
y=(4x/3)r.*[Wo(r.) ? are also required; these 
were obtained by numerical integration. Values 
of y and ev are listed in Tables I and II for Li and 
Na, respectively. Values of the Fermi energy F, 
as obtained from Eq. (28), are also given in these 
tables. It seems to be very difficult to estimate 
the magnitude of terms in k‘ and higher order in 
the Fermi energy. They are probably much more 
important for Li than for Na. 

The rapid decrease of the values of a@ for Li 
with decreasing 7, in the neighborhood of the 
normal lattice constant has already been noted 
by Seitz? and by Millman.’ The decrease con- 
tinues throughout the range we have considered. 
The Fermi energy itself decreases below 7r,=3. 
Since the Fermi energy does not change rapidly, 
the normal lattice constant is determined mainly 
by the minimum of Ey. The Fermi energy of Na 
is but slightly larger than the free electron value. 

If there were never more than one electron in 
each polyhedral cell, the energy per atom would 
be Ey +F. To this we must add the correction 
term, W,, resulting from the correlations between 
the positions of the electrons.!° This correction 
term is really the sum of three terms: (1) 
Electrostatic interaction energy of a given 
electron with other electrons in the same poly- 
hedral cell, assuming no correlation between the 
positions of the electrons, (0.6e?/r,) ; (2) Exchange 
energy resulting from the ‘‘Fermi hole’’ between 
electrons of parallel spin (—0.458e?/r.); (3) 
Energy gain resulting from correlations between 
the positions of electrons of one spin with those of 
opposite spin (~ —0.58/(r.+5.1) Ry). The values 
given in parentheses are the free electron values 





Cf. Mott and Jones, reference 1, p. 137. 





ENERGIES OF METALLIC Li AND Na 371 














TABLE III. 
HEAT OF 
LATTICE SUBLIMATION ComPRESSI- 
CONSTANT (KCAL. /MOLE) BILITY 
(exp.) 3.46A 391 7.4 
Li 
(calc.) 3.49A 34 8.4 
(exp.) 4.25A 26! 12.3? 
Wa 
(calc.) 4.53A 23 12.0 

















1 Bichowsky and Rossini, Thermochemistry of Chemical Substances, 
(Reinhold, 1936). 
2 From Bridgman’s data (cf. following paper), 


which we will use in the present work." Thus, in 
atomic units,” 


W.=[0.284/r,—0.58/(r,+5.1)] Ry. (29) 


The total energy per atom (experimentally, the 
ionization potential plus the atomic heat of 
sublimation) is 


W=£ot+ F+ W.. 


Values of W for different values of r, are also 
given in Tables I and II. 

A comparison of the theoretical and experi- 
mental values of the lattice constants, heats of 
sublimation, and compressibilities of Li and Na 
is given in Table III. The theoretical values for 
Li are very nearly the same as those obtained by 
Seitz.2 The lattice constant of Na is smaller than 
that obtained by Wigner and Seitz (4.75A) and is 
in better agreement with experiment. The smaller 
lattice constant is due to the fact that our values 
for the Fermi energy do not increase as fast with 
decreasing 7, as do those used by the above 
authors. Our calculations permit fairly accurate 
determinations of the compressibilities, which 
previously had only been roughly estimated. The 
variation of compressibility with pressure is 
discussed in the following paper. 

11 Seitz (reference 2) has shown that, even for Li, the free 
electron values are close to the actual values. 

12 Unit of length is the Bohr radius, 0.528 X 10-8 cm; the 


Rydberg unit of energy is equal to the ionization potential 
of hydrogen (13.5 ev). 
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As a further test of the method of Wigner and Seitz 
(for the calculation of the energies of the monovalent 
metals), the pressure variation of the compressibilities of 
the alkali metals is determined and compared with recent 
experimental results of Bridgman which extend to pressures 
of the order of 45,000 kg/cm*. Fréhlich’s approximate 
analytic expression for the dependence of the energies of 
the alkali metals on volume suggests the following semi- 
empirical form: 


W=a/v+b/vi-—c/v. 
If the parameters are determined from (1) the lattice 


constant, (2) the energy, and (3) the compressibility (all 
at zero pressure), the expression gives values for the 


relative change in volume with pressure in close agreement 
with Bridgman’s results. The constants so determined 
agree fairly well with Fréhlich’s theoretical values. Com- 
pressibilities of Li and Na, as obtained from the more 
exact and fundamental calculations of the energies given 
in the preceding paper, are compared with the experi- 
mental values. Theoretical values for Na are in good 
agreement with experiment, but the corresponding values 
for Li are about fifteen percent too high. It is suggested 
that the polymorphic transition in Cs, which occurs at a 
pressure of about 22,000 kg/cm?, is a transition from the 
normal body-centered to a face-centered form, resulting 
from the nonelectrostatic interaction energy of the ions, 
which favors close packing. 





I. INTRODUCTION 


RIDGMAN! has recently obtained rough 

experimental values of the compressibilities 
of a large number of the more compressible 
substances in the pressure range extending to 
about 45,000 kg/cm. Of particular theoretical 
interest are the results he obtained for the alkali 
metals. Thanks largely to the work of Wigner 
and Seitz,? it is possible to make fairly accurate 
calculations of a large number of properties of 
this particular group of metals. A calculation of 
the compressibility and its variation with pres- 
sure provides a very sensitive test of the general 
methods. 

The instantaneous compressibility, K, of a 
metal at the absolute zero of temperature is 
determined from the energy, W, by means of 
the following equation: 


1/K = —v(dp/dv) =v(d?W/dv"), (1) 


where v is the volume of the metal and # is the 
pressure. Our problem is thus the calculation of 
the energy of a metal as a function of volume. 
We have given, in the preceding paper,® a 
calculation of the energies of Li and Na. In 
Section III of the present paper, we discuss the 


* Society of Fellows. 

1P. W. Bridgman, Proc. Am. Acad. 72, 207 (1938). 

2 E. Wigner and F. Seitz, Phys. Rev. 43, 804 (1933); 46, 
509 (1934). E. Wigner, Phys. Rev. 46, 1002 (1934); F. 
Seitz, Phys. Rev. 47, 400 (1935). 

‘d) the present issue. This paper will be referred to 
as (I). 


compressibilities and their variation with pres- 
sure. Fréhlich* has given an approximate ex- 
pression for the energy of an electron in its 
lowest state. Using the free electron value for 
the Fermi energy, he has determined the initial 
compressibilities of the alkali metals, and ob- 
tained rough agreement with experiment. A 
modification of Fréhlich’s procedure is discussed 
in Section II. At the end of the paper, we give a 
brief discussion of the polymorphic transition in 
Cs, discovered by Bridgman. 


II. Sem1-EMpPIRICAL CALCULATION OF THE VARI- 
ATION OF THE COMPRESSIBILITY 
WITH PRESSURE 


In the present section we assume that the 
energy of the metal (per atom) is the sum of 
the energy of an electron in its lowest state, Eo, 
and the Fermi energy, F. We neglect the small 
correction term to the interaction energy of the 
valence electrons (given by Eq. (29), reference 3) 
and also the nonelectrostatic interaction energy 
of the ions. 

According to Frohlich,‘ for values of r, neat 
the energy minimum (or, in the neighborhood 
of the normal lattice constant), the energy £0 
(in atomic units®) is given by the following 


4H. Froéhlich, Proc. Roy. Soc. 158A, 97 (1937); Electro- 
nentheorie der Metalle (Berlin, 1936) p. 272. 

5’ Unit of energy is the Rydburg (13.5 ev) ; unit of length 
is the Bohr radius (0.528 X 10-8 cm). 
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approximate expression, 
Ey=a/r,2—3/r,, (2) 


where 7, is the radius of a sphere of atomic 
volume (s sphere). The repulsive energy is 
inversely proportional to the atomic volume, and 
the attractive energy, resulting from the inter- 
action of the valence electrons with the positive 
ions is inversely proportional to the lattice 
constant. 

The assumptions on which this formula is 
based are not very clear from Frohlich’s treat- 
ment. It is not difficult to show that (2) follows if 


7 = (42/3)r.*[Yo(r.) P=1, (3) 


or, in other words, if the value of the wave 
function at the boundary of the s sphere is 
equal to its mean value. This relation must hold 
for all values of r, for which (2) is to be applied. 
One would expect that (3) will be approximately 
satisfied if the wave function of an electron in 
its lowest state is fairly flat throughout most of 
the volume. There is reason to suppose that this 
will be the case for values of 7, near the actual 
atomic radius, since both (dyo/dr),=r, and 
(d*)o/dr?),-r Vanish when £ is a minimum.® 

To derive (2) from (3), let us set R=ryo, so 
that the Schrédinger equation for R (in atomic 
units) is: 


d°R/dr?=(V(r) —Ey)R. (4) 


The potential V(r) is that of a free (positive) ion. 
The boundary condition of Wigner and Seitz? is: 


dR/dr=R/r at r=f"z. (5) 


Both Ey and R are functions of 7,, the point at 
which the boundary condition (5) is to be 
satisfied. Denoting derivatives with respect to 
r, by a prime, we have, differentiating (4), 


d?R' /dr?=(V(r) —Ey)R' —Ev'R, (6) 


so that, multiplying by R, integrating, and 
making use of (4) to eliminate (V—E o)RR’, we 
obtain : 


Ey i) R°dr = J ‘LR’ (@R/dr*) — R(@R' /dr?) Jdr 
=[R’(dR/dr)—R(dR'/dr)},-,. (7) 


‘For a discussion of this point, see N. F. Mott and H. 
Jones, Theory of the Properties of Metals and Alloys, 


Oxford, 1936), p. . 
lated in’(I), ), p. 79. Values of y for Li and Na are tabu 





Differentiating the boundary condition (5) with 
respect to 7;, we have: 


[(@R/dr*) + (4R'/dr)]—1,=LR'/r}=1, (8) 


Making use of this relation, we can eliminate R’ 
from (7) and obtain 


Ey! f “Ride =[R(@R/dr) Joon 
=[V(r.)—Eo(rs) JR(r.)?. (9) 


We thus get the following differential equation 
for Eo, if condition (3) is satisfied: 


dE/dr,=(3/r.)[ V(r.) —Eo(rs)]. (10) 


Froéhlich’s result (Eq. (2)) is obtained by in- 
tegrating (10) with V(r,)=-—e?/r,. The latter 
will hold for values of r, larger than the radius 
of the ion core.’ 

It is, of course, impossible to tell a@ priori 
whether or not Fréhlich’s formula will hold 
without making an explicit calculation of the 
wave functions using the actual potential of the 
ion, as Wigner and Seitz have done for Na and Li. 
Eq. (2) does give fairly good values for these 
metals. We will use, for comparison with the 
experimental data (of all the alkali metals) an 
expression of the form: 


Ey=a/r,2—c/rs, (11) 


treating @ and c¢ as empirical constants. The 
value of c so obtained is then compared with the 
theoretical value (3) given by Eq. (2). 
The Fermi energy F, expressed in atomic 
units, is 
F=2.21a/r2, (12) 


where a is the effective number of free electrons ~ 
per atom. If the electrons are free, a is equal to 
unity. In general, a may be either greater or 
smaller than unity, and may be a function of r,. 
It is known that the value of a for Li is much 
less than unity, and decreases as r, decreases, 
while for Na, a is slightly greater than unity and 
increases slowly as r, decreases.® As a first approx- 
imation, we may assume that a=1, so that the 
Fermi energy will be inversely proportional to 
the square of the lattice constant. 


7 The author is indebted to Professor E. Wigner for a 
discussion of Fréhlich’s formula. 
8 Tables I and II, reference (1). 
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Equations (11) and (12) suggest the following 
type of equation for the energy of an alkali 
metal : 


W=A (v/v) +B(v0/v)'—C(vo/v)?, (13) 


where A, B, and C are independent of the volume, 
and vp is the atomic volume at zero pressure. 
The first and third terms come from the energy 
of an electron in its lowest state; the second 
from the Fermi energy. To determine the values 
of the three parameters, the following experi- 
mental data are used (all refer to 0°K and p=0): 

(1) The lattice constant ; 

(2) The energy, given by the sum of the 
ionization potential and the atomic heat of 
sublimation, and 

(3) The compressibility, Ko. 

The values of the constants so obtained are 
checked in two ways. First, the empirical values 
of B and C are compared with the theoretical 
values given by Eqs. (12) and (2). The semi- 
empirical equation is then used to determine the 
change in compressibility with pressure (or, what 
amounts to the same thing, the relative change 
in volume with pressure) and the calculated 
values are compared with experiment. 

Since the energy is a minimum when v=vo 
(or p=0) we have: 


2C/3=2A+4B/3. (14) 
The energy at zero pressure is given by: 
—W=2A-+B. (15) 


Finally, the compressibility at zero pressure 
gives the relation: 


3/K,=2A+2B/3. (16) 


To obtain (15) and (16) we have used the value 
of C given by Eq. (14). 

Table I gives the relevant data for the calcu- 
lation of A, B, and C for the alkali metals. 
Row (d) gives the energies, obtained from the 
sum of the ionization potential and the atomic 
heat of sublimation. The compressibilities, given 
in row (e) were obtained by extrapolating the 
experimental results, taken at room temperature, 
to T=0°K by means of the Debye equation of 
state, as discussed below. From the values given 
in these two rows, the constants A and B were 
determined by use of Eqs. (15) and (16); the 
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results are given in rows (f) and (g). Row (i) 
gives the values of $C obtained by substituting 
these values of A and B into Eq. (14). Row (h) 
gives the theoretical values of B as obtained 
from Eq. (12) under the assumption that a=1 
(free electron value). Row (k) gives the theo- 
retical values of 3C obtained from Eq. (2) of 
Frohlich. 

The empirical values of B, depending on the 
small difference between rows (d) and (e), are 
not at all accurate, but they are rather small for 
all metals except Li. The anomalous value for Li 
results from the fact that @ is not constant, but 
decreases as the volume decreases. The Fermi 
energy cannot be expressed as a single term, but 
is given in fair approximation in the region of 
interest by the difference between a large positive 
term inversely proportional to the square of the 
lattice constant, and a negative term inversely 
proportional to the volume. The negative value 
of B obtained for K may perhaps be explained 
in a similar fashion. The differences between the 
theoretical and experimental values of B for Na, 
Rb and Cs, are probably within the rather large 
experimental uncertainties. 

A much more significant test of the theory is 
obtained from a comparison of the empirical and 
theoretical values of C given in rows (i) and (h). 
The agreement is good for all metals except K, 
for which the empirical value is somewhat 
smaller than the theoretical. 

As a second (experimental) test of the semi- 
empirical equation (13), we compare the relative 


TaBLeE I. Calculation of the parameters in the semi- 
empirical equation for the energies of the alkali metals: 
W=A(v0/v) + B(v0/v)?/8— C(vo/v)"8. The theoretical values 
of B are obtained from the free electron formula for the Ferm 
energy; theoretical values of C are obtained from Frohlich’s 
formula for the energy of an electron in its lowest state. Units: 
row a, 10-** cm’; rows b and c, ev/atom; rows d to j, 10" 
erg/atom. 
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1 Bichowsky and Rossini, Thermochemistry of Chemical Substances, 
(Reinhold, 1936). 
2 From Bridgman’s data, extrapolated to T =0°K. 
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change in volume with pressure given by the 
theory, with the experimental values of Bridg- 
man.! It is first necessary to extrapolate the 
observed values to T=0°K. To do this, we use 
the Debye equation of state,’ although it 
probably does not apply accurately to the alkali 
metals. The Debye equation is: 


p= —(dW/dv) —3kT(d0/0dv)D(0/T). (17) 


Here W and v are the energy per atom and the 
atomic volume, respectively, and @ is the charac- 
teristic temperature. The function D(x) is de- 


fined by: — 
D(x) =— 
x3Jq e&—1 





(18) 


For small x (T large), D(x) <1. The value of 
d6/dv is obtained from Griineisen’s law: 


3a/K = — (d6/dv)(cv/6), (19) 


in which a@ is the linear coefficient of thermal 
expansion,* K is the compressibility, and c, is 
the atomic heat capacity : 


Cy=3kLD(0/T)—(0/T)D'(6/T)]. (20) 


Using (19) and (20), we have, for small values 
of 6/T, 
p= —(dW/dv)+Ap, (21) 


where 


Ap=(3aT/K){1—(30/8T) 
+(1/10)(0/T)?—--+-]. (22) 


The second term of Eq. (21), Ap, represents 
the pressure due to the thermal motion of the 
ions, which must be subtracted from the ob- 
served pressures (at room temperature) to obtain 
the pressure at T=0°K. Although this term 
depends on the volume, the change is probably 
not very great except at very high pressures. 
In this region, the correction term is but a small 
fraction of the total pressure. We have therefore 
assumed that Ap is independent of volume, and 
to determine its value we have used the values 
of a, K, and @ for zero pressure and room tem- 
perature. Table II gives the values of Ap for the 
alkalis, together with the relevant data used in 





_ *See, for example, J. K. Roberts, Heat and Thermodynam- 
as (Blackie, 1933). 

* Not to be confused with the effective number of free 
— per atom, for which the same symbol has been 
used, 


the calcilation. The total relative change in 
volume due to thermal expansion from room 
temperature down to T=0°K is also given. 

A plot of the relative change in volume of the 
different alkalis as a function of pressure is 
given in Fig. 1. The heavier lines represent the 
experimental values extrapolated to the absolute 
zero of temperature; the lighter lines are the 
theoretical values obtained from the semi- 
empirical formula (13). From this equation we 
find 


pro =y*(y—1)[2A+9B+A(y—1)], (21) 
with y = (v0/v)3. 


The values of A and B which were used are 
given in rows (f) and (g) of Table I. 

The agreement is good for all metals except 
Rb. Professor Bridgman informs me that the 
initial compressibility for this metal is less certain 
than for the other alkalis. The initial com- 
pressibilities were used in the determination of 
the parameters, so that the slopes of the theo- 
retical curves agree with the slopes of the 
experimental curves at the origin. It is the shape 
of the curve which is important for a test of the 
validity of the theory. The shape is determined 
mainly by the factor y(y—1) in Eq. (21). 
(The term A(y—1) in the brackets is small.) 

There is a tendency for the theoretical curves 
to be slightly above the observed at the highest 
pressures. The nonelectrostatic interaction en- 
ergy of the ions, which we have neglected, is not 
large enough to account for this discrepancy. 
The general agreement would have been better 
if the parameters had been adjusted: to make 
the experimental and theoretical curves match at 
some high pressure, rather than to fit the initial 
compressibilities. 


III. COMPRESSIBILITIES OF Li AND Na 


It is much more desirable to have a calculation 
which is based on fundamental principles, so that 
there are no parameters to be determined from 
experiment. In order to determine the com- 
pressibility and its variation with pressure, it is 
necessary to have an accurate calculation of the 
energy as a function of volume. Such a calcula- 
tion for Li and Na is given in the preceding 
paper.’ In Fig. 2, we show the theoretical curves 
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of the compressibilities as a function of volume, 
as determined from these energy values. For 
comparison, experimental curves derived from 
Bridgman’s data,! extrapolated to T=0°K, as 
indicated above, are also given. The abscissa 
gives the relative change in volume from the 
observed volume at zero pressure. The smallest 
volumes correspond to a pressure of about 
40,000 kg/cm?. The curves for Na are in very 
good agreement, but the theoretical values for Li 
are about 15 percent too high. The difference in 
the latter case is perhaps due to our neglect of 
terms of higher order than those in k? in the 
calculation of the Fermi energy. 

The reciprocal of the compressibility is given 
essentially by the second derivative of the 
energy with respect to volume. It is interesting 
to note that the second derivative of the Fermi 
energy of Li is negative, and that the main 
contribution to 1/K comes from the energy of 
an electron in its lowest state, EZ». On the other 
hand, the second derivative of Ey is small for 
Na, and the main contribution comes from the 
rapid increase of the Fermi energy with de- 
creasing volume. 


IV. POLYMORPHIC TRANSITION IN Cs 


At ordinary pressures, the alkali metals are 
body centered. It is not yet known why this 
structure is preferred over, for example, the 
close packed cubic structure (face centered). 
Fuchs” has shown that, for equal atomic volume, 
the energies of the valence electrons in the two 
structures differ by only about 10-* ev per 
electron. He attributes the face-centered struc- 
ture of the noble metals (Cu, Ag and Au) to the 
repulsion between the closed d shells of the ions, 
which naturally favors the close packing. The 
nonelectrostatic energy of the ions in the alkali 


TABLE II. Data used in the extrapolation of Bridgman’s 
results to T=0°K. 











METAL Li Na K Rb Cs 
a X 105 5.6 ton 8.3 9.0 9.7 
(1/K) X10-® (dynes/cm?) | 120 58 32 23 19 
6/r 1.5 0.6 0.40 0.29 0.20 
p X10-* (dynes/cm?) 3.5 2.8 2.0 1.6 1.5 
a (calc.) 0.031] 0.043} 0.060} 0.080] 0.084 























10K, Fuchs, Proc. Roy. Soc. A151, 585 (1935). 
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Fic. 1. Relative change in volume of alkali metals with 
pressure. Light lines are theoretical values obtained from 
the semi-empirical equation (13). Bridgman’s experimental 
data, extrapolated to T=0°K, are represented by the 
heavy lines. The initial compressibilities have been used to 
fix the values of the parameters in the semiempirical equa- 
tion, so that the slopes of the experimental and theoretical 
curves are in agreement at the origin. The break in the 
curve for Cs is due to a polymorphic transition. 


metals (resulting both from the exchange re- 
pulsion and the van der Waals energy) will also 
be lower in a face-centered lattice than in a 
body-centered lattice, but this energy is so small 
at ordinary pressures that it has little influence 
on the structure. As the volume is decreased, 
the energy of the ions will increase relatively to 
the energy of the valence electrons, so that one 
would expect that at very high pressures the 
alkali metals would take up a close packed 
structure. The rare gas solids are, of course, 
face centered, and it is likely that the alkali 
metals are also face centered at very high 
pressures. 

The metal for which the transition should take 
place at lowest pressure is Cs, since it not only 
has the largest ion of the alkali metals, but also 
has the highest compressibility. Bridgman has 
found a polymorphic transition in Cs at a pres- 
sure of about 22,000 kg/cm*®. The change in 
volume at the transition point is 0.00336 cm*/g, 
and the latent heat, determined from the change 
in transition pressure with temperature is 8.49 
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Fic. 2. Compressibilities of Li and Na as a function of 
pressure. Theoretical values are obtained from direct cal- 
culations of the energies (given in the preceding paper); 
— from Bridgman’s data, extrapolated to 


kg cm/g (0.2 cal./g). The small latent heat 
shows that the difference in entropy between the 
two phases is small. 

A rough test of the theory that the transition 
is due to the nonelectrostatic energy of the ions 
may be obtained from the change in volume at 
the transition point. Let the free energy of the 
low pressure phase be F,(v) that of the high 
pressure phase be F2(v), and the corresponding 
volumes at the transition point be 2; and 2. 
Expressing the condition that the pressures of 
the two phases are equal at the transition point, 
we have, for v;—v2 small, 


d(F,—F:) V1 —Ve 1 
———| sie (—) , (22) 
dv v=v} V1 K, v= 


where K, is the compressibility of the low 
pressure phase. The right-hand side may be 
evaluated from experimental data. The result is 
0.83 10° dynes/cm?." In order to obtain a 








4 This value is for room temperature. In order to com- 
pare with theory, we should take the value at absolute zero. 
Without a knowledge of the thermal expansion in the two 





rough theoretical value for the left hand side of 
the equation, we will assume that the main 
contribution comes from the difference in the 
repulsive energy of the ions in the two phases. 
The difference in energy of the valence electrons, 
which presumably favors the body-centered lat- 
tice, will not change rapidly with volume, so 
that the derivative will be small. If we assume 
that the van der Waals energy of the ions varies 
as 1/r®, the difference between its values for 
the body-centered and face-centered structures 
is very small. 

The exchange repulsion between the alkali 
ions has been estimated by Mayer and collabo- 
rators® from an analysis of the alkali-halide 
crystals. They find, for the mutual repulsive 
energy of two ions a distance 7 apart, 


W=A; exp [(2r:—1)/p], (23) 


where, if A; is taken to be 1.25X10-" ergs, the 
ionic radius 7;=1.455A for the Cs ion, and 
p=0.345A for all the alkali and halide ions. 
Using these values, we have estimated the con- 
tribution of the repulsion of the ion cores to (22) 
and find a value which is only about one-fifth of 
that observed. However, direct calculations" of 
the repulsion between the rare gas atoms He 
and Ne show that the energy can be given 
approximately by the form (23), with p=0.21A, 
a value much smaller than that of the above 
authors. The larger value probably results from 
the fact that it is a kind of a mean between the 
values for the alkali ions and the halide ions, 
and the latter are probably much larger than 
the former. If we take p=0.21A for thé Cs ion, 
and take the same value for A ;(1.25 X 10-” ergs) 
we estimate, from the lattice constant of the 
corresponding rare gas solid, that r;~1.75 
—1.80A. With these values, one obtains results 
of about the right order of magnitude for (22). 
Due to the difficulties in estimating the differ- 
ence between the energies of the valence electrons 
in the two different structures, and also to the 
difficulty in estimating the interaction energy of 


phases, the change with temperature cannot be deter- 
mined, but it is probably not very great. As we are inter- 
ested only in the order of magnitude, we use the room- 
temperature value. 

12 M. Born and J. E. Mayer, Zeits. f. Physik 75, 1 (1932); 
M. Huggins and J. E. Mayer, J. Chem. Phys. 1, 643 (1933). 

13 J. C. Slater, Phys. Rev. 32, 349 (1928) (He); W. E. 
Bleik and J. E. Mayer, J. Chem. Phys. 2, 252 (1934) (Ne). 
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the ions, one cannot be certain, but it does seem 
probable that the transition in Cs is from a 
body centered to a more close packed structure, 
which is probably face centered. 
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The author wishes to express his gratitude to 
Professors J. H. Van Vleck and E. Wigner, who 
have kindly read the manuscripts of both the 
present paper and the preceding one. 
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Electronic Structure of the First Row Hydrides 


G. W. Kinc* 
Princeton University, Princeton, New Jersey 
(Received April 27, 1938) 


Stehn’s formulae for the energies of the first row hydrides, BH to FH, in terms of Heitler- 
London exchange integrals have been applied to more recent data. The empirically derived 
values for the integrals are now smooth functions of the atomic number, Z. They are also 
considered as functions of the internuclear distance, 7., and taken in conjunction with data 
on BeH, and with the facts that r, and the force-constant, &, are smooth functions of Z, but 
invariant to electronic excitation, the results show that it is 2s?—1s exchange that predomi- 
nates in fixing the internuclear distance and force-constant. On the contrary, the dissociation 
energies cannot be expected to be smoothly dependent on Z, r, or k. The heat of dissociation 
of NH is calculated by Stehn’s formulae, after interpolating one exchange integral. 


INTRODUCTION 


RESENT theories of heats of dissociation of 
molecules with more than four electrons are 
quite crude, and any correlation between theory 
and experiment, even for the simplest molecules 
would have value in giving theoretical thermo- 
chemistry a better quantitative basis. Stehn ' 
made a considerable contribution in deriving 
formulae for the energies of the first row hydrides, 
BH, CH, NH, OH, FH, in terms of exchange and 
Coulomb integrals which can be determined em- 
pirically and tested for self-consistency. Since 
his publication new data have appeared, and in 
this paper the integrals are recalculated. The 
three classes of integrals are now smooth func- 
tions of the atomic number, and combined with 
facts known about the internuclear distances and 
force-constants give a more detailed picture of 
the orbitals in these molecules. 


REDETERMINATION OF THE INTEGRALS 


Stehn applied the usual Heitler-London theory, 
building the molecules from a hydrogen 1s elec- 
tron and various configurations of the heavy 
atom, to obtain expressions for the electronic 


* National Research Fellow in chemistry. 
1J. R. Stehn, J. Chem. Phys. 5, 186 (1937). 


energy 
W=A+a,B+a2K +03J.+asJoq, (1) 


where A is the atomic level of the heavy atom 
into which the molecule adiabatically decom- 
poses, B is an atomic exchange integral, K the 
exchange of the hydrogen 1s electron with s 
electrons plus Coulomb terms, J, and J, the 
exchange with 7 or op electrons. The a’s are 
numerical coefficients of the matrix elements. 
Stehn noticed that the internuclear distance re- 
mains practically the same in all excited states of 
any one molecule, so that the J’s and K remain 
the same, permitting simultaneous equations in 
them to be set up and solved. In this application 
then the principal defect of the Heitler-London 
method, that of nonorthogonality, cancels out; 
but other corrections, such as higher order terms 
destroy the linearity of Eq. (1). 


DISCUSSION OF THE DATA 


The excitation energies Ymo1 and Vatom are 
usually known quite accurately, and sometimes 
the J integrals can be calculated from them alone. 
But a complete solution also depends on one dis- 
sociation energy, which is invariably the weakest 
datum. However, for any one molecule only one 
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ELECTRONIC STRUCTURE OF HYDRIDES 


such energy need be known, by virtue of the 
well-known cycle connecting upper and lower 
level dissociation energies, D,’ and D,”’, 


Ymoit D,.’ = Vatom +D,”. (2) 


The literature on the dissociation energies of 
this series of hydrides has been investigated 
thoroughly. In our calculations we have not used 
any approximate D’s, e.g., from the formula 
w?/4x.e, or from extrapolation to the limiting 
vibrational quantum numbers, but only those 
obtained from dissociation through rotation, 
which gives an upper and a lower limit to D.,, 
(BH, CH), or from thermal data, (OH, FH). 


Further details 


An error in D for CH or NH will not change the values 
of their J, or Jz, which are obtained from the v’s. The 
CH datum is reliable, so K of CH is a good point, Fig. 1. 
There is not enough data for the evaluation of all three 
integrals in the case of BH. However, if J, of BH is 
taken to fit in with J, of C, N, OH, then J, automatically 
will fit in with J, of C, N, OH, and K will fit with the 
K’s of CH and OH. Since J, and J; of NH fall smoothly 
on their respective curves, we can expect to interpolate 
K of NH. 

There is no value for the heat of dissociation of NH; 
but in the appendix we show this interpolated value of K 
is probably very good. 

The heat of dissociation of FH can only be obtained 
with the help of thermochemical data. Three investigators 
report? a heat of formation of +63.0 kcal. To get the 
dissociation energy of FH, the dissociation energy of F2 
is also needed. This has been evaluated as —62.6 kcal. 
by Franck’s method from the onset of continuous absorp- 
tion in the spectrum.*? This value can be verified from 
spectroscopic constants‘ as follows: 

The upper state of F2, ‘Il, has an w.=837 cm™ and 
anharmonicity xwe=140. The formula w2/4x, gives a 
value of D, that is usually too high, but the absolute 
error in this case is very small, of the order of magnitude 
of D,’ itself, 0.155 ev. Mulliken® has shown that the 
upper state dissociates into normal atoms. Hence the 
dissociation energy of the lower state, D,’’,(which cannot 
be reliably determined by w,.’?/4x,"w,'"), is D.’+vmol 
=0.155+2.575=2.73 ev or 63 kcal, as above. The D, 
for HF then is 6.35 ev. 

Kirkpatrick and Salant® obtained bands with high 

*Bichowski and Rossini, Thermochemistry of Chemical 
Substances (1937). 

*von Wartenburg, Sprenger and Taylor, Zeits. f. physik. 
Chemie, Bodenstein Festschrift, 61 (1931). 

*Sponer, Molekulspektren I, p. 18. 


*R. S. Mulliken, Phys. Rev. 36, 701 (1930) 


1938) E. Kirkpatrick and E. O. Salant, Phys. Rev. 48, 945 
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Fic. 1. Values in electron volts of the exchange integrals 
J, and J,, the total Coulomb terms plus 2s?— 1s exchange 


K, and the 2—1s Coulomb term Kz, as functions of the 
atomic number. 


vibrational quantum numbers in the infra-red, and extra- 
polated to get D,’’=6.1 ev.” 


DISCUSSION OF THE INTEGRALS 


The solutions obtained for J,, J, and K are 
plotted against the atomic number Z in Fig. 1, 
and against the internuclear distance in Fig. 2. 
The principal remark is that the integrals are 
fairly constant for the series and that the curves 
are quite smooth. The values for K especially 
have become less erratic than in Stehn’s figure, 
and show less variation. The fact that fairly con- 
stant values are obtained, even though the inter- 
nuclear distance has changed markedly (from a 
constriction of the wave functions with increas- 
ing Z) shows that overlapping occurs to nearly 
the same extent and in the same relative 
positions. 

The integral J, is large, as to be expected. As 
will be shown below the equilibrium distance is 
fixed by 2s?—1s interchange. The minimum of 
J, is beyond this point, because the radius, and 
also the maximum value, of a 2p electron in a 
“dumbbell” orbital is larger than that of the 2s 
electron. The relative positions of the minima 
of 1s—2s and is—2p exchange can also be 

7 Landolt-Bérnstein, Tabellen, Vol. 3, p. 2873, state 
thermochemical data gives D,’’=120 kcal. or 5.2 ev, but 
this probably was old thermal data, ibid. Vol. 2, p. 1489, un- 


corrected for polymerization, ibid., Vol. 3, p. 2746 and also 
Bichowski and Rossini, reference 2. 





G. W. 











? 
085 09 








0.5 


Fic. 2. Values of the integrals J;, J, and K as functions 
of the internuclear distance. The dotted curve represents 
Jz as calculated from Rosen’s expression for this integral. 


estimated from the integrals derived by Rosen.*® 
The J, integral is quite flat, so that the points on 


Figs. 1 and 2 are very close to the maxima of 
this integral. 


The second-order effect, the slight increase of J, with 
increasing 7, can also be accounted for by simple exchange 
integral theory. Rosen evaluated the integral between a 
2p and a 1s electron, which is, except for 1/ri2 terms, 


Je= NZ°7.(Zire+1)e~71"e, (3) 


N being unity minus a nonorthogonality function, fairly 
independent of Z;. In this derivation both nuclei are 
supposed to have the same charge Z;. Since the heavy 
nucleus is heavily shielded in these molecules, (3) should 
not be a bad formula for comparisons. It can in fact be 
fitted to the straight line of J, versus r. of Fig. 2 by two 
parameters, the ‘‘constant” N and a shielding constant 
of 0.63 for each p electron on the # electron that is ex- 
changing. Usual values of shielding of 1s?2s* were taken,° 
but the calculations are insensitive to them. Shielding 
constants are discussed again later. The same value is 
obtained assuming Z, to be the average of the effective 
charge of the heavy nucleus and proton, indicating (3) 
may be fairly good even when there is asymmetry. 


The integral J, decreases with increasing 7”, 
and quite suddenly farther out, as would be 
expected of overlapping of the hydrogen 1s with 

8 N. Rosen, Phys. Rev. 38, 2009 (1931). 


®Slater and Frank, Introduction to Theoretical Physics, 
p. 432. 
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m electrons. As r, decreases to form the united 
atom, J, becomes an atomic exchange integral 
between a 2p, and a 2p, electron, equal to 3F, 
(Shortley-Condon notation). An extrapolated 
value of 3F2; through C, N, O, F, is just about the 
value of 0.94 ev to which J, of Fig. 2 converges. 

The quantity K is actually a sum of several 
terms; the exchange and Coulomb integrals of 
2s*—1s binding, say K,, and Coulomb integrals 
of the hydrogen 1s with the p, and , electrons, 
K, and K, in Stehn’s notation. To solve the 
equations Stehn made the reasonable assumption 
that K,=K,=K,. Thus 


K=K,+0;K,+a¢K, 


aK, 40%. (4) 


where Z, is the number of # electrons. 

In view of the constancy of J,, the 2s?—1s 
exchange K, should also remain reasonably con- 
stant (as contrasted with the variation in K), 
and should be equal to the dissociation energy of 
BeH, which has no p electrons. The difference 
K,=(K-—K,)/Z, is calculated and plotted on 
Fig. 1. This, then, is the Coulomb energy (other 
than the attraction —e?/r, of each electron to the 
other nucleus, which is canceled by the nuclear 
repulsion term in the energy expression) of the 
hydrogen 1s and a 2 electron. As the series ap- 
proaches more to a united atom K, falls off toa 
small value of 0.2 ev, of the right order of magni- 
tude for an atomic electrostatic repulsion of two 
2p electrons. 
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Fic. 3. The linear dependence of 1/(re—0.34) on the 
atomic number Z, and the effective nuclear charge 4: 


For most excited states the internuclear distance is the 
same as in the ground state. 
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Fic. 4. Force-constant as a function of the atomic number, showing the invariance to 
electronic excitation. 


BOND STRUCTURE 


These results can be supplemented with the 
following facts (see Figs. 3 and 4, and Table I): 


(1) the internuclear distance, 7., and 

(2) the force-constant, k, are smooth, but rapidly 
varying functions of Z (including BeH which 
has no p electrons), but 

(3) and (4) are practically the same in the excited 
states of any one molecule. 


On the other hand the integrals J, and J, rep- 
resenting exchange of the hydrogen 1s electron 
with p, or p, electrons of the heavy atom are 
much less dependent on Z, and of course do not 
always appear in the energy expressions (for the 
coefficients a3 and a, in Eq. (1) are far from being 
smooth functions of Z). 

Apparently, then, the exchange integrals with 
the p electrons do not determine 7, or k to any 
significant extent, which is especially surprising 
in the case of the o electrons inasmuch as the 
7 integral contributes most to the energy." It is 
to be concluded that the equilibrium distance 
uieeeenionn 

” Since K is the sum of the Coulomb forces (positive) and 


exchange with 2s?, the latter will be somewhat larger than 
, In fact equal to K,, and be equal or larger than Jz. 


and force-constant are predominantly deter- 


mined by exchange with 2s electrons, and that 
the exchange energy with #,, although large, has 
a flat minimum, as is indeed the case as shown by 
evaluation of Rosen’s expression for this integral. 


Relative effect of different components of a potential curve 


An integral that gives a large contribution to the 
stabilization energy does not necessarily affect r, or & 
very much, for it may be flat. This phenomenon was 
emphasized by Pauling" in the case of the carbon-carbon 
bond distance. Even a small introduction of double or 
triple bond character shortens single bonds, whereas a 
similar amount of single bond wave function does not 
affect (i.e., lengthen) double or triple bonds. For consider 
a molecular potential curve that is made up of a linear 
combination, with coefficients x; of component molecular 
orbital potential curves of different curvatures (i.e., 
force-constants) &;, each with its own minimum at 7, 


V= Lixsk; (r%; - r)?, 


The equilibrium distance is the weighted mean 


(5) 


(6) 


with weight factors &;. That is, the most curved compo- 
nents contribute the most. A potential curve resulting 
from the usual type of exchange integral will have a 
curvature k; roughly proportional to 1/r;2. Hence a curve 


.= Lixikir;/ Lixirs 


1! Linus Pauling, Princeton lectures, 1938. 
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second differential of (1), viz., 
k=a2K"+a;J," 


Spectroscopic data show 2& is in 


so the quantities J,’ and J,’’ must be small relative to 
K"’; that is, the energy curves for exchange with # electrons 


are relatively flat. 


RELATIONS BETWEEN INTERNUCLEAR DISTANCE, 
FoRCE-CONSTANT AND NUCLEAR CHARGE 

With this picture of strong 1s — 2s? bonding it is 

possible to derive laws relating r., k and Z by a 


with a minimum relatively far out, as our J’s, will not 
affect appreciably the 7, of the whole molecule. 
The force-constant of the molecule is given by the 


G. W. KI 








+agJq”. (7) 


dependent of a3 and ay, 





TABLE I, 


NG 


simple theory of the dependence of the dimen- 
sions of the 2s? orbits on Z. These relations are 
smooth functions, which fit the data (including 
the excited states) very well for the series 
F, O, N, C, B, BeH. The values for LiH show a 
discontinuity, as is to be expected at the transi- 
tion from the configuration 2s? to 2s. But there 
is no discontinuity apparent on the change of 
number of # electrons, either 7 or o. These elec- 
trons have an effect only through shielding. 

If then, the potential curve has a component 
in the form of 2s?—1s exchange which predom- 
inates because of a sharp minimum at the short- 
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est r,, the whole series Be to FH should have 7, 
and k of the same potential curve, variation being 
due solely to contraction with increasing Z. 


Internuclear distance 


It is well established from measurements of the 
interatomic distances and from the simple idea of 
maximum overlapping that the internuclear dis- 
tance is roughly the sum of the atomic radii. For 
2s*—1s binding 


re=(n?/Z—S) +a. (8) 


Here S is the shielding constant for the heavy 
atom, and 7? is a little more than 2? in the case 
of s electron binding.'? The distance a can be 
taken either as the hydrogen radius, of 1 atomic 
unit, or as the position of the maximum of the 1s 
hydrogen wave function from the hydrogen 
nucleus, 0.34 atomic units. 

Curve I of Fig. 3 is a plot of 1/(r,—0.34) versus 
Z, showing linearity in agreement with (8). 

Since the 2s interchange fixes r., S refers to the 
shielding of a 2s electron. It is equal to A—BZ,, 


2 Reference 9, p. 445. 


where A is the 1s?2s shielding, B the shielding 
constant for a 2p on a 2s electron, and Z, is the 
number of 2 electrons. Thus 
1 zZ* 
=——, (9) 
re—a on 


Z*=Z(1—BZ,)—A. (10) 


For A the usual? 0.85 for each 1s electron, 0.35 
for the other 2s electron were used. The constant 
B then plays the role of a parameter which is 
adjusted to make the graph of 1/(r,—0.34) versus 
Z* pass through the origin in accordance with 
(9). It turns out that the points still lie on a 
straight line, Fig. 3, with a slope of 4.28, which is 
a good value for n?. 

The required value of 0.77 for the shielding of 
p electrons (along the internuclear axis) is larger 
than the values found in the central fields of 
atoms. This is to be expected as p, electrons in 
cylindrical fields of molecules are well inter- 
mingled with the bonding electrons, and even p, 
electrons have orbits located along the inter- 
nuclear axis, not “‘dumbbells’’ at the end of the 


TABLE II. 
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molecule. There is some indication that 7 elec- 
trons have a smaller shielding effect than the 
o type. 

This value of 0.77 for 2s—2p shielding is in 
agreement with the 2p—2p shielding of 0.63 ob- 
tained from our discussion of the exchange inte- 
gral J,, and further, is the only one that fits the 
following equations for the force-constants. 


Force-constants 


From the theoretical calculations of the poten- 
tial energy curves for Hz, Het, HeH*+ and other 
studies of interatomic forces, it is seen that the 
attractive forces fall off as — exp (—ar), or as 
—a’'/r®, while repulsive forces have a range about 
the square of this, exp (— 2a7) or b/r'®. Following 
Slater and Frank” we can estimate the constant 
a as Z*/r., where Z* is the effective nuclear 
charge for the principal binding electrons (10), 
giving a potential curve of the Morse type, 


De-22*r—r)Ir¢ ae 2De—2*(r—- 0)! 6, (1 1) 
The force-constant for this function is 
k=(2Z2*?/r2)D. (12) 


If it is true that 2s?—1s exchange fixes k in 
the series BeH to FH, the D in (11) and (12) will 
be practically constant, and equal to K,; there- 
fore the quantity kr?2/2Z*? should be roughly 
constant, equal to K,. Column (11) in Table II 
shows that the data are in fair agreement with this 
simple theory. The trend in K, is almost the 
same as for J,; the average value is 2.94 ev, 
whereas the dissociation energy of BeH, also 
approximately equal to K, is 2.22 ev. 





Heat of dissociation of NH 


There is no known value of the dissociation energy of 
NH. [The value of 102.2 kcal. quoted by Bichowski and 
Rossini? from Bates" is an upper limit. ] 

If we can rely on the smooth curves of Fig. 1 we can 
calculate, by means of Stehn’s formulae, this dissociation 
energy as D.=3.61 ev, or 83 cal. The thermal dissociation 
energy, Do will be D.—3w.=79 kcal. 

13R. M. Badger, J. Chem. Phys. 2, 128 (1934). 


4 ©=Bates, Zeits. f. physik. Chemie, Bodenstein Fest- 
schrift, 329 (1931). 


W. 
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Equations (9) and (12) can be combined to express the 
force-constant & as a function of Z* alone, 


2K .Z*4 





k= . (13) 
(n?)?+-0.68(n2)Z*+0.122" 


The best fit, column (9) of Table II, is with values of 
K,=2.94 ev and (n?) =4.29, values in excellent agreement 
with those obtained with the simpler formulae. At first 
sight this might be expected, but had there been a real 
trend from (9) and (12), the constants from (13) would 
not have been the same. 

As a comparison, Badger’s constant! is calculated in 
column (12), to show the agreement with the theoretical 
is as good as with the empirical formulae. 


SUMMARY 


By the aid of Stehn’s energy expressions we 
have made some analysis of the potential curves 
of the molecules BeH, BH, CH, NH, OH, FH. 
The force-constant and the internuclear distance 
are determined by 2s?—1s exchange, and are not 
influenced by the # electrons. This accounts for 
their regular variation with Z, and invariance to 
electronic excitation. On the other hand this 
regularity is entirely dispersed in the case of the 
dissociation energies. That is, we cannot expect 
graphs of D versus Z, r. or k to be simple curves, 
because in the formulae for D the p electrons con- 
tribute largely, in a fashion quite irregular with 
respect to Z, r, or k, as displayed explicitly by 
Stehn’s formulae. 

The author takes pleasure in acknowledging 
the help of Professor J. H. Van Vleck and of Dr. 
J. R. Stehn in preparing this paper, and also the 
hospitality of the physics, chemistry and mathe- 
matics departments of Princeton University. 


TABLE III. Bond energies of simple and saturated hydrides. 








CH 78.9 kcal. 4 CH, 87.6 kcal. 
NH 79 kcal. 4 NH; _ 83.9 kcal. 
OH = 100.2 kcal. 4 OH, 110.4 kcal. 





— 





This seems to be a good value in comparison with } of 
the heat of formation of NH; from atoms, namely 84 kcal., 
a calculation suggested by the theory of directed valency. 
Table III is a summary of such a computation for CH. 
NH, and OH. 
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Dielectric Absorption in Polar Media and the Local Field 
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The phenomena of dielectric dispersion and absorption in polar liquids subject to applied 
high frequency electrical fields are discussed on the basis of a modified Onsager theory for the 
local field. The results of this theory agree quite as well with experiment as do those from 
Debye’s application of the hypothesis of hindered rotation with which the present work is 
compared. The application of the modified theory to the static case is also considered. A law 
of corresponding states is found to exist for the alcohols and an empirical equation proposed 
by Van Arkel and Snoek for ‘‘normal” polar liquids is shown to have an explanation in terms 
of the local field. The desirability of more complete measurements of dielectric constants and 
absorption as a function of temperature and frequency is stressed. 





I. INTRODUCTION 


HE familiar equations of Debye for the real 
and imaginary parts (e’, e’’) of the dielectric 
constant of a polar liquid in an alternating elec- 
tric field E= Ey exp (iwt) of frequency v (v=w/27) 

are 
e' = €a+(€0—€0)/(1+2?), 


fele~edalled®. (1) 


Here €9 is the dielectric constant at zero fre- 
quency, due to the permanent moment plus in- 
duced polarization, and «., is the dielectric con- 
stant at infinite frequency due to the induced 
polarization only. The frequency dependent 
variable x is given by 


x= (€9+2)wr/(€o.+2), (2) 


where 7 is the so-called relaxation time. 

Equations (1) are derived on the assumption 
that the local field at a molecule can be repre- 
sented by the Lorentz expression : 


E,voc=E+42rP/3=3(€+2)E. (3) 


It is further supposed that the orienting effect 
of this field on the molecule of dipole moment y is 
opposed by a viscous force proportional to the 
angular velocity of rotation of the molecule. For 
spherical molecules of volume V in a liquid of 
Viscosity 4 the parameter 7 is then given by! 


t=3nV/kT, 


I being the absolute temperature. 
LS 


*P. Debye, Polar Molecules, Chap. V. 


The Clausius-Mossotti expression for ¢€ may 
be written 


e—1 e€—1 
e+2 €.+2 


where B=47Ny?/9kT and N is the number of 
molecules per cc. It is well known that this ex- 
pression, which is derived on assumption of the 
validity of (3), is not an adequate representa- 
tion for polar media.? In particular, if induced 
polarization is neglected (4) predicts an infinite 
value of ¢ for 8=1. Such values of 8 are found at 
temperatures of the order of room temperature or 
higher for substances with a large dipole moment 
but such a behavior of ¢ is, needless to say, not 
found experimentally. For nonstatic fields, 
Egs. (1), although of the proper form to account 
for many of the observed facts, require in many 
cases unreasonable values of V in the parameter 
7 to give the characteristic frequency region of 
dispersion and absorption correctly, especially 
for liquids of high dielectric constants. 

Debye and his co-workers have achieved con- 
siderable success in explaining these discrepan- 
cies by assuming that the dipole axis of a given 
molecule in a liquid is bound to an axis defined by 
the configuration of its neighbors by a ‘“‘hinder- 
ing”’ potential —yF cos 3’ where F is a potential 
energy and #’ the angle between the axes.® 

Van Vleck ? has, however, shown on the basis 
of statistical-mechanical calculations of dipole- 


2 J. H. Van Vleck, J. Chem. Phys. 5, 556 (1937). 
3 P. Debye, Physik. Zeits. 36, 193 (1935); P. Debye and 
W. Ramm, Ann. d. Physik 28, 28 (1937). 
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logz 





Fic. 1. Real and imaginary parts (¢’, e’’) of the dielectric 
constant plotted against the frequency dependent variable 
z. Solid curves are for the approximate equations, dashed 
curves for the Onsager theory as explained in the text. The 
scale for ¢”’ is twice that for ¢’. 


dipole coupling that the local field as given by 
(3) is not an adequate representation of the state 
of affairs. Onsager ‘ has obtained an expression 
for the local field on the view that the molecule 
in virtue of its electric moment has a modifying 
effect on the surrounding medium. This expres- 
sion Van Vleck has shown to be in agreement with 
his results to a second approximation in the di- 
pole-dipole coupling (the extent to which the 
calculations were carried). This agreement is, 
however, obtained only if one assumes in the 
dipole calculations that the molecules surround- 
ing a given one may be represented as a continu- 
ous distribution. Van Vleck has also shown that 
the excessive saturation curvature of ¢ as a func- 
tion of field strength found by using (3) is not 
predicted by the Onsager field. In view of all this 
it has seemed worthwhile to see what results are 
given by the Onsager theory as applied to the 
dielectric constant of polar liquids in radio- 
frequency fields. 


II. DISPERSION AND ABSORPTION FOR THE 
ONSAGER LOCAL FIELD 


A. Derivation ——Onsager’s formula may be 
written 


e—1 €—1 
e+2 €.+2 
where f= 3¢e(€.+2)/(2e+€.)(e+2). This formula 








4L. Onsager, J. Am. Chem. Soc. 58, 1486 (1936). 


COLE 


is derived on the assumption of a continuous 
medium in which a given molecule occupies a 
spherical cavity. For a more realistic picture of 
the environment as a discrete assemblage agree- 
ment with Van Vleck’s results can be obtained in 
an empirical way by multiplying the right side of 
(5) by a factor g (0< q< 1) as suggested by him in 
the paper referred to. Such a procedure can 
hardly be justified in any rigorous manner at the 
present time, but introduction of such a factor 
does not complicate matters unduly, and it is of 
interest to examine the behavior of the values of 
q required to obtain experimental agreement with 
(5) thus modified as will be done in Section III. 

In Debye’s analysis of the effect of a nonstatic 
field he shows that one obtains a complex factor 
multiplying the moment of the form 1/(1+7w7). 
Exactly the same considerations apply in the 
Onsager case and the net result is simply that 
the terms involving 8 are multiplied by 
1/(1+iw7). The dielectric constant is then com- 
plex : «= «’ —ie’’. As it stands (5) is quadratic in « 
and an exact solution is too unwieldy to be readily 
handled. However, if the quantity 3¢«/(2e+..) in 
f be replaced by $ the equation is linear in e: 


e—1 En — 1 3(€a+2) : 
“i -a=q45—— 1), (6) 
e+2 €o+2 2(€+2) 








and the essential difference between the Onsager 
and Lorentz fields is retained. Making this ap- 
proximation, inserting the complex modifying 
factor and solving for ¢’, e’’, one finds: 


€’ = €a + (€0— €0)/(1+27), 
e’ = (€9— €x)2/(1+27), 
where 
z=wr/(1—84), 


5=38(€.+2)(1—g). 


The curves for ¢’, «’’ are plotted against log 2 
in Fig. 1 and for comparison the curves obtained 
by a rigorous solution of the Onsager equation 
with ¢9=10, €.=2 are indicated by dashed lines. 
The error in the approximate analysis is ob- 
viously not serious even for this small a value 
of €0. 

B. Comparison with experiment.—It must be 
remembered that in applying the Debye Eqs. (1) 
the value of ¢ is taken from the experimental 
data on the assumption that it satisfies the 


(8) 
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Clausius-Mossotti relation (4). This is, of course, 
not the case for polar substances, a fact which 
must be borne in mind when the results of the 
hypothesis of hindered rotation are compared 
with the present work. Since ¢9 does satisfy the 
modified Onsager Eqs. (6) derived on plausible 
assumptions, this latter relation is evidently 
preferable as a basis for consideration of disper- 
sion and absorption. 

Equations (7) are quite general expressions for 
é, ein that they contain the results of the Debye 
theory and the unmodified Onsager theory as 
special cases for g=0 and 1, respectively. The 
various cases are seen from Eqs. (7) and (8) to 
differ only in the form of the variable z defined in 
(8). From Fig. 1 the center of the region of dis- 
persion and absorption occurs for z=1. Writing 
w=wmax for this case one has from (8) : 


(a) (€04+2)wmaxt/(€0.+2)=1, g=0, 

Debye, 
(b) wmax7/(1—6)=1, Modified Onsager, (9) 
(c) wmaxt=1, g=1, Onsager. 


Most liquids for which data on €9 over a wide tem- 
perature are available require values of g between 
0.8 and 1.1 (see Section III) if Eq. (6) is to fit 
the facts. The corresponding values of 6 do not 
depend markedly on the temperature and lie be- 
tween 0.3 and —0.15. For the alcohols 6 has 
values between 0.7 and 0.9 but the values of q 
are strongly temperature dependent as discussed 
in Section III. It is seen from these values of 6 
and the fact that (€9+2)/(€.+2) is a large num- 
ber that the values of 7 found for a given wmax 
from the Debye case (a) of (9) are considerably 
smaller than those required by (0d) or (c). Ex- 
perimentally, the values of 7 for polar liquids 
found assuming the Debye theory are much too 
small for liquids of large static dielectric constant 
. That is, they require, on the assumption of 
viscous damping, values of the molecular volume 
V to which 7 is then proportional, of a smaller 
order of magnitude than those inferred from in- 
dependent data, e.g. the van der Waals volume 
for gases. The prediction of larger values of 7 
and hence V on the modified Onsager basis is thus 
in the right direction. 

In general the experimental verification of (1) 
has consisted in seeing whether for a reasonable 
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value of the molecular volume V, assuming 7 
known, the data for ¢’, e’’ as a function of w and 
T can be reproduced in their general behavior. It 
is unfortunate that almost the only nearly ade- 
quate data on dispersion and absorption of simple 
liquids are for the alcohols which are abnormal in 
their static dielectric properties and are other- 
wise known to be highly associated. Debye in 
his book finds very good agreement between the 
data of Mizushima on n-propyl alcohol for three 
different frequencies over a temperature range 
20° to —60°C and calculated values of ¢’, ¢’’ 
from Eq. (1) using a value V=34X10-*4 cm’ 
(radius 2.0 10-8 cm). As good agreement can be 
obtained from (7) and (8) using the value V = 60 
X 10-** cm’ and the values of g required to give 
€9 correctly in (6) for the static case. 

Schmelzer ° has made absolute measurements 
of absorption for several alcohols at a single fre- 
quency and temperature and calculated the 
value of wax on the basis of an expression for e’’ 
as in (1) or (7). The values of V required on the 
basis of formulas (9) are given in Table | 
together with those deduced from the van der 
Waals equation of state. 

A detailed agreement with other estimates of I” 
is hardly to be expected when one considers the 
crudity of assuming viscous damping of spherical 
molecules but in view of how well such an hy- 
pothesis works for the problem of ionic mobilities 
one should expect a reasonable degree of correla- 
tion with values of V found in other ways. It is 
seen that the derivation with the modified local 
field does lead to quite reasonable estimates. 

C. Comparison with the hypothesis of hindered 
rotation.—Debye and Ramm ? in their paper on 
the effect of hindered rotation consider three 


TABLE I. Molecular volumes from the various assumptions 
about the local field and from the van der Waals coeffi- 
cient (allX 10™*4 cc). 











Mobpi- 
FIED VAN DER 
SUBSTANCE DEBYE | ONSAGER} ONSAGER| WAALS 
Water | V=0.5| 4.1 12 13 
Methyl Alcohol 16 38 132 28 
Ethyl Alcohol 16 30 111 35 
n-Propyl Alcohol 40 49 227 42 
Buty! Alcohol 72 76 350 47 




















5 W. Schmelzer, Ann.’d. Physik 28, 35 (1937). 
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cases : 


(a) The hindering energy F < kT, 
(b) The hindering energy F is comparable with kT, 
(c) The hindering energy F > kT. 


The effect in case (a) is to replace the quantity 
1/(1+7w7) multiplying 6 by a series of terms of 
the form: 


1/(1+¢wr/dn), (10) 


where the \, are constants. These reduce for w 
small to 


[(1—y?/9) —iwr(1 —11y?/54) ], 


where y= F/kT. The absorption is thus appar- 
ently reduced by the factor (1—11y?/54). How- 
ever, y is to be so taken that the difference be- 
tween ¢9 and ¢. is represented by (1—y?/9)8 
when one uses the correct value of the dipole 
moment yu in B rather than a spurious one. (The 
latter is implicitly done if one accepts the 
Clausius-Mossotti Eq. (4) as an expression for 
€o.) The absorption equations are written in terms 
of ¢9—€. rather than 6 (Cf. Eqs. (1) and (7).) 
Hence the absorption written in this way is re- 
duced by a factor (1—11y?/54)/(1—y?/9) which 
has obviously much less effect than the factor one 
might at first expect to be correct. 

From the series of terms in (10) one should 
expect a series of absorption regions and as a 
result a composite absorption curve somewhat 
broader, displaced to higher frequencies and with 
a slightly smaller maximum. This type of absorp- 
tion is actually found for many complex sub- 
stances but the curve is still much broader in 
most cases than such a series of terms seems 
capable of accounting for and in addition is in 
many cases symmetrical about the maximum on 
a logarithmic scale of frequency. It is difficult to 
see how this latter behavior can be given by terms 
of the form (10) with the constants \, varying in 
a reasonable way. There seems to be no evidence 
that the maximum of absorption is less or the 
frequency spread greater for simple liquids where 
we can have the greatest hope for the adequacy of 
either picture. If the region of dispersion and ab- 
sorption is merely shifted on a frequency scale as 
compared to the prediction of the original Debye 
theory with a reasonable value of 7 the explana- 
tion proposed here would seem to be superior. 
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and water at different temperatures plotted against 7/T, 
(=9kT/4x Np’). 














The foregoing remarks apply only if the hinder- 
ing energy F is small compared with kT. For 
many polar liquids much larger values of F are 
necessary to compensate for the inadequacy of 
the Clausius-Mossotti formula. Unfortunately, 
in the intermediate case that F is comparable 
with kT the series development of Debye and 
Ramm breaks down and no definite conclusions 
can be drawn. 

If F>>kT the effect of temperature is small 
compared with that of F. Debye and Ramm show 
that in this case kT is replaced by F/2 where it 
appears in 8 and 7. F is again to be so taken that 
this modified 8 gives ¢€9 correctly in Eq. (4) with 
the true value of the moment yp. The maximum 
value of the absorption is the same as that given 
by the original equation with the experimental 
value of €9; the only difference is that there is now 
a justification for using this value. However the 
absorption has a different frequency dependence 
because of the new significance of r. In particular 
a given value of the frequency of maximum ab- 
sorption and hence 7 requires a larger value of 
the molecular volume V if 7 is given by 3n V/(F, 2) 
rather than 3nV/kT. This is, as already discussed, 
in the right direction. If one computes JV’ for 


strongly polar liquids with an appropriate F the- 


values are found to be comparable with those 
obtained on the modified Onsager hypothesis but 
in general slightly smaller. The latter is thus 
quite as satisfactory in the way of quantitative 
agreement and in addition there is no need to 
account for a large hindering potential F. 

The inadequacy of the experimental data for 
simple polar liquids makes a more detailed discus- 
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Fic. 3. The local field parameter g as computed from Wy- 
man’s data at room temperatures plotted against T/T... 


sion hardly worthwhile at present. Enough has 
been said to show that the Onsager hypothesis 
leads to results which have the right trend from 
the original Debye theory and seems to work 
quite as well as Debye’s hypothesis of hindered 
rotation. This lends additional confirmation to 
the view that a great deal of the difficulties with 
Eqs. (1) and (4) is due to the failure of Eq. (3) 
for the local field in polar media. A survey of the 
experimental data emphasizes the necessity for 
measurements over as wide a frequency and tem- 
perature range as possible, detailed consideration 
of isolated data being comparatively useless as 
the above discussion makes evident. 

In the case of solids the interpretation of the 
relaxation time in terms of an ordinary viscosity 
is obviously impossible. One is then forced to 
introduce the relaxation virtually as an empirical 
parameter in the absence of a reasonable theory 
and there is therefore no such criterion for reason- 
ableness of values found for 7 as exists in the 
case of liquids. 


III. THe Locat Fretp IN THE Static CASE 


A. The temperature dependance of q.—It re- 
mains to consider briefly the values of the param- 
eter g required to give agreement of (6) with ex- 
perimental results. A value g=0 corresponds to 
the Lorentz local field and g=1 to the Onsager 
case. Intermediate values might be expected to 
simulate a local field of the Onsager type for dis- 
crete molecular configurations. In the case of 
cubic structures, agreement with the calculations 
of Van Vleck requires q~ }. 

One can hardly attach too much significance to 
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detailed behavior of the factor for specific cases 
but its general behavior should be of interest in 
indicating the extent to which dielectric con- 
stants of liquids can be explained in terms of 
dipole-dipole coupling. The significant parameter 
in such considerations is the temperature 
(8=4rNy?/9kT) and unfortunately there are 
comparatively little data for dielectric constants 
of polar liquids over a wide range of temperature. 
Perhaps the most complete data exist for various 
alcohols and water which are known to the chem- 
ist as highly associated and hence are not ideally 
suited for the present purpose. One might expect 
q to be more or less independent of 7, its value 
to be between 3 and 1. Actually g is for these 
substances strongly temperature-dependent. 
However, as pointed out to the writer by Profes- 
sor Van Vleck, one should expect a law of corre- 
sponding states if dipole-dipole coupling is the 
important factor. That is, the values of g as a 
function of a reduced temperature 7/7. should 
all lie on the same curve. Here 7, is the ‘‘Curie”’ 
temperature, having the value 47Ny?/9k if in- 
duced polarization is neglected.? How well such a 
condition holds for water and the alcohols is 
seen in Fig. 2. This rather striking behavior sug- 
gests the possibility of a gradual transition from 
coupling of the dipoles to the free rotation of 
Debye’s original theory as the temperature in- 
creases. However, it is unwise to generalize too 
much on the basis of data for such a limited num- 
ber of abnormal substances and in addition it 
should be remembered that the dipole moment yu 
is very nearly the same for all these liquids, being 
due primarily to the OH group. The behavior 
for acetone and nitrobenzene, also plotted in 
Fig. 2, serves further to emphasize this need for 
caution. These comprise the substances for 
which the International Critical Tables give 
fairly complete temperature data. 

B. The data of Wyman.—Wyman ° has col- 
lected an extensive set of data on electric polariz- 
abilities at room temperatures and dipole mo- 
ments from which it is possible to compute gq. 
Such values have been calculated for various re- 
lated chemical series in his table for which the 
data seem reasonably reliable and they have 
been plotted against 7/7. in Fig. 3, which thus 


5 J. Wyman, Jr., J. Am. Chem. Soc. 58, 1482 (1936). 
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shows that the required values are in general not 
at all unreasonable on the Onsager picture. It 
would be highly desirable to have more complete 
data on the dielectric constants of polar liquids 
as a function of temperature to see whether such 
a behavior as indicated by Fig. 2 for the alcohols 
holds for other series of related substances. Such 
data should be very valuable in any attempt at 
understanding the electrical properties of polar 
liquids generally. 

C. The local field and the work of van Arkel and 
Snoek.—Another line of evidence which lends 
confirmation to the modified Onsager hypothesis 
in regard to the local field is found in a study by 
van Arkel and Snoek ’ of the experimental data 
for solutions of polar liquids in nonpolar sol- 
vents. They found that for a considerable num- 
ber of ‘‘normal’’ polar liquids and over a wide 
range of concentrations the dielectric constant 
can be represented by a formula of the form: 


e—1 e«—1 4rN we 
€+2 Cot2 3 3kT+(40C/3)NiE 





C being a constant characteristic of the polar 
liquid (V=number of polar molecules per cc) 
having values between 1 and 2. In the notation of 
the present work this may be written: 


e—1 e€—1 B 


_ == , (11) 
e+2 e.+2 1+CB 





Now if one rewrites the modified Onsager formula 
(6) in terms of 8 on the right side, it assumes the 
form (taking ¢, to be the same for solute and 
solvent) : 


“1 a~t 2+¢ 


ai — B, 
€+2 €o+2 2+9(€.2+2)p 





' which is seen to be very similar to the expression 
(11). Namely, (12) differs from (11) only by 
having 2+4q instead of 2 in the numerator on the 
right. Furthermore, the values of C found to work 
in (11) over a considerable range of concentration 
are such as to make the factor on the right agree 
quite well with the right side of (12) for most of 
the cases considered by these authors. Most of 


7A. E. v. Arkel and J. L. Snoek, Trans. Faraday Soc. 30, 
707 (1934). 
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the data given are at one temperature or over 
only a limited range so the difference in functional 
dependence on 8 of the two expressions cannot be 
used to decide which of them is preferable. Here 
again it would be very desirable to have measure- 
ments at different concentrations made over as 
wide a temperature range as possible because of 
the dependence of 8 on T. 

D. Dielectric constants of solids——It should 
also be of interest to examine the values of g 
required to account for the dielectric properties of 
solids at various temperatures. There is, however, 
one thing to bear in mind: measurements of the 
dielectric constant as a function of T for solids 
have usually been made at a low frequency but 
one which is kept the same at all temperatures of 
measurement. Under such conditions many sub- 
stances show a dielectric constant decreasing 
quite rapidly with 7,° a behavior totally inex- 
plicable in terms of the local field. However, it 
should be remembered that the relaxation times 
found experimentally for solids are very long as 
compared to those for liquids and may be ex- 
pected to become increasingly so at lower tem- 
peratures. Thus the characteristic region of dis- 
persion of the dielectric constant and absorption 
occurs at lower and lower frequencies and the di- 
electric constants found experimentally as above 
may be very far from the static dielectric con- 
stants at low temperatures, particularly if this 
region of dispersion is much broader than that 
predicted by the simple theories. Because of this 
difficulty and the incompleteness of dispersion 
and absorption measurements as a function of 
frequency no attempt has been made in the 
present work to carry out an analysis of the data 
for solids. It would be highly interesting to see 
how much of the apparent decrease is real and 
how much is due to the relaxation time by means 
of dispersion and absorption measurements over 
a wide frequency range. 


IV. CONCLUSION 


From the considerations above it appears rea- 
sonable to conclude that the difficulties found in 
attempts to apply the Clausius-Mossotti equa- 
tion and the Debye equations for dispersion and 
absorption to measurements on dielectric media 


8 See, for instance: R. H. Fowler, Statistical Mechanics, 
2nd Ed., p. 824. 
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can be traced in large part to inadequacy of the 
Lorentz expression (3) for the local field. The 
hypothesis of hindered rotation then seems much 
less convincing as an explanation of the general 
behavior of such substances unless one considers 
it as an expression of the dipole-dipole coupling in 
the fashion that the Onsager local field appears to 
be. This latter with the modifying factor q is 
seen from several points of view to be a very satis- 
factory representation of the general behavior of 
dielectrics when one considers the crudeness of 
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some of its basic assumptions. Because of the im- 
portance of temperature and frequency as param- 
eters affecting the electrical properties of polar 
substances it would seem best to make a more 
detailed, critical study only when measurements 
are available over a wide range of these quanti- 
ties for a reasonable number of different sub- 
stances. 

The writer is indebted to Professor J. H. Van 
Vleck, who suggested the problem, for his 
interest and helpful advice. 
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The equation for the fluidity of a liquid based on the 
reaction rate theory of viscous flow is modified to interpret 
the viscosity data for liquid films. The free energy of 
activation for viscous flow in the surface films, as calculated 
from the viscosity data of Harkins and Meyers, and 
Langmuir and Schafer, is found to be about twice that for 
flow in the liquid in bulk. This may be explained by the 
occurrence of various types of association in the liquid 
films, e.g. lateral association of the surface molecules with 


HE general equation for the fluidity of a 
liquid has been formulated as follows: ! 
AV 


kT 2 (eae) 
a—— <0 =i CLD 
kT 





fh fh 


—AFt = (Aed3A/2)f 
— ( kT )| “ 





\=distance between equilibrium positions in 
direction of flow, 

\i=perpendicular distance between adjacent 
molecules in direction of flow, 

he=distance between adjacent molecules in 
direction of flow, 

\3= distance between molecules in plane of flow 
and normal to direction of flow, 


'Eyring, J. Chem. Phys. 4, 283 (1936). Eyring and 
Ewell, ibid. 5, 726 (1937). 


one another and vertical association with the solvent 
(solvation of the polar groups). Equations are given for 
the pressure and temperature effect on surface viscosity. 
It is shown how the knowledge of these effects may be 
used to interpret the structure of the films, and to separate 
and evaluate the contributions to the free energy of 
activation due to lateral association, solvation and hydro- 
gen bonds. 


f=applied external force per square centi- 
meter, 
AV=relative velocity of adjacent layers a dis- 
tance \, apart, ; 
AFt=free energy of activation for viscous flow 
per molecule. 


Recently, Harkins and Meyers,’ and Langmuir 
and Schafer * have determined the absolute vis- 
cosities of films of the liquid type. Whatever the 
mechanism for viscous flow in the surface layer, 
it can be treated by the same general reaction 
rate theory, since it is a process involving an 
activation energy. The three-dimensional equa- 
tion may thus be easily modified to interpret 
these new data. 


2 Harkins and Meyers, J. Chem. Phys. 5, 601 (1937); 6, 
53 (1938). 

3 Langmuir and Schafer, J. Am. Chem. Soc. 59, 2400 
(1937). 
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TABLE I. 
AF t For AFt For AFt (FILM) / 
CoMPOUND Liguip FILms | Pure Liquips AFt (LIQ.) 





Myristic acid | 10360 cal.®| 5770@70°C 1.80 


Palmitic acid 10900 5920" * 1.84 
Stearic acid 11030 6300“ “ 1.75 
Oleic acid 10360 6220@20°C 1.67 
Cetyl alcohol 15300 5800@60°C 2.62 

















The shear force acting on a single molecule at 
the interface is f\2 and it acts through a distance 
\/2. Thus in the two-dimensional case 


kT —AFI+(A/2)dof 
AV= —| exp ( ) 
h kT 


—AFt—(n/2)reo 
—exp ( : was *)}. (2) 
kT 








Therefore the fluidity 


kT rr —) 
pe ariitl exp (—) —exp ( ~) ; 
difh 2kT 2kT 





This, on expansion, when Af\»<kT, becomes 


AAAz2 
o= e-AFt/kT 
Ah 





To a first approximation \,;=), and 
$= (a/h)ersrviet (3) 


where a is the cross-sectional area of the units 
which flow. Although the correctness of our gen- 
eral equations (2) and (3) does not depend on 
the exact mechanism of flow, we may most easily 
interpret the viscosity of liquid films by assuming 
the same model as was used in the case of the 
three-dimensional liquid. Therefore we picture 
the molecules oriented at the surface as rolling 
around one another whenever there is a neighbor- 
ing hole in the surface which permits this 
rotation. 

Using the available viscosity data we have cal- 
culated AFt for the viscous flow of the liquid 
films, and compared the values with those for 
the pure liquids obtained from the three-dimen- 
sional equation and viscosity data available in the 
literature. We have assumed in the calculation 
that a is the cross section of a single molecule; 
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doubling the value of a increases AF} by only 
about 300 cals. The results are given in Table I. 

Langmuir and Schafer found that the surface 
viscosities showed large variations with the age 
of the water on which the films were spread, the 
age of the film, the pH of the water, etc. With 
new water and newly spread monolayers, the 
surface viscosity of stearic acid was found to 
vary from 0.002 surface poises at pH 5.8 to 2.6, 
at pH 9.0. This corresponds to a variation in AF} 
from 12,100 cal. to 16,500 cal. 

It is evident from these considerations that the 
free energy of activation for viscous flow in the 
surface layer is about two to three times that for 
flow in the liquid in bulk. In viscous flow, as in a 
chemical reaction, the process will go principally 
by the path requiring the least expenditure of free 
energy. Since the units of flow in the liquid are 
probably single molecules, there is evidently 
some degree of association in the liquid films. 
This may be due to lateral association of the sur- 
face molecules to give dimers, trimers, etc., or to 
vertical association or solvation of the polar 
groups. In either case the free energy of activa- 
tion required for the flow process will be greater 
than it would be if the flow units were single 
molecules. 

The widely varying viscosities found by Lang- 
muir and Schafer suggest that no single explana- 
tion will cover all the observed data. Thus the 
high values of AFt may be due to a combination 
of such factors as vertical association, lateral as- 
sociation, and salt formation. The picture is akin 
to the situation in reaction kinetics in a system 
where there are a great many possible simultane- 
ous reactions, with widely different reaction 
rates. The viscous flow process may be con- 
sidered as a bimolecular reaction, consisting of 
the collision of a pair of molecules, rotation of 
the activated complex formed, and then dissocia- 
tion of the complex. 

If there are two different molecular species 
present, e.g. double molecules and single mole- 
cules, we may formulate the fluidity for the 
binary system as follows 


n; N\N2 n? 
= ¢1— +2612 +¢r—. (4) 


No" NoiN 02 Nor 





Here ¢; is the specific fluidity for the rotation of 
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two single molecules, ¢12, for the rotation of a 
single around a double molecule, and ¢2, for the 
rotation of two double molecules; 2; and nz are 
the number of molecules of each species per cc in 
the mixture, and mo; and mo2 are the numbers per 
cc in a system containing only the one type of 
molecule. We see that the presence of large mole- 
cules will decrease the fluidity by cutting down 
the number of pairs of small molecules which are 
able to rotate. It is to be noted, however, that 
even if a large proportion of the surface molecules 
have such a large AF} as to be practically im- 
mobilized, the over-all fluidity, ¢, may still be 
appreciable since it is the processes of low AFt 
that mainly determine the reaction rate. In other 
words ¢; is much larger than @¢y2 or ¢o. 

We may predict the pressure effect on viscosity 
as in the three-dimensional case. If 6s is the area 
of the surface hole necessary for viscous flow, the 
free energy required to make a hole this size in 
the surface is (p+ ;)6s, per molecule. Here p is 
the applied external pressure; p; is the internal 
pressure, defined as (0E/ds)r where s is the sur- 
face area; p and ; are actually pressures per unit 
length per film thickness. Because film thickness 
cannot usually be varied arbitrarily, this dimen- 
sion is customarily ignored in specifying surface 
forces and energies. If we put Nés=¥%, the free 
energy of making one mole of holes becomes 
(p+pi)=. 


We may now write 
AFt=RT In f+ (p+ p:)2+AF st. (5) 


Here f is the partition function for the transla- 
tional degree of freedom of the normal liquid 
which, in the activated state, is included in the 
factor RT/h. RT Inf amounts to about 1 kcal. 
for most substances. AF,t is the structural free 
energy of activation. For nonpolar molecules 
itis negligible, but for systems having directional 
bonds, e.g. hydrogen bonds, the satisfactory 
saturation of these bonds in the activated state 
may be prevented because of their short range, 
directional nature. 

It was shown in earlier treatments for liquids 
that p;V (V=Nz», where v is volume of the hole 
required for flow process) was approximately one- 
third the energy required to vaporize the unit 
which flowed. Since the energy required to make 
a hole the size of a molecule equals the energy of 
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vaporization, this was interpreted to mean that a 
hole one-third the size of the molecule is required 
for viscous flow. This picture then led auto- 
matically to an explanation of the pressure effect 
on viscosity. If we assume that 9, is one-third 
the energy required to vaporize the flow unit 
from the surface, we may apply this same treat- 
ment to the liquid films. Eq. (5) thus becomes 


AFt=RT Inf+AE,/3+AF.t+p>. (6) 


We might also use the well-known thermo- 
dynamic relationship p+;=(a/8)T, where a 
and £ are the coefficients of thermal expansion 
and compressibility of the liquid films, respec- 
tively. We then obtain from Eq. (5) 


AFt{t=RT In f+(a/B)T=Z+AF,t. (7) 


Eq. (7) suggests that if we knew the variation of 
a and 8 with temperature and pressure, we should 
have a method of predicting the temperature and 
pressure variation of the viscosity, provided that 
changes in AF, were negligible. 

The only experimental measurement of the 
variation of surface viscosity with pressure is that 
of Talmud, Lubman, and Suchowolskata.t They 
found well defined maxima in the pressure vs. vis- 
cosity curves, which they attributed to a varia- 
tion in the degree of hydration with pressure. 
Such a variation is likewise indicated by Lang- 
muir’s theory of the structure of liquid films. 
Since these authors gave no absolute viscosities, 
we were unable to calculate AFt from their data. 
We shall see below however that we can estimate 
the effect of hydration on AFt, so that a further 
determination of the variation of viscosity with 
pressure might enable us to determine whether or 
not hydration varies with pressure. 

A knowledge of the temperature coefficient of 
the surface viscosity would enable us to separate 
the entropy and energy terms in AFt. This would 
furnish a method of investigating further the 
mechanism of the flow process. Any mechanism 
tending to decrease the number of molecules 
which flow readily would result in a decrease in 
the entropy of activation. A change in the type 
of molecule or of its interaction with the solvent 
would be evidenced by a change in the energy 
term. 


4 Talmud et al., Zeits. f. physik. Chemie A151, 401 (1930). 
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Although we do not know the temperature 
variation of the surface viscosities we may 
roughly separate the entropy and energy terms 
by the following approximations. From the tem- 
perature coefficient of the viscosity of liquid 
stearic acid,® we find that at 70°C, AFt=6300 
cal., AHt=3580 cal.; and ASt=8 e.u. It is 
probably safe to assume that the principal differ- 
ence in the AFt terms in liquid in bulk and liquid 
film is in the terms for AHt. Making this assump- 
tion we find that AHt¢ for the film is 4830 cal. 
greater than for the fluid in bulk. This difference 
would be explained if each stearic acid molecule 
flowing in the surface had to drag with it one 
water molecule, since AHt for the viscous flow of 
water is about 4800 calories.*® 

Since the chief factor determining the viscosity 
of unimolecular films is A//t, it is of interest to 
consider another method of estimating the con- 
tributions which make up this quantity. If the 
unit which flows is a single acid molecule with 
one or more solvent molecules attached, the 
energy required to make a hole can be written as 
the sum of two terms. 


AHt=AH;{+ Apt. 


The term AH;jt is the energy required to make a 
hole in the fatty acid surface,, while AHot is the 
energy required to extend the hole into the 
solvent. We expect AM; to be approximately one- 
third the heat of vaporization of a liquid film toa 
gaseous film. We have already seen that AHot 
can be estimated from either the viscosity or heat 
of vaporization of the solvent. 

On the basis of his duplex film theory Lang- 
muir found that the two-dimensional vapor pres- 
sure of films of myristic acid on dilute HCI obeyed 
the Clapeyron equation.’ He gives the following 


5 Bingham and Fornwalt, J. Rheology 1, 372 (1930). 
6 Eyring and Ewell, reference 1, p. 733, Table 6. 
7 Langmuir, J. Chem. Phys. 1, 770 (1933). 
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equation, calculated from Adam’s data: 


logio p’ = 8.058 —1950/T. 


From this we find that AE,,, of liquid to gaseous 
film=8750 cal. Therefore AH,t for myristic 
acid = 8750/3 or 2920 cal. Considering the ap- 
proximations involved, this agrees reasonably 
with the 3580 cal. found for the energy of activa- 
tion for viscous flow of stearic acid. 

It should be emphasized that the gaseous film 
in equilibrium with a liquid film cannot be con- 
sidered as analogous in all respects to a three- 
dimensional gas-liquid system. Since even for 
gaseous films there probably is some activation 
energy for viscous flow due to the solvated polar 
groups, we shall not, in general, have the situa- 
tion in the real gas, where viscosity increases with 
temperature. 

We have distinguished in the case of the liquid 
films two types of association : lateral association 
of the fatty acid molecules with themselves, and 
vertical association between the polar groups and 
the underlying water molecules. Thermodynamic 
consideration of the equation of state for the uni- 
molecular films measures the lateral association. 
Kinetics provides a method of measuring both 
types of association. Thus a gaseous film will ex- 
hibit thermodynamic properties ‘similar to those 
of a real gas. In regard to rate processes, such as 
viscosity and diffusion, the gaseous films will 
probably have the properties of a liquid. 

We have outlined briefly the possible ap- 
proaches to the viscosity data for liquid films on 
the basis of the reaction rate theory of viscous 
flow. As further data become available, the ap- 
plication of such methods should tell us more of 
the actual mechanisms involved, and of the actual 
structure of the surface layers. 

We wish to acknowledge our indebtedness to 
Dr. Irving Langmuir and to Professor Hugh S. 
Taylor for helpful discussions on this subject. 
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The dielectric constant of methyl alcohol vapor has 
been measured at 14 temperatures between 25.0°C and 
206.0°C. At each temperature measurements were made 
at an average of 23 different pressures ranging from zero 
to near saturation. The slope, obtained by least squares, 
of the resulting (K—1)/(K+2) versus pressure curve 
gives one point on the Debye line. The 14 such points fall 
accurately on a straight line whose equation as obtained 
by least squares, and after making proper allowance for 
probable errors of calibration, is: [(K—1)/(K+2)]RT/p’ 
=(7.72+0.15)+(17,730+71)/T. Using recently estab- 
lished values of the Boltzmann gas constant and Avo- 
gadro’s number, this yields a value (1.706-+0.005) x 10718 
es.u. for the electric moment of the methyl alcohol 
molecule. Taking the sum of the electronic and atomic 


polarization as 8.6 leads to an electric moment of (1.688 
+0.006) x 10-18. The mean (1.698+0.005) x 10~!8 is taken 
as the electric moment of the methyl alcohol molecule. 
Data on ethyl alcohol previously published by Knowles 
have been recalculated, using a least squares treatment 
throughout, and making a minor correction to the value 
used for the absolute zero of temperature. The data 
conform to the equation [(K—1)/(K+2) ]RT/p’ =(13.06 
+0.20)+(17,725+87)/T, wherein probable errors of cali- 
bration have again been allowed for. This yields a moment 
(1.706+0.006) x 10-18, Taking the sum of electronic and 
atomic polarization as 13.6 leads to a moment (1.696 
+0.008) x 10-18. The mean (1.702+0.007) x 10-8 is taken 
as the electric moment of the ethyl alcohol molecule. 





ONTINUING an eariier reported! effort to 
fix more definitely the values of the electric 
moment of several of the lighter alcohol mole- 
cules, the author has made extensive dielectric 
constant measurements on methyl alcohol vapor. 
Dielectric constants were measured by the 
heterodyne beat method, using apparatus previ- 
ously described.2 The vapor condenser with 
quartz insulators, and the vapor system, were the 
same as those used in previously reported work 
on water vapor*® and isopropyl alcohol vapor.! 
Dielectric constants were measured at 14 different 
temperatures between 25.0°C and 206.0°C. At 
each temperature, measurements were made at 
an average of 23 different pressures, ranging from 
the lowest measurable up to near saturation when 
such saturated vapor pressure did not exceed one 
atmosphere. Two different samples of methyl 
alcohol were used, each of the highest purity 
obtainable and carefully dried. 

Values of (K —1)/(K+2) for a given tempera- 
ture were plotted to a large scale against p’, 
where ~’ represents the pressure the vapor would 
exert were it an ideal gas; the quantity p’ is 
obtained from the observed pressure p and the 
van der Waals constants a and b, thus: 


b’=p[1+(p/RT)(a/RT—D)]. 


'Stranathan, J. Chem. Phys. 5, 828 (1937). 


; Stranathan, Rev. Sci. Inst. 5, 315, 334 (1934). 
Stranathan, Phys. Rev. 48, 538 (1935). 


At the higher temperatures, for which the highest 
pressures used were far below saturation, the 
resulting curves are accurately linear over their 
entire range. At lower temperatures, slight 
curvature is apparent at the higher pressures, 
though a large scale plot is necessary to detect it 
except very close to saturation. The only use 
made of such plots was to determine the extent of 
linearity. Barely perceptible curvature became 
apparent at pressures approximately two-thirds 
saturation. For each temperature the slope of the 
best least squares line through all points on the 
linear section was found; the average number of 
points on these linear sections was 19. The several 
advantages of securing extensive data and of the 
least square treatment of these data have been 
discussed previously.' 
The molecular polarization is given by 


(K—1)M/(K+2)d 
=(K—1)RT/(K+2)p’=A+B_T, 


the last equality being a result of the Debye 
theory. The molecular polarization at each 
temperature was obtained from the slope of the 
(K—1)/(K+2) versus p’ curve for that tempera- 
ture. The accuracy of the data can be judged 
from the calculated molecular polarizations and 
their probable errors as shown in Table I. 
Figure 1 shows the products of molecular 
polarization and absolute temperature plotted 
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against 7. The data fall accurately on a straight 
line, no point being off by more that 0.4 percent. 
Previous data by Miles* and by Kubo’ are also 
shown. In calculating the best straight line 
through the 14 points of the present work, the 
question arises as to how the individual points 
should be weighted. While the author doubts 
that all variations are truly accidental, he feels 
that the largest are. Under such circumstances it 
appears advisable to weight each point inversely 
proportional to the square of its probable error. 
With such weighting, a least square solution for 
the best straight line through the author’s data 
gives: 
(K—1)RT/(K+2)p’ 

= (7.72+0.15)+(17,730+59)/T. 


Giving all points equal weights results in 
essentially the same equation, the constants so 
obtained being (7.55+0.13) and (17,797+52). 

According to the Debye theory, the permanent 
electric moment uw is given by: 


u=(9kB/4rN)}, 


where & is the Boltzmann gas constant per mole- 
cule and N is Avogadro’s number. Recent work of 
Bearden,* Kellstrém,’ and many others on a new 
value for the electronic charge leads to new values 
of N and k. These works have been so numerous 
and so convincing that there appears no reason 
for delaying longer the use of these new values in 
calculations of electric moments. Taking the gas 
constant per mole as (8.3136+0.0010) X10’ ergs 
per C°,® and N as (6.0230+0.0005) X10”? mole- 


TABLE I. Methyl alcohol. Experimental values of molecular 
polarization, with associated probable errors. Values 
of T are absolute temperatures in C°. 








MOLECULAR 
POLARIZATION 


53.633 40.073 
51.9344 .063 
50.2354 .037 
48.5654 .035 
47.055+ .039 
46.007+ .029 
44.735+ .030 


TEMPERA- 
TURE T 


MOLECULAR 


TEMPERA- 
a POLARIZATION 


TURE 7 





298.2 67.29 +0.25 386.5 
313.7 64.39 + .12 401.7 
313.7 64.47 + .11 417.1 
326.6 61.86 + .12 433.4 
331.2 61.32 + .18 449.0 
346.7 58.6844 .083 464.2 
367.4 55.9274 .097 479.2 




















4 Miles, Phys. Rev. 34, 964 (1929). 
5 Kubo, Inst. Phys. Chem. Res., Tokyo, Sci. Papers 26, 
242 (1935). 
° Bearden, Phys. Rev. 48, 385 (1935). 
a Nature 136, 682 (1935); Phil. Mag. 23, 313 
§ Birge, Rev. Mod. Phys. 1, 1 (1929). 
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300 340 380 420 460 500 


cules per gram mole,’ then k= (1.3803+0.0002) 
X 10-1; hence, 


u = (1.2811+0.0002) x 10-*°B?. 


This is larger than the value based upon con- 
stants previously accepted® by a factor 1.0069; 
hence all previous moments obtained from 
dielectric studies should be increased by this 
factor. While the difference is not great, it is 
definitely larger than the probable errors associ- 
ated with a number of moments. The electric 
moment of the methyl alcohol molecule obtained 
from the present data is (1.7058+0.0039) x 10- 
e.s.u.l° 

The probable errors quoted above are those 
calculated from accidental variations of observed 
data. Two other appreciable errors enter in the 
calibration of apparatus. Numerous and varied 
measurements of the effective capacity of the 
evacuated vapor condenser indicate that 1 part 
per 1000 would be a reasonable estimate of the 
probable error associated with this measurement. 
Similarly, an allowance of 2 parts per 1000 might 
be associated with the intercalibration of the 
large (1500 uuf) condenser used in measuring the 
effective capacity of the evacuated vapor con- 
denser and the small (15 uuf) condenser used in 
measuring the increase in capacity occasioned by 
the presence of the vapor. Error introduced 
through measurement of temperature is neg!i- 


9 Birge, Nature 137, 187 (1936). 

10 This moment was reported as 1.694107! at the 
Indianapolis Meeting of the Am. Phys. Soc. (1937). When 
corrected to the new values of constants used here this 
corresponds exactly to the value here reported. 
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ELECTRIC MOMENTS OF ALCOHOLS 


gible. Error due to adsorption of material on the 
insulator surfaces* may not be completely elimi- 
nated even by the procedure used here, though 
the residual must be small. Some evidence" 
indicates an appreciable error due to failure of 
the van der Waals gas law to represent accurately 
the density of the vapor in terms of its pressure. 
As far as is known, however, this error is quite 
small for measurements made at low pressures. If 
allowance is made for the probable errors of 
calibration, it is found 


(K—1)RT/(K+2)p’ 
=(7.7240.15) +(17,730471)/T 


and “u=(1.7058+0.0048) X10-'*. The present 
results are compared with those of other extended 
works on vapors in Table II. While results of 
dilute solution studies are in general agreement 
with those on vapors, they are not as reliable; 
therefore they are not included. 

The term A should be the sum of the electronic 
and atomic polarizations. A fair estimate" of this 
sum is a value 5 percent greater than the molecu- 
lar refraction for the Na D line. Hence, one ex- 
pects a constant term A about (1.05 X8.23) =8.6, 
whereas the experimental value found here is 
(7.72+0.15). The difference may be experimental 
error, though it is considerably larger than one 
would suspect from the associated probable error. 
As Groves and Sugden" point out, it may be pref- 
erable in the calculation of moments to use the 


TaBLE II. Comparison of results. A, B, and moments under 
(1) are calculated entirely from dielectric data. Moments 
under (2) are calculated by taking A =8.6 for methyl alcohol 
and 13.6 for ethyl alcohol. All moments are on the basis of 
the new values of Avogadro’s number and Boltzmann's gas 
constant, 








; u X 1018 
AUTHOR B (1) (2) 


Methyl Alcohol 








17,730+71 
17,430 
17,710 


Stranathan 
Milest 
Kubo 


7.72 40.15 
8.13 
9.54 


1.706 +0.005 
1.691 
1.705 


1.688 +0.006 
1.684 
1.724 











Ethyl Alcohol 





1.706 +0.006 
1.698 
1.708 
1.680 


Stranathan 
iowles!4 

Miles* 

Kubo 


13.06 +0.20 
13.6 


11.80 
14.25 


17,725 +87 
17,550 
17,780 
17,200 


1.696 +0.008 
1.697 


1.675 
1.691 

















Stuart, Zeits. f. Physik 51, 490 (1928). 
® Sugden, Trans. Faraday Soc. 30, 738 (1934). 
(193 . and Sugden, J. Chem. Soc. London, 158, Jan. 
/). 
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term A estimated from optical data. If A is taken 
as 8.6, then statistical weighting of data 
at different temperatures leads to a value B 
=(17,354+11) and u=(1.6877+0.0008) X 10-8. 
This appears much better determined, and the 
method is probably preferable for comparison of 
moments where calibration errors cancel out. 
However, if the same allowance as before is made 
for calibration errors, together with an allowance 
of 2 percent in the value chosen for A, we find B 
= (17,354+80) and u=(1.6877+0.0055) x 107". 
This moment is smaller than that found entirely 
from dielectric data by about twice the sum 
of the probable errors. The mean of these 
(1.698+0.005) X 10-18, may be taken as the most 
probable value for the moment of the methyl 
alcohol molecule. Results of similar calculations 
from data of other workers are shown in Table II. 

The data of Knowles“ on ethyl alcohol vapor 
appears sufficiently extensive and reliable to 
warrant a least squares treatment. Knowles 
obtained data at 10 different temperatures be- 
tween 24.55°C and 177.1°C. Observations at an 
average of 17 widely different pressures were 
made at each temperature. An average of 14 of 
these points fell on the straight line section of 
each curve. The author, having the original data 
available, has therefore made a recalculation, 
using least squares throughout, and making a 
minor correction for the previous use of the 
approximate value —273.0°C for the absolute 
zero of temperature. The resulting molecular 
polarizations with their associated probable 
errors are shown in Table III. 

Products of molecuiar polarization and abso- 
lute temperature are shown plotted against 
absolute temperature in Fig. 1. The 10 points fall 
accurately on a straight line, no point being off by 
more than 0.4 percent. The previous points 


TABLE III. Ethyl alcohol. Experimental values of molecular 
polarization, with associated probable errors. Values 
of T are absolute temperatures in C°. 








MOLECULAR 
POLARIZATION 


60.58 +0.23 
57.844 .11 
56.29+ .09 
53.02+ .16 
52.334 .34 


TEMPERA- 
TURE T 


MOLECULAR 
POLARIZATION 


72.30+0.46 373.3 
68.28+ .34 395.5 
66.394 .11 410.9 
64.214 .17 443.6 
61.834 .34 450.2 


TEMPERA- 
TURE T 





297.7 
319.7 
332.0 
347.9 
364.3 




















4 Knowles, J. Phys. Chem. 36, 2554 (1932). 
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obtained graphically by Knowles" from the 
same data, as well as previous data by Miles* and 
by Kubo are shown. If the present points are 
weighted statistically, a least squares solution for 
the best straight line through them gives: 


(K-1)RT/(K+2)p' 
= (13.06-0.20) +(17,725-+78)/T. 


This yields an electric moment (1.7056+0.0052) 
X10-'* for the ethyl alcohol molecule. Giving 
equal weights to all points leads to constants 
(13.35+0.23) and (17,610+85) in the above equa- 
tion, thus yielding a moment (1.7000+0.0058) 
X10-'8. The probable errors of calibration in 
Knowles work probably approximated those 
estimated in the present work on methyl alcohol. 
If the same allowance is made for them, we 
obtain from a statistical weighting of points: 


(K—1)RT/(K+2)p' 
= (13.06-+0.20) + (17,725487)/T 


from which u=(1.7056+0.0059) X 10-8. These 
results are compared with other extended works 
on the vapor in Table II. 

The constant term A =(13.06+0.20) is again 


LEROY AND A. R. GORDON 


slightly smaller than the value (1.05 12.93) 
= 13.6 which we might reasonably expect. Taking 
A = 13.6, statistical weighting of data at the vari- 
ous temperatures leads to values B = (17,517+8) 
and u= (1.6956+0.00055) X 10-8. If allowance is 
made for probable errors of calibration, together 
with a 2 percent probable error in the value 
chosen for A, then B=(17,517+110) and 
u=(1.6956+0.0076) X10-'8. This moment is 
smaller than that found entirely from dielectric 
data by somewhat less than the sum of the proba- 
ble error. The mean (1.702+0.007) X10-!8 may 
be taken as the most probable value for the 
moment of the ethyl alcohol molecule. Results of 
similar calculations from data of other workers 
are shown in Table II. 

Judging from the general agreement of results 
of various workers, together with the probable 
errors associated with these, it appears that little 
real improvement can be made in the moments of 
the methyl and ethyl alcohol molecules until 
several possible sources of error are eliminated. 
The author is attempting several improvements 
in this direction before proceeding to a study of 
other moments. 
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The Transference Numbers of Potassium Acetate in Aqueous Solution 


D. J. LE Roy Anp A. R. GoRDON 
Chemistry Department, University of Toronto, Toronto, Canada 


(Received April 25, 1938) 


The transference numbers for aqueous solutions of potassium acetate at 25°C have been 
determined: by the moving boundary method for the concentration range 0.01 N to 0.25 N. 
Two different types of cell and two different methods of measuring the number of coulombs 
were employed. Up to 0.16 N the transference numbers can be represented by the Shedlovsky- 
Longsworth equation within +0.05 percent, and lead to a transference number for the cation 
at infinite dilution of 0.6428, in close agreement with that obtained from the limiting ionic 


mobilities. 


N recent years the determination of the trans- 

ference numbers of ions in solution by the 
moving boundary method has been developed 
(thanks primarily to the work of MaclInnes, 
Longsworth and their associates) until it will 
yield results of the highest precision ; combined 
with accurate conductivity data, such numbers 
permit the calculation of limiting ionic mobilities 


with an accuracy of a part in several thousand. 
In the course of a study of the transference 
numbers for several electrolytes, Longsworth’ 
has determined the transference numbers for 
solutions of sodium acetate; in this paper we 
report the numbers for solutions of potassium 
acetate. 
1 Longsworth, J. Am. Chem. Soc. 57, 1185 (1935). 
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TRANSFERENCE OF POTASSIUM ACETATE 


Since the moving boundary method has been 
adequately discussed in the literature,” no general 
description is needed here. The cells used were 
of two types; the first (referred to as I) is essen- 
tially that of Smith* and is similar to that used 
by Reevely and Gordon ;‘ here a lithium acetate 
solution serves as indicator solution with a 
descending boundary. The electrode on the 
closed side (the cathode) was a platinum wire 
covered with fused silver chloride, and the 
general technique of filling the cell and forming 
the boundary was essentially theirs. In our 
measurements, the graduated tube (A in their 
Fig. 1) had an internal diameter of 4.3 mm. 
The second type (II) was a simplified autogenic 
boundary cell with rising boundary, illustrated 
in Fig. 1. The anode (of chemically pure 
cadmium) was carefully machined to fit the 
lower end of the graduated Pyrex tube, and 
before each experiment was sealed in the tube 
with De Khotinsky cement. To do this the lower 
end of the tube was heated by a small wire 
resistance bound around it, and the anode, 
whose sides had been covered with melted 
cement, was pressed into place; the lowest 
graduation on the tube is far enough from the 
heater so that there is no appreciable rise in 
temperature of the graduated part of the tube 
during the brief heating. The tube B was of 
sufficient length to ensure that none of the 
products of electrolysis at the cathode (a plati- 
num wire covered with fused silver chloride) 
could reach the graduated tube during a run. 
In filling the cell, solution was added at C, any 
air that might be trapped in the graduated tube 
being allowed to escape through a narrow 
capillary inserted at A, the capillary being 
removed after filling. Two cells of type II were 
used (internal diameter of the graduated tube 
3.4 and 2.5 mm). The volumes between gradua- 
tions in all three cells were determined by 
mercury weighings as recommended by Mac- 
Innes, Cowperthwaite and Huang.® The cells 
were mounted in a glass-sided water bath, 
electrically controlled to +0.01°C, and the 


2 MacInnes and Longsworth, Chem. Rev. 11, 171 (1932). 

*Smith, Nat. Bur. Stand. J. Research 6, 917 (1931). 
on and Gordon, Trans. Electrochem. Soc. 26, 261 

5 MacInnes, Cowperthwaite and Huang, J. Am. Chem. 
Soc. 49, 1710 (1927). 
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movement of the boundary was followed by 
the illuminated slit method recommended by 
MacInnes and Longsworth.? 

To determine the number of coulombs neces- 
sary to displace the boundary through a known 
volume, two entirely different methods of 
measurement were employed. The first of these 
involved the use of microcoulometers as de- 
scribed by Reevely and Gordon.‘ The cathodes 
were small platinum crucibles, and the anodes 
were silver wires surrounded by glass ‘shields 
perforated at the sides to permit the passage of 
the current; each coulometer contained 5 cc of a 
10 percent silver nitrate solution tested as 
recommended by Rosa and Vinal.* The only 
change in the technique from that of Reevely 
and Gordon was that the removal of the electro- 
lyte and wash water from the crucible after the 
experiment was effected by suction using a 
small filter stick of a type familiar in micro- 
chemical work; the stick was always weighed 
with its coulometer, and thus there was no 
chance of loss of minute silver crystals. Duplicate 
determinations of the number*of coulombs by 
this method as a rule are self-consistent to 0.03 


6 Rosa and Vinal, Nat. Bur. Stand. Bull. 13, 479 (1916). 
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percent or better when the total weight of silver 
deposited is from 10 to 40 mg; numerous trials 
have shown that such determinations agree 
within the same limit of accuracy with the 
result of current-time measurements. 

The second method is what might be called 
a semi-constant current device, which permits 
the number of coulombs to be measured from 
current time readings without employing the 
elaborate apparatus developed by MaclInnes, 
Cowperthwaite and Blanchard.’ The significant 
features are shown in Fig. 2; a somewhat similar 
circuit has been reported by Hartley and 
Donaldson.’ The 1 B 4 vacuum tube has the 
characteristic that the plate current is almost 
independent of the plate voltage; consequently 
when the tube is connected as shown, the current 
is practically independent of the internal re- 
sistance of the cell, which may increase by as 
much as 100 percent during the course of a run. 
The potential EZ, is supplied by 45-volt B 
batteries, and the filament current by a 2-volt 
storage cell; for satisfactory operation E£, can be 
varied up to 630 volts depending on the resistance 
to be expected in the cell. The screen voltage 
FE, is kept at 62} volts, and the grid bias re- 
sistance R is about 20,000 ohms; small C 
batteries E; are used to bring the grid back to 
2.5-4.0 volts negative. During an experiment, 
there is a slight but steady fall in current,® but 
this can be followed without difficulty by read- 
ings at 5- or 10-minute intervals; a simple 
graphical integration then gives the mean 
effective current during any part of the run. 
The current was measured by means of a 500- 


7 MacInnes, Cowperthwaite and Blanchard, J. Am, 
Chem. Soc. 48, 1909 (1926). 

( 8 ed and Donaldson, Trans. Faraday Soc. 33, 457 
1937). 

* For example, in one experiment with 0.06627 N solu- 
tion, the current when the boundary crossed the first of the 
division marks was 1.4668 ma, while when it crossed the 
last, the current had fallen to 1.4639 ma. 








ohm certificated standard resistance and a Leeds 
and Northrup type K potentiometer, while time 
readings to the nearest half-second were effected 
by an electrical chronograph. 

The potassium acetate solutions were made by 
essentially the method employed by MacInnes 
and Shedlovsky’ in preparing their sodium 
acetate solutions. B.D.H. “analar’’ potassium 
bicarbonate, recrystallized from conductivity 
water, was heated in a platinum crucible to 
constant weight at 350°C in an atmosphere of 
CO.; a weighed amount of the resulting normal 
carbonate was dissolved in water, a weighed 
quantity of glacial acetic acid, prepared after the 
method of Bousfield and Lowry," was added in 
slight excess of the amount needed to neutralize 
the carbonate, the solution was boiled to remove 
COs, and then diluted gravimetrically to form a 
stock solution of convenient strength. All solu- 
tions were prepared by gravimetric dilution of 
such stock solutions. It was found essential to 
evacuate any solution before placing it in the 
cell, but weighings before and after evacuation 











TABLE I. 
SOL. VoL. 

CELL "94 ttobs | Corr. Corr. ty tteale 

II |0.009963/0.6493 | 0.0007 |—0.0000} 0.6497 | 0.6498 
.6499 
.6478 

II |0.021849} .6525 | 0.0003 |—0.0001} .6529 .6527 
.6529 

II |0.036317| .6557 | 0.0003 |—0.0002} .6554 .6552 
.6548 

I |0.05003 | .6569 | 0.0002 }—0.0002| .6569 .6569 
.6565 
.6574 

II |0.06627 | .6592 | 0.0001 |—0.0003} .6588 .6585 
.6588 

I [0.09990 | .6613 | 0.0000 |—0.0003) .6606 .6609 
.6605 

II |0.10747 | .6618 | 0.0000 |—0.0006| .6610 .6614 
.6614 

I |0.16964 | .6643 | 0.0000 |—0.0004| .6641 .6639 
.6647 

I |0.2521 .6690 | 0.0000 |—0.0006} .6687 6652 
.6696 


























10 MacInnes and Shedlovsky, J. Am. Chem. Soc. 54, 1429 
(1932). 
1 Bousfield and Lowry, J. Chem. Soc. 99, 1432 (1911). 
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TRANSFERENCE OF POTASSIUM ACETATE 


measured any loss in water from this cause. 
The lithium acetate solutions, required for cell I, 
were made by the same procedure, the lithium 
carbonate being twice recrystallized from hot 
conductivity water. The conductivity water used 
had a specific conductivity of from 0.4 to 
1.0 10-*. 

The pH of the acetate solutions was always 
adjusted to 5.5—-5.8 (by the addition of weighed 
amounts of acetic acid after dilution if necessary) 
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since by so doing the effect of dissolved carbon 
dioxide is minimized; the pH was determined 
colorimetrically with brom-cresol purple or by 
means of a quinhydrone electrode. The con- 
centration of potassium acetate in the solution 
was calculated from its potassium content, and 
the effect of the slight excess of acetic acid was 
taken into account in computing the solvent 
correction; see below and also reference 10. In 
computing the volumetric concentration of the 
solutions from the gravimetric concentration, 
density data” from International Critical Tables 
were employed. 

Table I summarizes the results; tors is com- 
puted by the relation t:.»,=1000VC/f where C 
is the concentration of the solution in equivalents 
per liter and f is the number of faradays neces- 
sary to displace the boundary through V cc. 
When f was measured by means of current time 
readings, the number recorded in the table is 
obtained by averaging the results of at least 6 
independent measurements of the number of 
coulombs required to displace the boundary 
between two of the division marks on the 


ES 
® International Critical Tables, Vol. 3, p. 90. 





+ 


64 


Fic. 4. {263 


62 








_t l 
0 0-3 





graduated tube; when f was measured by 
coulometer, the entry is the average of three 
independent measurements, each corresponding 
to the f obtained from one coulometer. In 
general, the values of t,o»; thus obtained from 
one experiment are self-consistent to 2 or 3 
units in the fourth decimal place or better, 
although the variation from one experiment to 
the next is somewhat greater. 

The numbers of column 3 are subject to two 
corrections. The first, the solvent correction, 
which must be added to f1.1s, has been shown" 
to be #,.- (conductivity of solvent) /(conductivity 
of solute). The conductivity of the solvent in our 
solutions is from two distinct causes, from the 
excess acetic acid added in making up the 
solution, and from residual impurities (apart 
from carbon dioxide) present in the conductivity 
water. We estimate the specific conductivity of 
the solvent due to this latter factor as 0.6 X 10-°; 
the correction due to the former can be calcu- 
lated from the measured hydrogen ion concen- 
tration of the solution by using the limiting 
mobilities listed by MacInnes, Shedlovsky and 
Longsworth" for the ions involved. The resulting 
solvent corrections for the various solutions are 
given in the fourth column of the table. 

The second, or volume correction is given by 


ty. =tyops-t CAV/1000 


where AV is the change in volume in cc brought 
about at the electrode on the closed side by the 
passage of one faraday. The upper sign is 


13 For a general discussion of the solvent correction, see 
reference 2, and for the special case of an acetate solution 
containing free acetic acid, see reference 10. 

14 MacInnes, Shedlovsky and Longsworth, J. Am. Chem. 
Soc. 54, 2758 (1932). 
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appropriate for cell I and the lower for cell II; 
AV is given by 


AV=Vag— Vageit Vxci—t_* Vac, (I) 


Here ¢, and ¢_ are the transference numbers of 
potassium and acetate ion, Vag, Vagci and Vea 
are the molal volumes of silver, silver chloride 
and cadmium, and Vxci, Vxac and Veaac, are 
the partial molal volumes of potassium chloride, 
potassium acetate and cadmium acetate, respec- 
tively; from the numerical values of these 
quantities listed by Longsworth,' the volume 
corrections of the fifth column are obtained. 
The resulting corrected values of the trans- 
ference number for potassium ion are given in 
the sixth column, and are shown graphically in 
Fig. 3. 

Shedlovsky' has shown that his semi-empirical 
extension to the limiting Onsager equation for 
conductivity fits the experimental data for many 
electrolytes remarkably closely, and Longs- 
worth!® has used an analogous equation to 
represent his transference numbers. For a 1-1 
valent electrolyte, Longsworth gives 


1— 2t..° 


2BV/C 
——ayc+ac(1+ ). (1) 
A’ A’ 





t.=t,°- 


Here C is the concentration, ¢,° is the trans- 
ference number of the cation at infinite dilution, 
A’=An—(aAo+28)\/C, Ao is the equivalent con- 
ductance of the solution at infinite dilution 
(114.37 for potassium acetate at 25°), and a and 
B have their theoretical values for 1-1 electro- 
lytes in aqueous solution at 25°C viz. 0.2274 and 
29.90, respectively; A is a disposable constant. 
If t,.°’ be written for t,°+AC, Eq. (1) is equiva- 
_ lent to 


1° = (tye M+ BYC)/(A'+28VC). (2) 


or 


Hence, a plot of t,° against C should be a 
straight line whose slope is A, and whose inter- 
cept at C=0 is ¢,°. Fig. 4 shows the result of 
such a plot; except for the highest concentration, 


8 Shedlovsky, J. Am. Chem. Soc. 54, 1405 (1932). 
16 Longsworth, J. Am. Chem. Soc. 54, 2741 (1932). 
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the points lie on a straight line for which 
A = —0.105, and #,° is 0.6428, as determined by 
the method of least squares. The last column of 
Table I and the continuous curve of Fig. 3 give 
the values of ¢; computed from Eq. (1) with 
these values of A and ¢t,°; Table II gives values 


TABLE II. 








C 001 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 
t, 0.6498 0.6523 0.6556 0.6579 0.6596 0.6609 0.6620 0.6629 0.6636 








of t, for a few round values of the concentration, 
and the straight line of Fig. 3 shows the limiting 
Debye slope. The discrepancy between the calcu- 
lated and observed values of t, for the highest 
concentration is not surprising, since this con- 
centration lies somewhat above the range where 
the Shedlovsky-Longsworth equation would be 
expected to hold; on the other hand, with such 
a strong solution, a relatively high current must 
be used, and consequently part of this dis- 
crepancy may be simply experimental error due 
to heating of the solution. 

MacInnes, Shedlovsky and Longsworth" give 
as the mobilities of potassium and acetate ion at 
infinite dilution at 25°, 73.50 and 40.87, values 
which they obtained from a consideration of the 
transference and conductance data for a number 
of electrolytes; the corresponding value of /,° is 
0.6427. The fact that the transference number 
obtained here is in such very close agreement is 
further evidence, if any were needed, of the 
essential correctness of their assignment of 
limiting ionic mobilities.!7 In conclusion, we wish 
to express our thanks to Mr. R. W. Allgood for 
assistance in calibrating one of the cells. 

"17 MacInnes, Shedlovsky and Longsworth’s value for 
acetate ion was obtained by subtracting from MaclInnes 
and Shedlovsky’s limiting equivalent conductance for 
sodium acetate 90.97 (see reference 10) the limiting mobil- 
ity of sodium ion, 50.10, which is known from their work to 
a high degree of precision; the extrapolation of the con- 
ductance to infinite dilution was effected by the Onsager- 
Shedlovsky equation. A recalculation of the disposable 


constants by the method of least squares from MaclInnes 
and Shedlovsky’s data shows, however, that the equation 


A=90.935—80.47C!2+ 101.1C(1—0.2274C1?) 


fits their experimental values more closely than does the 
equation they employ. With Ao=90.94, the limiting ionic 
mobility for acetate ion will be 40.84, and the value of t- 
for potassium acetate 0.6428; this change in the transport 
number is of course quite insignificant. 
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The Density of Energy States in Solids 


S. H. BAvER 
Pennsylvania State College, State College, Pennsylvania 


(Received June 4, 1938) 


OR crystals or very complex molecules the 
energy levels are numerous and are either 
inherently or instrumentally unresolvable. As 
yet little progress has been made in devising 
methods of locating these levels although such 
information is of importance to any theory of 
molecular structure. In this note I wish to point 
out that a very general approach is available 
through the use of thermodynamic functions of 
the aggregate; further, to demonstrate that solu- 
tions exist for such equations which are experi- 
mentally indistinguishable from the ones speci- 
fied by the data. 
The sum-over-states may be readily com- 
puted from the entropy, 


1 7 
log Z(T)=— S(T)dT, 
og maf (T) 


while Z(T) = Ddogme—2~!*T (1) 
isa fundamental definition of quantum statistical 
mechanics. The problem is reduced to an inver- 
sion of (1) so that gn(E») could be found once 
a Z(T) is given. Clearly a rigorous solution is 
impossible even if the sum-over-states were 
known for all temperatures since, as T—>~@, 
Z(T)— very rapidly (in the order of T¥). 
Several of a number of possible approximations 
are briefly discussed below. 

A. The function Z(T) is continuous, finite, 
single valued and specified over the region (0, 71). 
However, it may have a discontinuous first de- 
rivative at temperatures at which the substance 
undergoes a change of phase. One may therefore 
approximate it to within any prescribed accuracy 
by a polynomial in T in the above range. 


Z(T)=Ca.T, Z(0)=0. 


The change of variable 7 =1/s results in 


Z(1/s) = Y(s) = Xan(1/s)"** 


so that Y(s) is analytic in the interval (s!, ~). 


Now, Eq. (1) may be rewritten in the form 


V)= f g(B)e*MdE (2) 
0 

where g(£) is the desired density function. The 
change from the sum to the integral is permissible 
because of the closely spaced energy levels. Since 
for all s whose real part is positive 


co 


(1/syri= f eE/k En /(nik"+1)dE, 
0 


Eq. (2) is satisfied by 


g(E)=d0a,E"/(n!k"*1) for s> 54. 


Therefore a solution for (2) exists for all s2 5. 

B. If the range OS s< & is desired, an equa- 
tion approximating (2) must be treated. Intro- 
duce the function 


Y(s+s1) = Da,(1/s+si)"** 


which is analytic for s20. The necessary and 
sufficient conditions that a solution for 


vs+s)= f kg(E)e-*/d(E/k) (3) 
0 


exists are then satisfied (inversion by means of 
the Laplace transform!). It follows that if the 
maximum temperature for which Z(7) is known 
is large—i.e. s; is sufficiently small—so that 
Y(s+s,) may be identified with Y(s) to within 
the accuracy of the data, that solution could be 
accepted also for (2). One should keep in mind 
that often for T>8, a safe extrapolation of S(T) 
to larger T may be made. 
C. Write Eq. (1) in the form 
BkT 
Z(T)= g(E)e-*!*TdE; 

£0 
Eo describes the ground state. The approximation 
enters in replacing the upper infinite limit by 


1D—D. V. Widder, Trans. Am. Math. Soc. 36, 116 (1934). 


403 














































404 


BkT where £ is a very large number, of such mag- 


co 


nitude as to make f g(E)e-#/*7dE negligible. 
BT 
A change of variable BkT=¢ transforms the 


above into a Volterra integral equation of the 
first kind, 


t 
Z(1/6)= | (Be **'aB (4) 
Eo 
with Eyg EX t< BRT,. Integration by parts or a 
differentiation leads to an equation of the second 
kind for which a unique continuous solution 
exists, in the latter case provided Z(E)/8) =0.? 
This will be almost true if 8 is taken sufficiently 


2G. Kowalewski, Integralgleichungen (Gruyter, 1930). 


R. BRDICKA 














equation 


t 


Z(t/8) —Z(Eo/B) = f e(E)e-8'dE, 


Eo 


where Z=Z except in that the discontinuities in 
the derivative of Z if there are any, have been 
smoothed over. 

I have undertaken the application of these 
suggestions to the case of diamond for which 
accurate specific heat data have recently been 
published.’ Pitzer found that the experimental 
curve lies between the one predicted using the 
continuous solid model and that deduced from 
the three-dimensional lattice theory of Blackman. 


3K. S. Pitzer, J. Chem. Phys. 6, 68 (1938). 
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R. BroiéKa 
Physico-Chemical Institute of the Charles University, Praha, Czechoslovakia 





INCE in the recent paper of J. N. Pearce and 
Lyle R. Dawson,! the inadequacy of existing 
explanations of the well-known color changes in 
cobalt chloride solutions has been mentioned, the 
present author wishes to show that his explana- 
tion put forward eight years ago? cannot be 
regarded as an unsatisfactory one for the experi- 
ments given by the above named authors. 

The explanation of the present author is based 
on a polarigraphic study of current-voltage 
curves obtained in electrolysis of cobalt chloride 
solutions with the dropping mercury cathode, as 
well as on spectroscopic research. Through these 
investigations it has been proved that cobalt 
deposits at the dropping mercury cathode from 
blue solutions with a large concentration of 
calcium chloride at a more positive potential 


1J. N. Pearce and Lyle R. Dawson, J. Chem. Phys. 6, 
128 (1938). 

2R. Brdi¢ka, Coll. Czechoslov. Chem. Comm. 2, 489 
and 545 (1930). 





(Received May 16, 1938) 


than it does from the pure pink water solutions. 
Further, while the current-voltage curves in the 
latter case do not conform to the expected equa- 
tion for the deposition of a divalent metal, the 
deposition from the blue cobaltous solutions 
proceeds according to this equation. The devia- 
tion from this equation has been explained as 
being due to the slowness of dehydration of 
cobaltous ions, which takes place at the surface 
of the cathode in order to deposit metallic cobalt. 
The, easier deposition from blue solutions in 
which a blue complex [CoCl, ]%-”” has to be 
assumed, is imagined to be facilitated by the 
ease of splitting of such a complex into free 
cobaltous ions. 

These facts led the present author to the con- 
clusion that the effect of the addition of a large 
excess of calcium chloride in the hydrated pink 
cobaltous ions is twofold: 1. The activity of 
water is lessened and thus the concentration of 
dehydrated cobaltous cations is increased, 2. the 


large; actually, the solution found will be for the 
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NOTE ON THE PAPER BY J. 


dehydrated cations on the other hand associate 
with chloride anions to form complexes like 
[CoCl, ]”-®” according to the equilibrium 


Co’ *+nCl-2[CoCl, J”. 


This can serve as a satisfactory explanation of 
the findings of Pearce and Dawson, according 
to which equimolal high concentrations of added 
chlorides do not produce equal color changes. 
In accordance with the author’s view quoted 
above, the concentration of the complex depends 
primarily on the dehydrating power of the cation 
added (i.e., on its ionic charge and its ionic 
volume) and the complex formation proceeds 
thus parallel with the dehydration of cobaltous 
ions. This would correspond to some extent with 
the explanation given recently by Pearce and 
Dawson who stated: ‘‘As strong hydrating ions 
are added, they attract and orient solvent mole- 
cules. At high concentration of the added salt 
much of the solvent envelope about the cobalt 
ions is removed. Cobalt and chloride ions are 
then drawn more closely together by electro- 
static attraction.” 

No doubt however, can be cast upon the theory 
that complex ions of the type CoCl,-~ are 
formed, since the existence of such particles has 
been sufficiently justified by the electrolytic 
transfer at which cobalt present in blue solutions 
migrates towards the anode.*-® 

As to the mechanism of the light absorption 
of the cobalt chloride solutions in the green and 


ness G. Donnan and H. Bassett, J. Chem. Soc. 81, 939 
). 
*H.G. Denham, Zeits. f. physik. Chemie 65, 641 (1909). 
° J. Groh, Zeits. f. anorg. allgem. Chemie 146, 305 (1925). 
J. Gréh and R. Schmid, Zeits. f, anorg. allgem. Chemie 
163, 321 (1927). 
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red part of the spectrum, no experimental evi- 
dence has been brought forward by the named 
authors which could prove that the absorption 
is due to the electronic transition in the unfilled 
third quantum group of cobalt. In this respect 
more satisfactory measurements have been given 
by W. R. Brode and R. A. Morton.” ® 

The statement of Pearce and Dawson that 
chloride anions distort the electronic system of 
the cobalt ion in such a way that the electrons 
in the third quantum group are capable of ab- 
sorbing smaller quanta of energy, is by no means 
less hypothetical than the original conclusion of 
the present author, who ascribes the color change 
to the bond of chloride anions in the coordinative 
sphere of cobalt by which the outer electronic 
shell of Cl- ions is supposed to be deformed. 

If one takes into account that the two elec- 
trons of such a deformed outer electronic shell 
of chloride anions are really a part of the elec- 
tronic system belonging to cobalt, one could 
equally accept that the absorption is due just 
to this electron pair supplied to the cobalt by 
Cl- ion. 

Finally, the present author wishes to correct 
the following statement used by Pearce and 
Dawson wrongly ascribed to him, viz.: ‘“That 
such color changes could be due to the deforma- 
tion of the chloride ions (cited by R. Brditka) 
by strongly hydrating cations. .. .’’ The de- 
formation of chloride anions was explained in 
the original paper of the present author as being 
due to the bond with the cobalt central atom, 
but never as being effected by strongly hydrating 
cations. 


7 W. R. Brode, Proc. Roy. Soc. A118, 286 (1928). 
8 W. R. Brode and R. A. Morton, Proc. Roy. Soc. A120, 
21 (1928). 
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Note on the Heat Capacities and Energies of SiCl,, TiCl, and SnCl, 


RosBert C. HERMAN 
Palmer Physical Laboratory, Princeton, New Jersey 


(Received June 10, 1938) 


HE interest in thermodynamic functions 

calculated from spectroscopic data is the 
reason for tabulating some information that has 
been available to us for some time. The heat 
capacities and energies of SiCl,, TiCl, and SnCl, 
were calculated by means of well-known for- 
mulae.! The frequencies (liquid) employed were 
taken from Kohlrausch? and are given in Table I 
along with their degeneracies in parentheses. 
Since the molecules investigated are nonpolar 
there should be no great difference in the calcu- 
lated functions if we employ the frequencies 


TABLE I, 








| vi(1) | v2(2) ¥3(3) ¥4(3) 


SiCl, 422 cm | 148 cm | 608 cm | 220 cm 
TiCl, 386 119 491 139 
SnCl, 367 104 401 136 

















TABLE II. Heat capacities in cal./mole deg. 








SiCl TiCla SnCl 





Cvib Cp 


13.19 | 21.14 
13.41 | 21.36 
13.63 | 21.58 
13.83 | 21.78 
14.02 | 21.97 
14.20 | 22.15 
14.37 | 22.32 
14.52 | 22.47 
14.68 | 22.63 
14.81 | 22.76 
14.94 | 22.89 
15.92 | 23.87 
16.48 | 24.43 


Cvib Cp 


14.50 | 22.45 
14.68 | 22.63 
14.85 | 22.80 
15.01 | 22.96 
15.16 | 23.11 
15.29 | 23.24 
15.42 | 23.37 
15.54 | 23.49 
15.65 | 23.60 
15.75 | 23.70 
15.85 | 23.80 
16.55 | 24.50 
16.94 | 24.89 


Cvib Cp 


15.24 | 23.19 
15.39 | 23.34 
15.53 | 23.84 
15.66 | 23.61 
15.78 | 23.73 
15.90 | 23.85 
16.00 | 23.95 
16.09 | 24.04 
16.19 | 24.14 
16.26 | 24.21 
16.34 | 24.29 
16.89 | 24.84 
17.19 | 25.14 























TABLE IJ]. E°—£,° cal./mole. 








SiCls 


3576 
3962 
4557 
5583 
7726 
9936 


TiCls 


4031 
4444 
5075 
6152 
8377 
10650 
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TABLE IV. 








T°C 


90-234 
136-271 
149-273 


Cp cal./mole deg. 


22.5 
24.5 
24.5 





SiCl, 
TiCl, 
SnCl, 








obtained from the liquids rather than the gases. 
In fact, for gaseous SnCl, the frequencies v; and 
v3 are 367 and 400 cm-!? as compared tc 367 
and 401 cm for the liquid. The heat capacity 
in cal./mole ‘deg. for each gas in its standard 
state is tabulated in Table II and in Table III 
we have tabulated the energies in cal./mole. 
Since these molecules have low frequencies and 
anharmonicity corrections are not available, the 
calculation was not carried beyond 300°C. Table 
IV gives the only experimental heat capacity 
data‘ that could be found for these tetrachlorides 
in the gaseous phase and as can be seen the 
calculated values agree fairly well. The ice point 
on the Kelvin scale has been taken as 273.1°K 
and the fundamental constants employed were 
those given by Birge.® 








1For instance, see Fowler, Statistical Mechanics (Mac- 
millan, 1936). 
? Kohlrausch, Der Smekal-Raman-Effekt (Berlin, 1931). 


3 Braune and Engelbrecht, Zeits. f. physik. Chemie B19, 
303 (1932). 

4 Regnault, Mém. de 1’Acad. 26, 1 (1862). 

5 Birge, Phys. Rev. Sup. 1, 1 (1929). 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


Hindered Internal Rotation of Ethane 


We have further improved the ‘‘hot wire’’ gaseous heat 
capacity apparatus of Kistiakowsky and Nazmi,! elimi- 
nating possible errors due to gas adsorption and ‘‘end 
conductance.”’ Measurements on C2Hs and C,D¢, near 
liquid oxygen temperature show that the observed heat 
capacity is independent of the nature of the wire, of the 
gas pressure, and of the temperature difference between 
the wire and the surrounding wall. The vibrational heat 
capacity of C.H¢ at these temperatures is negligible, and 
that of C2D¢ is certainly less than 0.05 cal./mole deg. By 
subtracting the vibrational, over-all rotational, and trans- 
lational heat capacities from our measured values, we find 
for C2Hg the contribution of the hindered internal rotation 
to be 0.59 and 0.67 cal./mole deg. at 94° and 100°K. For 
C:D¢ similarly we find 1.03 and 1.13 cal./mole deg. at 94° 
and 100°K. The values for light ethane are to be compared 
with Hunsmann’s? data: 0.40 at 95° and 0.56 at 101°. 

The theoretical dependence of the hindered internal 
rotational heat capacity upon the assumed magnitude of 
a restricting potential of the form (1/2)V(1—cos 38) is 
shown in Fig. 1 for the temperatures in question. It will 
be seen that a barrier of about 3000 cal./mole is the only 
assumption consistent with the experimental data, pro- 
vided this method of calculation® is not invalidated by 
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Fic. 1. Internal rotational heat capacity offfethane as a _functionTof 
Potential barrier. Filled circles, this work. Open circles, Hunsmann. 
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not later than the 15th of the month preceding that of the 
issue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


quantum effects due to the existence of nuclear spin 
isomers—a point raised by Hunsmann.? In the following 
letter Professor Wilson reports that these quantum effects 
are completely negligible at these temperatures. It follows 
that the low temperature gaseous heat capacity of ethane 
cannot be interpreted in terms of a low potential barrier 
and these quantum effects as suggested by Hunsmann, 
but requires instead a high restricting barrier as indicated 
above. It also follows that the agreement of the entropy 
calculated semi-statistically by Hunsmann (i.e., statis- 
tically for over-all rotation and translation but graphically 
for the internal degrees of freedom) with the third law 
value is not accidental, as he suggests, but is a direct 
consequence of the existence of the high barrier. On the 
other hand the poor agreement of Hunsmann’s calculated 
heat capacity of C2H¢ for a high barrier with the experi- 
mental points between 100° and 200°K is due to an in- 
adequate method of calculation. Using the hindered rotator 
model with V=3000 cal./mole, satisfactory agreement 
between theory and experiment is obtained over the entire 
temperature range for which the experimental data exist, 
provided the uncertain vibrational frequency in C2Hg is 
taken to be between 1000 and 1200 cm™. This frequency is 
not in accordance with Bartholomé and Karweil’s* vibra- 
tional analysis, a point which will be discussed in papers 
soon to appear from this laboratory. 

Thus, in addition to the other thermodynamic® and 
spectroscopic® evidence’ for the existence of a fairly high 
barrier hindering internal rotation in ethane, the low 
temperature gaseous heat capacity leads directly to the 
conclusion that this barrier must be of the order of magni- 
tude of 3000 cal./mole, without any other alternative. 

G. B. KisTIAKOWSKY 
J. R. LACHER 
*FRED STITT 
Mallinckrodt Chemical Laboratory, 
Harvard University, 


Cambridge, Massachusetts, 
June 13, 1938. 


* National Research Fellow in chemistry. 

1G. B. Kistiakowsky and F. Nazmi, J. Chem. Phys. 6, 18 (1938). 

2 W. Hunsmann, Zeits. f. physik. Chemie B39, 23 (1938). 

3K. S. Pitzer, J. Chem. Phys. 5, 469 (1937). 
(1938) Bartholomé and J. Karweil, Zeits. f. physik. Chemie B39, 1 

5 Kistiakowsky, Ruhoff, Smith and Vaughan, J. Am. Chem. Soc. 58, 
137 (1936); J. D. Kemp and K. S. Pitzer, ibid. 59, 276 (1937). 

6 J. B. Howard, J. Chem. Phys. 5, 451 (1937). 

7 This evidence has not been generally regarded as conclusive. See 

ssel, J. Am. Chem. Soc. 59, 2745 (1937). 
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Absence of Spin Effects in the Heat Capacity of Ethane Gas 


Hunsmann! has suggested that the low value of the heat 
capacity of ethane gas at 100°K does not necessarily prove 
the existence of a high potential barrier in this molecule 
because of the possibility that the calculations have been 
incorrectly made, since no consideration was taken of 
the fact that there are several noncombining nuclear spin 
species, analogous to ortho- and parahydrogen. Because of 
the relatively large moment of inertia of ethane, this 
suggestion seemed inherently unreasonable, but on account 
of the importance of this problem, a detailed calculation of 
the partition function and heat capacity of each of the spin 
species of ethane has been made for several different 
barriers and temperatures, in which full account was 
taken of the statistical weights of each individual rotational 
energy level. The method used to find these weights was an 
extension of that described earlier for the general problem.? 
The rotational group for ethane has 18 elements and 9 
irreducible representations. Some of these occur as con- 
jugate pairs, however, so that there are only six possible 
spin species of ethane. At 90°K the over-all rotation of the 
molecule is classically excited with the result that these six 
species can be divided into two sets, the members of each 
set having the same partition functions. The partition 
function for one set involves the state sum over those 
levels of the internal motion which are nondegenerate, 
while the partition function for the other set involves the 
sum over the degenerate levels of the internal torsional 
motion. For very high barriers the levels of the two species 
become identical and for very low barriers the levels go 
over into those for free rotation. In each of these limiting 
cases the two partition functions are negligibly different 
at 90°K. A detailed calculation for a barrier of 300 cal. 
yields the same result. This means that the equilibrium 
proportions of the two species have practically attained 
their high temperature limiting values at 90°K. Further- 
more the heat capacities of the two species are identical 
at this temperature or higher. Consequently, the usual 
method? of calculation which involves the neglect of spin 
and symmetry effects is entirely justified. 

This eliminates one possible explanation of the heat 
capacity of ethane gas at low temperatures and seems to 
leave the existence of a high barrier as the only alternative. 

E. BricHt WILSON, JR. 

Mallinckrodt Chemical Laboratory, 

Harvard University, 


Cambridge, Massachusetts, 
June 13, 1938. 
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Carbon Isotope Effect in the Raman Spectrum of 
Dimethylacetylene 


For dimethylacetylene Glockler and Davis! reported a 
pair of strong Raman frequencies, 2235 and 2313 cm=! 
which were accompanied by two overlapping pairs of weak 
lines, 2201 and 2280 cm™, 2179 and 2258 cm~, all excited 
by Hg 4358A. (Fig. 1.) They divided these lines into two 
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Fic. 1. Raman frequencies of the triple bond 
in dimethylacetylene (cm™). 


groups, one consisting of the higher members of the pairs, 
the other of the lower frequency members. Each group 
was considered to represent a strong vibration accompanied 
by weaker components which were the first members of P 
rotation series. Such a lopsided rotation-vibration band 
had been proposed by Bonino and Cella? to account for a 
similar group of lines in a-pinene. In both substances 
excellent checks on observed separations were obtained by 
assuming that the rotating portions of the molecules were 
simple rotators with one degree of freedom. In dimethyl- 
acetylene it was assumed that the six hydrogen atoms were 
rotating about the figure axis. 

Of course there were objections to the above interpreta- 
tion. The restricted nature of the rotator calculations and 
the absence of R branches caused doubt. Recently in 
studies conducted here on other low molecular weight 
hydrocarbons? measurements were made on the weak 
Raman shifts due to molecules containing the C! isotope. 
In continuing this study reference was made to the avail- 
able plates on dimethylacetylene in a search for lines of 
isotopic origin. The first pair of weak lines listed, giving a 
separation of 33 cm~! from the corresponding members of 
the strong pair, were considered to be due to molecules 
containing C®, 

The elimination of the first component in each rotation 
series necessitated another explanation for the other 
members. It was then found that the 2179 and 2258 cm™ 
weak lines were exactly in position to be the 2235 and 
2313 cm~ strong shifts excited by the 4347.5A member of 
the mercury triplet. In the paper mentioned no further 
weak lines were reported; however, in an enlargement of 
one of the plates by Davis a faint line was discernible which 
was in position to be the 2235 cm~ shift excited by the 
4339.2A member of the triplet; the 2313 cm~ shift thus 
excited would not be visible because of coincidence with a 
strong line. In consequence of these measurements the 
rotational interpretation of the faint lines is to be dis- 
carded. 

Subsequent to the studies of Glockler and Davis the 
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pair of strong lines, 2235 and 2313 cm™ were explained by 
Badger‘ as manifestations of a single fundamental vibration 
interacting with an overtone of a lower fundamental in 
Fermi resonance due to accidental degeneracy. As a result, 
this carbon-carbon triple bond vibration frequency of the 
normal molecule appears split into two components, the 
actual fundamental frequency lying near the mean value 
of the pair. Similarly the weak pair of lines, 2201 and 
2280 cm=, are now considered to result from the inter- 
action of the carbon-carbon triple bond fundamental of 
isotopic molecules with an overtone of a lower fundamental 
of these isotopic molecules. 

By considering the normal molecule to be a simple 
harmonic two-body oscillator with reduced mass of 
1/u=1/25+1/25 and the isotopic molecule to have similar 
properties with a reduced mass of 1/u=1/25+1/26 it is 
possible to calculate a separation of 20 cm™ for the carbon- 
carbon triple-bond frequencies of the normal and the 
isotopic molecules. This is not checked very well by the 
experimental measurement of 33 cm™ separating the mean 
values of the resonance split pairs. Of course this calcula- 
tion is only an approximation. However, it must be 
emphasized that a close check is not to be expected since 
the separation of mean values of split frequencies does not 
give an accurate measure on separations of the ‘‘unsplit”’ 
carbon-carbon triple bond frequencies. The resonance split 
values observed depend not only on the ‘‘unsplit’’ funda- 
mentals but also on the overtone of the lower fundamentals 
likewise engaged in the resonance interaction.® Because the 
isotopic shift for the lower frequency will differ from that 
for the higher frequency, the separation of mean values of 
split frequencies necessarily will vary from the separation 
of the true, ‘‘unsplit’”’ fundamental vibrations in normal 
and isotopic molecules. 

In the present approximate treatment no consideration 
was given to location of the C!. In dimethylacetylene 
roughly four percent of the molecules are isotopic with 
half of these having isotopic triple bond carbons, the other 
half possessing C!* in the methyl groups. 

GEORGE GLOCKLER 


Institute of Technology, 
University of Minnesota, M. M. RENFREW® 
Minneapolis, Minnesota, 


June 13, 1938. 
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Raman Effect in Liquid Ethylene and Ethane 


In the Raman spectrum of liquid ethylene two lines, 943 
and 3075 cm! have been observed! which have not been 
found in the gas on very long exposures at high pressure.” 
Since the work of Bonner*® on Raman effect in the liquid 
and infra-red absorption in the gas, it has been customary 
in theoretical treatments‘ to use these lines (950 and 3069 
cm™) as fundamental vibrations along with other fre- 
quencies including 940, 950, and 3107 cm~ which are 
infra-red active. According to the selection rules derived 
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for gases,® in molecules having a center of symmetry fre- 
quencies which are infra-red active should not appear in 
Raman effect. Of course there exists the possibility of an 
accidental coincidence of values. However, it is our belief 
that the appearance of the two lines mentioned in the 
Raman spectrum of liquid ethylene constitutes a violation 
of the selection rules since it seems likely that similar 
values in the infra-red absorption spectrum arise from 
the same vibrations. 

The 3075 cm™ line is somewhat diffuse but is of such a 
strength that it certainly should be visible on reasonable 
exposures in the gas, if permitted there. The line cannot 
be interpreted readily as a combination or overtone. Hence 
it seems convenient to relate this line to the strong 3107 
cm! vibration observed in the infra-red. Of course, some 
increase in value is expected in going from liquid to gas 
for high frequencies. The increase here seems large, but 
does not appear excessive. 

The 943 cm™ shift in Raman effect for the liquid is very 
diffuse, and it seems natural to identify this with the 940 
and 950 cm™ lines in the infra-red. For low frequencies 
no great drop in value is observed generally in passing 
from gas to liquid. 

Moreover, a similar comparison of data is possible with 
the Raman spectrum of liquid ethane and its infra-red 
absorption spectrum for the gas. In pure liquid ethane the 
Raman line 1463 cm™, though diffuse, is well established.® 
However, in the gas Lewis and Houston’ could not find it 
and aithough Bhagavantam® measured a shift in this 
neighborhood he was using an impure sample and con- 
sidered the line observed to be due to a contaminant. 
Now in infra-red absorption experiments® strong bands 
are observed at 1465 and 1495 cm™. Perhaps the mean 
value of 1480 cm™ for the gas corresponds to the 1463 cm™ 
Raman shift observed in the liquid. 

In addition there are indications of a faint, diffuse 
Raman shift of 813 cm~ in liquid ethane which may be 
related to the infra-red band 827 cm™. 

Now it is possible that the failure to observe these 
Raman iines in the gas is due to their low intensities. 
But other lines of comparable strength (frequently of 
lower intensity) have been recorded in the gas. For this 
reason it seems likely that the Raman lines discussed are 
permitted only in the liquid state. The shifts are of nearly 
the same magnitude as infra-red active frequencies (see 
Table I), which fact may not be accidental. The presence 
of the C'* isotope has been shown to be responsible for 
slightly altered Raman shifts of low intensity accom- 
panying corresponding strong lines excited in the normal 
molecules. The presence of this isotope contributes a small 
degree of asymmetry to the molecules which would cause 
infra-red active frequencies to be permitted in the Raman 
effect. However, such manifestations should be visible on 
long exposures in the gas as well as in the liquid. Since this 
condition has not been observed, it seems more reasonable 
to suppose that if the Raman and infra-red frequencies do 
belong to identical fundamentals, then a breakdown in 
selection rules in the liquid state for certain lines is indi- 
cated and must be due to the more strenuous environment 
existing there. Thompson and Linnett* have offered this 


























LETTERS TO THE EDITOR 














RAMAN EFFECT 
INFRA-RED 
SUBSTANCE (GAS) (LIQUID) (GAS) 
Ethylene 940 
>945 943 not found 
0 
3107 3075 not found 
Ethane 827 813 not found 
465 
> 1480 1463 ? (8) 
1495 

















suggestion in the case of ethylene, but for purposes of 
comparing calculated values with experimental data they 
considered the Raman lines in question to be bona fide 
Raman active fundamentals. Our experience adds weight 
to their tentative objections for the ethylene lines. More- 
over, data given for ethane are substantiating. Perhaps 
theoretical treatments of the ethylene molecule will profit 
by including among the Raman active fundamentals one 
or both of the high frequency values, 3230 and 3263 cm“, 
whose existence we have verified.! 

Recently we have made more accurate evaluations on 
our plates for the Raman spectrum of liquid ethylene, 




















resulting in slight alteration of previously published 
values.! The new measurements are: 943 (1d), 1341 (10), 
1602 (0.5), 1621 (7), 1653 (1), 2870 (2), 2997 (0.5), 3008 
(10), 3075 (2d), 3230 (0.5), 3263 (0.5). The estimated in- 
tensities are in parentheses. 

The 1602 cm™ shift is the carbon-carbon vibration in 
the isotopic molecule H,C%=C!H,. The shift 2997 cm™ 
probably also is due to this isotopic molecule, the separa- 
tion from the 3008 cm™ vibration of the normal molecule 
being checked roughly by calculation. However:, the accu- 
racy of measurement has been diminished by the lack of 
clear resolution of the weak and strong line. Our value for 
the isotope line probably is low. 

GEORGE GLOCKLER 

Institute of Technology, M. M. RENFREW"® 
University of Minnesota, 


Minneapolis, Minnesota, 
June 13, 1938. 
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